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(2)  determination  of  die  internal  triaxial  stresses  due  to  the  surface  pressure; 
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failure  modes.  In  studying  the  impact  response  of  composite  cylinders,  a 
computer  program  capable  of  handling  stress  wave  effects  was  used.  The 
analytical  techniques  that  were  developed  were  used  to  study  the  influence  of 
the  following  parameters  on  impact  response;  fiber  and  matrix  properties, 
fiber  layup  and  stacking  sequence,  thickness  effects,  curvature  effects  and 
hybridization..  At  the  end  of  Phase  II  a test  plan  was  prepared.  In  Phase  III 
the  teat  plan  was  implemented,  which  involved  fabrication  and  impact  testing 
of  specimens.  Types  of  graphite  fibers  used  in  the  experimental  program 
included  Thornel  300,  Modmor  II  and  Celion  GY-70.  Types  of  resins  investi- 
gated included  5208  epoxy,  ERLA  4617  and  polysulfone.  Fiber  layups  investi- 
gated included  unidirectional,  2:1  bidirectional,  1:1  bidirectional  and 
tridirectional  (pseudo-isotropic).  In  the  case  of  tridirectional  composites, 
three  different  stacking  sequences  were  used  to  establish  how  die  latter  influ- 
ences the  impact  response.  Two  types  of  hybrid  composites  that  were  investi- 
gated included  Thornel  300/Fiberglass/S208  and  Thornel  300/Kevlar  49/5208. 

The  investigation  of  the  influence  of  specimen  geometry  included  fabrication 
and  impact  testing  of  composite  plates  of  three  different  thicknesses  (0.  066  in. , 

0.- 139  in.,  and  0.  260  in. ) and  fabrication  and  testing  of  3-in.  - and  6-in.  -diameter 
bidirectional  Thornel  300/5208  composite  cylinders.'  The  important  results 
if  the  theoretical  and  experimental  studies  were  as  follows:  (1)  Resistance  to 
impact  damage  increases  as  the  strength  of  the  graphite  fibers  increases  and 
the  modulus  of  the  fibers  decreases;  Of  die  three  fibers  investigated,  composites 
made  with  Thornel  300  showed  best  resistance  to  impact  damage.  (2)  Resistance 
to  imp ict  damage  increases  as  the  Young's  modulus  of  the  matrix  decreases 
and  the  strength  of  the  matrix  increases.  The  energy  level  required  to  cause 
visible  damage  in  composites  made  with  polysulfone  matrix  and  Thornel  300 
fibers  was  almost  an  order  of  magnitude  higher  than  that  for  composites  made 
with  other  resins;  the  mechanical  properties  of  the  Thornel  300/polysulfone 
matrix  were  poor,  however.,  (3)  Bidirectional  layup  was  more  efficient  in 
resisting  impact  damage  than  tridirectional  or  unidirectional  layup.  (4)  Con- 
struction employing  complete  dispersion  of  layers  (having  different  fiber 
orientation)  through  the  thickness  was  more  resistant  to  impact  damage  than 
construction  m which  the  layers  were  not  dispersed.  (5)  Effects  of  Hybridization 
on  impact  damage  were  inconclusive  because  of  poor-quality  test  specimens. 

(6)  Impact  of  a spherical  impactor  into  cylindrical  target  or  a flat  target  having 
unidirectional  fiber  layup  resulted  i>.  an  elliptical  area  of  contact;  whereas  in 
plates  having  1:1  bidirectional  or  tridirectional  fiber  layup,  the  area  of  contact 
was  circular.  (7)  The  visible  damage  resulting  from  impact  consisted  of  trans- 
verse cracks,  running  parallel  to  the  fibers,  on  the  back  face  of  the  plates  or 
cylinders.  (8)  Theoretical  prediction!,  of  impact  damage  in  composites  showed 
fair  correlation  with  test  results;  however,  further  refinement  of  the  theory 
appears  desirable. 
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INTRODUCTION 


During  the  past  decade,  significant  progress  has  been  made  in 
developing  advanced  composites  possessing  high  strength,  high  mod- 
ulus, and  low  density,  and  in  understanding  their  behavior  under 
certain  types  of  loading.  Asa  result,  composites  such  as  boron- 
epoxy  and  graphite  epoxy  have  been  successfully  employed  as  struc- 
tural materials  in  aircraft,  missiles,  and  space  vehicles,  and  have 
satisfactorily  demonstrated  their  performance  through  extensive 
ground  testing  and  in  flight.  Moreover,  the  application  of  these  com- 
posite materials  to  various  structural  components  of  aerospace  vehicles 
has  demonstrated  significant  (10  to  ~50%)  weight  savings  over 
comparable  components  made  of  conventional  metals. 

Despite  the  tremendous  advantages  that  advanced  composites  have  over 
metals  in  applications  requiring  high  strength,  high  stiffness,  and  low 
weight,  in  applications  where  impact  by  foreign  objects  is  a design  con- 
sideration, the  advantages  inherent  in  composites  are  overshadowed  by 
their  poor  resistance  to  impact  loading.  The  severity  of  this  problem 
\ | was  demonstrated  when  a major  aerospace  company  attempted  to  use 

| graphite-epoxy  turbine  blades  in  one  of  its  production  engines.  During 

| a test  flight  through  a rain  storm,  the  turbine  blades  were  severely 

] damaged  and  eroded  and,  although  the  plane  landed  safely,  the  near- 

j disaster  created  an  impetus  for  an  in-depth  understanding  of  the  prob- 

| lem  and  a search  for  means  to  overcome  it.  Numerous  other  reports 

! in  the  open  literature  dealing  with  the  response  of  advanced  composites 

I to  various  types  of  impact  have  further  increased  the  need  for  a better 

understanding  of  the  problem  so  that  the  survivability  of  composites 
! under  various  types  of  impact  loading  can  be  increased. 

! The  studies  presented  in  this  report  are  oriented  toward  this  need;  that 

is,  toward  understanding  brittle  fracture  in  graphite -epoxy  composites 
subjected  to  foreign-object  impact.  The  approach  used  in  studying  the 
problem  consists  of  experimental  and  theoretical  efforts  and  is  dis- 
cussed in  detail  in  the  sections  that  follow. 
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SECTION  I 

LITERATURE  SURVEY  RELATIVE  TO  IMPACT  BEHAVIOR 


AND  BRITTLE  FRACTURE  IN  GRAPHITE -EPOXY  COMPOSITES 


A literature  survey  was  conducted  on  brittle  fracture  in  graphite - 
epoxy  composites  subjected  to  impact  loading,  and  a bibliography  was 
compiled  on  the  subject.  In  order  for  the  results  obtained  from 
the  literature  survey  to  be  of  maximum  benefit  to  this  program,  papers 
and  publications  dealing  with  the  following  areas  were  sought,  reviewed, 
and  evaluated: 

1.  Brittle  fracture  of  composites  under  static  loading,  including 
strength  prediction  from  given  loading. 

2.  Behavior  of  composites  under  impact  loading  (experimental 
and  theoretical  work). 


3.  Brittle  fracture  of  composites  under  impact  loading  (experi- 
mental and  theoretical  work). 

4.  Test  methods,  specimens,  and  condition,  for  assessment  of 
impact  response  of  brittle  composites. 

5.  Material  design  concepts  for  improving  impact  resistance  of 
composites  (hybrid  composites,  protective  coatings,  etc. ). 

6.  Other  pertinent  work  on  the  subject  such  as  stress -wave 
response  of  composites  and  failure  in  brittle  materials 
under  impact. 

The  results  of  the  literature  survey  are  presented  under  three  main 
subheadings:  (1)  Behavior  of  Composites  Under  Impact  Loading, 

(2)  Failure  Criteria  for  Brittle,  Orthotropic,  Anisotropic,  and  Composite 
Materials,  and  (3)  Stress- Wave  Response  of  Composite  Materials, 

BEHAVIOR  OF  COMPOSITES  UNDER  IMPACT  LOADING 


Material  dealing  with  the  behavior  of  composites  subjected  to  impact 
loading  has  only  recently  begun  appearing  in  the  literature.  The  amount 
of  literature  in  this  area  is  still  quite  limited.  Brief  summaries  of  the 
more  pertinent  papers  related  to  impact  response  of  composites  and 
anisotropic  materials  are  presented  below. 
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A number  of  papers  on  the  behavior  of  composites,  including  graphite- 
epoxy  composites,  subjected  to  various  types  of  impact  loading  were  pre- 
sented at  a recent  ASTM  Symposium  on  Foreign  Object  Behavior  of  Com- 
posites, held  in  Philadelphia,  Pennsylvania,  in  September  1973.  The 
papers  that  were  presented  will  appear  as  an  ASTM  Special  Technical 
Publication  (STP),  which  will  also  include  some  papers  that  were  not 
presented  at  the  symposium.* 

The  symposium  focused  on  the  understanding  of  the  foreign-body  impact 
failure  mechanism  of  filamentary  composites,  identification  of  suitable 
tests  for  impact  response  studies,  experimental  studies  on  performance 
of  composite  laminates  and  sandwich  construction  under  impact  condi- 
tions, and  the  selection  and  evaluation  of  suitable  protective  coating 
materials  or  design  techniques  for  increasing  the  survivability  of  com- 
posites in  impact-type  environments. 

Several  of  the  papers  that  contained  data  and  results  directly  applicable 
to  this  program  are  briefly  described  below. 

In  a paper  by  Oplinger  and  Slepetz^,  test  results  are  presented  on  the 
behavior  of  sandwich  panels  with  faces  made  of  different  composites  and 
subjected  to  falling -weight  impact  tests.  The  variables  in  this  study 
were  different  types  of  composite  materials,  fiber  orientation  and  layup, 
sandwich  core  density,  and  loading  parameters.  Strength  measurements 
were  made  on  the  initial  test  specimens  and  also  on  specimens  that  were 
subjected  to  impact  loading.  Sandwich  panels  made  with  graphite -epoxy 
faces  were  shown  to  be  more  susceptible  to  impact  damage  than  similar 
panels  made  with  S-glass-epoxy  faces.  The  poor  impact  strength  of 
graphite -epoxy  composites  was  attributed  to  their  low  strains  to  failure. 

The  investigation  described  by  Broutmen,  Rotem,  and  Zych^  involved 
impact  testing  of  various  composites  using  an  instrumented  drop-weight 
impact  machine  and  a more  conventional  instrumented  Charpy  impact 
machine.  Various  failure  modes  were  discussed,  as  well  as  the  influ- 
ence of  composite  materials  and  lamination  on  the  failure  modes  caused 
by  impact.  Instrumented  Charpy  testing  was  also  used  in  the  studies 


♦References  1 through  12,  listed  on  page  201,  are  the  papers  that  were 
presented  at  the  ASTM  Symposium.  References  13  and  14  are  the  papers 
that  were  not  presented.  Individual  references  will  be  cited  throughout 
where  discussed. 


Oplinger,  D.  W.,  and  Slepetz,  J.  M. , IMPACT  DAMAGE  TOLERANCE 
OF  GRAPHITE/EPOXY  SANDWICH  PANELS,  U.S.  Army  Materials  and 
Mechanics  Research  Center,  Watertown,  Mass.,  September  1973. 

Broutmen,  L.  J.  Rotem,  A.,  and  Zych,  J. , IMPACT  ENERGY  OF 
GLASS  FIBER  AND  GLASS -CARBON  FIBER  HYBRID  COMPOSITES, 
Illinois  Institute  of  Technology,  Chicago,  Illinois,  September  1973.’ 
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reported  by  Beaumont,  Riewald,  and  Zweben?.  The  instrumentation  was 
such  that  both  the  impact  loads  and  the  impact  energy  histories  were 
measured.  Types  of  materials  investigated  included  unidirectional 
E- glass -epoxy,  PRD-49 -epoxy,  and  graphite- epoxy  composites.  Impact 
behavior  of  hybrid  composites  consisting  of  high-modulus  graphite  and 
PR.D-49  fibers  in  epoxy  matrix  was  also  investigated.  Hybridization  was 
shown  to  increase  the  "impact  strength"  by  a factor  of  10  in  comparison 
to  that  of  all- graphite  composites. 

The  paper  by  Greszczuk***  presents  theoretical  studies  on  the  response 
of  isotropic  and  composite  materials  to  particle  impact.  Solutions  are 
given  therein  for  the  magnitude  and  distribution  of  the  time-dependent 
pressure  resulting  from  particle  impact,  intern  *1  triaxial  stresses  in 
the  target  caused  by  the  surface  pressure,  and  the  target  failure  modes 
caused  by  the  internal  stresses.  Types  of  composite  materials  for 
which  results  are  presented  include  glass-,  boron-,  and  graphite- epoxy. 
It  is  shown  that  the  impact  resistance  of  materials  is  governed  by  elastic 
and  strength  properties  of  the  target  materials,  and  that  the  sequence  of 
failure  modes  changes  as  the  mechanical  properties  of  the  composite  are 
changed. 

In  addition  to  papers  described  above,  a number  of  other  papers  and 
publications  were  found  concerning  the  behavior  of  composites  under 
impact  loading. 

A paper  by  Chamis,  Hanson,  and  Serafim* 5 contains  theoretical  and 
experir  iental  studies  on  the  impact  resistance  of  unidirectional  com- 
posites made  with  S-glass,  PRD-49,  and  various  types  of  graphite  fibers, 
as  well  as  on  the  impact  resistance  of  hybrid  composites.  Miniature  Izod 
and  longitudinal  impact  tests  were  used  in  the  experimental  portion  of 
this  investigation  to  evaluate  the  impact  resistance  of  unidirectional 
composites  in  the  longitudinal  and  transverse  directions.  Impact 
resistance  for  the  shear  loading  is  also  investigated  experimentally. 
Experimental  and  theoretical  results  are  also  presented  on  failure 
modes  associated  with  impact  loading  of  specimens  described  above. 
Types  of  failures  observed  in  various  types  of  composites  include 


^Beaumont,  P.  W. , Riewald,  P.  G. , and  Zweben,  C. , THE  IMPACT 
FRACTURE  CHARACTERISTICS  OF  SOME  FIBER-STRENGTHENED 
EPOXY  RESIN  SYSTEMS,  E.  I.  duPont  de  Nemours  and  Company, 
Incorporated,  Wilmington,  Delaware,  September  1973. 

10Greszczuk,  L.  B. , RESPONSE  OF  ISOTROPIC  AND  COMPOSITE 
MATERIALS  TO  PARTICLE  IMPACT,  McDonnell  Douglas  Astronautics 
Company,  Huntington  Beach,  California,  September  1973. 

15Chamis,  C.  C. , Hanson,  M.  P. , and  Serafini,  T.  T. , IMPACT 
RESISTANCE  OF  UNIDIRECTIONAL  FIBER  COMPOSITES,  IN  COM- 
POSITE MATERIALS:  TESTING  AND  DESIGN,  American  Society  for 
Testing  and  Materials,  ASTM  STP  497,  February  1972,  pp.  324-349. 
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cleavage,  fiber  debonding  and  pullout,  delaminations,  matrix  fracture, 
and  fiber  splitting.  Impact  strength  in  the  transverse  direction  is 
shown  to  increase  with  increasing  fiber  tensile  strength  and  decreasing 
fiber  modulus. 

A paper  by  Novak  and  DeCrescente  presents  Charpy  impact  strength 
for  unidirectional  giaphite-,  boron-,  and  glass-resin  composites  tested 
in  the  fiber  direction.  The  impact  strength  of  glass-resin  is  shown  to 
be  significantly  higher  than  that  of  either  the  boron  or  the  graphite  com- 
posites. The  mechanisms  by  which  resin-matrix  composites  absorb 
energy  are  studied,  and  it  is  found  that  the  tensile  stress- strain  charac- 
teristics of  the  fibers  are  of  primary  importance  in  determining  the  level 
of  composite  impact  resistance.  The  toughness  of  the  resin  matrix  is 
shown  to  represent  a small  contribution  to  composite  impact  energy. 
Increasing  the  impact  strength  of  boron  and  graphite  composites  by  the 
inclusion  of  glass  fibers  is  discussed. 

A recent  report  by  Friedrich  and  Preston*?  presents  studies  on  impact 
resistance  of  composite  blades  used  in  aircraft  turbine  engines.  Impact 
damage  caused  by  steel,  gelatin,  and  ice  projectiles  was  investigated 
experimentally.  Types  of  composites  investigated  included  S-glass- 
epoxy,  boron- epoxy,  PRD-49-epoxy,  various  types  of  graphite- epoxy, 
and  hybrid  composites.  Impact  tests  were  performed  on  cantilevered 
flat  plates  and  simulated  fan-blade  specimens  consisting  of  cantilevered 
double- tapered  specimens.  Other  variables  investigated  included  mass 
and  size  of  the  projectile,  impact  velocity,  impact  angle,  and  material 
layup.  Photomicrographs  showing  internal  damage  caused  by  impact 
onto  various  composite  materials  are  also  presented.  The  test  results 
indicate  that  the  most  important  factor  in  improving  the  impact  resistance 
is  to  increase  the  transverse  tensile  strength  of  composites.  A paper 
based  on  the  above- described  work  was  presented  at  the  A5TM  Symposium 
by  Preston  and  Cook’2, 


l6Novak,  R.  C.  , and  DeCrescente,  M.  A.,  IMPACT  BEHAVIOR  OF 
UNIDIRECTIONAL  RESIN  MATRIX  COMPOSITES  TESTED  IN  FIBER 
DIRECTION,  IN  COMPOSITE  MATERIALS:  TESTIN'".  AND  DESIGN, 
American  Society  for  Testing  and  Materials,  A3TM  SIP  497, 
February  1972,  pp.  311-323. 

*?Friedrich,  A.  L.  , and  Preston,  J.  L.  , Jr.,  IMPACT  RESISTANCE 
OF  FIBER  COMPOSITE  BLADES  USED  IN  AIRCRAFT  TURBINE 
ENGINES,  Pratt  and  Whitney  Aircraft,  Division  of  United  Aircraft 
Corporation;  NASA  Technical  Report  CR- 134502,  National  Aero- 
nautics and  Space  Administration,  NASA  Lewis  Research  Center, 
Cleveland,  Ohio,  May  1973. 

12Preston,  J.  L. , Jr.,  and  Cook,  T.  S. , IMPACT  RESPONSE  OF 
GRAPHITE-EPOXY  COMPOSITES  TO  SIMULATED  AIRCRAFT 
ENGINE  FOREIGN  OBJECT  DAMAGE  ENVIRONMENTS,  Pratt  and 
Whitney  Aircraft  Corporation,  Division  of  United  Aircraft,  East 
Hartford,  Conn.  , Sept.  1973. 
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Another  publication  in  the  area  of  the  response  of  graphite-epoxy 
turbofan  blades  subjected  to  impact  is  a report  by  Coppal8.  Coppa 
describes  experimental  studies  on  the  effects  of  simulated  bird-carcass 
impact  on  graphite- epoxy  turbofan  blades,  including  test  facilities  and 
instrumentation  that  were  utilized  and  the  results  that  were  obtained. 

Such  quantities  as  energy  and  momentum  transferred  to  the  blades, 
blade  impact  forces,  vibration  frequencies,  damping  factors,  dynamic 
tip  deflections,  and  bending  moments  were  evaluated. 

A paper  by  Sayers  and  Harris^  is  another  recent  publication  related  to 
application  of  graphite -epoxy  composites  to  turbofan  engine  blades.  The 
paper  presents  test  results  on  the  interlaminar  shear  strength  of 
graphite -epoxy  composites  subjected  to  impact  loading.  Both  static 
and  drop-ball-impact  interlaminar  shear  tests  were  conducted  on  uni- 
directional and  cross -plied  graphite-epoxy  composites.  The  interlaminar 
shear  strength  from  the  impact  tests  was  found  to  be  about  70  percent  of 
the  static  interlaminar  shear  strength.  The  static  and  impact  tests 
lasted  approximately  50  and  0.  0007  seconds,  respectively. 

A paper  by  Bradshaw,  Dorey,  and  Sidey^®  describes  the  various  types 
of  fracture  observed  in  graphite -epoxy  composites  subjected  to  impact. 
These  are  discussed  in  terms  of  the  fracture  energies  and  the  stresses 
to  cause  fracture.  Experiments  are  described  showing  how  component 
geometry  and  material  variables  affect  both  the  elastic  3train  energy 
before  fracture  and  the  type  and  extent  of  impact  fracture  damage. 

Failure  and  design  criteria  for  different  applications  are  discussed. 
Suggested  impact  tests  include  threshold  drop-weight  or  ballistic  tests 
for  determining  the  resistance  against  handling  or  repeated  low-energy 
impact,and  Izod  or  Charpy  tests  for  containment.  An  assessment  is 
made  of  the  various  material  modifications  for  optimizing  impact  per- 
formance. These  include  fiber,  matrix,  interface,  mixed  carbon- 
glass  fiber  j,  and  three-dimensional  reinforcement. 


18Coppa,  A.  P.  , MEASUREMENT  OF  TURBOFAN  BLADE  RESPONSE 
TO  SIMULATED  BIRD-CARCASS  IMPACT,  General  Electric  Report, 

71  SD  266,  December  1971. 

19Sayers,  K.  H. , and  Harris,  B.  , INTERLAMINAR  SHEAR  STRENGTH 
OF  A CARBON  FIBER  REINFORCED  COMPOSITE  MATERIAL  UNDER 
IMPACT  CONDITIONS,  J.  Composite  Materials,  Vol.  7,  January  1973, 
pp.  129-133. 


^Bradshaw,  F.  J.  , Dorey,  G. , and  Sidey,  G.  R.  , IMPACT  RESIST- 
ANCE OF  CARBON  FIBER  REINFORCED  PLASTICS,  Royal  Aircraft 
Establishment  Technical  Report  72240,  March  1973. 
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A more  recent  paper  by  Dorey2*  discusses  fracture  processes  by  which 
graphite -epoxy  composites  fail  under  impact  loading  and  also  gives  the 
residual  strengths  and  stiffnesses  of  several  composites  that  have  been 
subjected  to  drop-weight  and  ball -gun  impact.  Material  modification 
aimed  at  minimizing  certain  types  of  impact  damage  are  also  described 
therein. 


A recent  report  by  Moon22  presents  calculations  for  stresses  and 
displacements  induced  in  anisotropic  plates  by  short- duration  impact 
forces.  The  theoretical  model  attempts  to  model  the  response  of  fiber 
composite  turbine  fan  blades  to  impact  by  foreign  objects  such  as  stones 
and  hailstones.  In  this  model,  Hertzian  impact  theory  is  used  in 
determining  the  impact  force.  The  plate  response  treats  the  laminated 
blade  as  an  equivalent  anisotropic  material,  using  a form  of  Mindlin's 
theory  for  crystal  plates.  The  analysis  makes  use  of  a computational 
tool  called  the  "fast  Fourier  transform.  " Results  are  presented  in  the 
form  of  stress  contour  plots  in  the  plane  of  the  plate  for  various  times 
after  impact.  Examination  of  the  maximum  stresses  due  to  impact 
versus  ply  layup  angle  reveals  that  the  ±15  layup  angle  gives  lower 
flexural  stresses  than  0 , ± 30  , and  ± 45  cases. 

In  another  publication.  Moon  presents  a survey  of  papers  on  wave  propa- 
gation and  impact  response  in  composite  materials.^  The  major 
portion  of  the  references  presented  therein  are  on  various  wave  propa- 
gation problems  in  composite  materials. 


21Durey,  G.  , FRACTURE  BEHAVIOR  AND  RESIDUAL  STRENGTH  OF 
CARBON  FIBER  COMPOSITES  SUBJECTED  TO  IMPACT  LOADS, 
Paper  B3,  presented  at  the  NATO  AGARD  Meeting,  Munich, 

Germany,  October  7-9,  1974., 

22Moon,  F.  C.  , THEORETICAL  ANALYSIS  OF  IMPACT  IN  COMPOSITE 
PLATES,  National  Aeronautics  and  Space  Administration  Report  NASA 
CR-121110,  1972. 

23Moon,  F.  C.,  A CRITICAL  SURVEY  OF  WAVE  PROPAGATION  AND 
IMPACT  IN  COMPOSITE  MATERIALS,  National  Aeronautics  and 
Space  Administration  Report  NASA  CR-121226,  May  1973. 
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Studies  of  Hertzian  impact  in  anisotropic  bodies  have  been  made  by- 
Chen24  and  by  Willis2*.  Willis  shows  that  the  area  of  contact  in 
anisotropic  materials  is  elliptic,  in  contrast  to  its  being  a circle  in 
isotropic  materials.  Earlier  work  (1964)  by  Svelko26  considered  the 
Hertzian  contact  problem  in  transversely  isotropic  materials.  The 
solutions  to  the  contact  problem  can  also  be  obtained  from  the  general 
theories  presented  by  Lekhnitskii27  and  by  Green  and  Zerna28.  It  is 
noted  here  that  by  combining  the  solution  for  the  contact  problem  with 
the  solution  for  impact  between  two  bodies,  the  solution  for  stresses  in 
a target  subjected  to  impact  loading  can  be  obtained^. 


Other  papers  dealing  with  behavior  of  composites  subjected  to  impact 
loading  are  given  as  References  24-51^  and  will  be  cited  where  discussed 
in  this  report. 

FAILURE  CRITERIA  FOR  BRITTLE,  ORTHOTROPIC.  ANISOTROPIC, 
AND  COMPOSITE  MATERIALS" 


To  understand  the  brittle  failure  of  graphite -epoxy  composites  requires 
an  understanding  of  how  such  materials  fail  under  multiaxial  and 
combined  loading.  The  failure  criterion  used  in  predicting  impact- 
induced  failures  is  therefore  one  of  the  key  factors  in  understanding 
the  material  response.  For  this  reason,  and  as  a part  of  the  literature 
survey,  a bibliography  has  been  compiled  on  the  failure  theories  for 
brittle,  orthotropic,  anisotropic,  and  composite  materials.  The 
various  publications  dealing  with  failure  of  such  materials  are  given 
in  alphabetical  order.  Most  of  the  papers  are  on  failure  of  anisotr  >pic 
and  composite  materials.  Several  pertinent  references  are  also 
included  on  the  failure  of  brittle,  isotropic  materials  under  static 
and  impact  loading. 

Papers  by  Franklin,  Kaminski  and  Lantz,  Sandhu,  and  Sendeckyj  con- 
tain reviews  and  evaluations  of  a number  of  failure  theories  that  are 
given  in  the  literature.  Many  of  the  failure  theories  given  in  the 
bibliography  are  reviewed  in  a paper  by  Sandhu.  Of  the  various  failure 


24Chen,  W.  T. , STRESSES  IN  SOME  ANISOTROPIC  MATERIALS  DUE 
TO  INDENTATION  AND  SLIDING,  International  Journal  of  Solids 
and  Structures,  Vol.  5,  1969,  p.  191. 

2 5 Willis , J.  R. , HERTZIAN  CONTACT  OF  ANISOTROPIC  BODIES, 
i . Mechanics  and  Physics  of  Solids,  Vol.  14,  1966,  pp.  163-176. 

26Svelko,  V.  A.,  BOUSSINESQ-TYPE  PROBLEM  FOR  ANISOTROPIC 
SEMISPACE,  Prikladnaia  Mathematika  i Mekhanika,  Vol.  28,  196-4, 
pp.  908-913. 

27 Lekhnitskii,  S.  G.  , THEORY  OF  ELASTICITY  OF  AN  ANISOTROPIC 
ELASTIC  BODY,  Hoiden-Day  Inc.-,  San  Francisco,  1963. 

28Green,  A.  E. , and  Zerna,  W.  , THEORETICAL  ELASTICITY, 
Oxford  at  the  Claredon  Press,  19  54. 
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theories  that  exist,  no  one  theory  has  been  shown  as  THE  theory  for 
predicting  the  strength  and  failure  modes  of  composites  subjected  to 
multiaxial  and  combined  loading. 

Stress -Wave  Effects  in  Composite  Materials 

For  certain  impact  problems  involving  composites,  stress-wave  effects 
have  to  be  considered  in  studies  of  impact  response.  This  is  particularly 
true  for  high  impact  velocities  (v  > 1000  ft/sec).  Although  no  studies 
of  stress -wave  effects  were  conducted  in  the  program,  nevertheless  a 
brief  survey  was  made  of  the  problem. 

A thorough  review  of  studies  on  stress -wave  propagation  in  composite 
materials  is  given  by  Moon23,  The  review  covers  the  period  up  to 
December  1972.  The  major  properties  of  waves  in  composites  are 
discussed,  and  the  major  experimental  results  in  this  field  are 
summarized.  Various  theoretical  models  for  analyzing  wave 
propagation  in  laminated,  fiber-  and  particle -reinforced  composites 
are  surveyed.  The  anisotropic,  dispersive,  and  dissipative  properties 
of  stress  pulses  and  shock  waves  in  such  materials  are  reviewed.  A 
brief  review  of  the  behavior  of  composites  under  impact  loading  is 
presented  together  with  the  application  of  wave  propagation  concepts 
to  the  determination  of  impact  rtresses  in  composite  plates.  An 
extensive  bibliography  containing  211  entries  is  given  in  the  report. 
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SECTION  II 

| THEORY  DEVELOPMENT  FOR  PREDICTING  THE  RESPONSE  OF 

I COMPOSITES  TO  IMPACT  LOADING 

f The  approach  used  in  studying  the  response  of  composite  targets  to 

| _ impact  by  foreign  objects  was  similar  to  the  approach  described  in 

I References  10  and  46.  Whereas  the  latter  references  dealt  primarily 

with  the  response  of  isotropic  and  planar  isotropic  materials  subjected 
to  impact  by  spherical  impactors,  the  results  presented  here  apply  to 
generally  orthotropic  targets  (which  can  be  solids  of  revolution) 
impacted  by  an  isotropic  or  composite  impactor  in  the  form  of  a body  of 
revolution.  The  approach  that  was  employed  is  shown  in  schematic 
form  in  Figure  1,  and  consists  of  three  main  steps:  (1)  determination 
of  time- dependent  surface  pressure  distribution  under  an  impacting 
particle,  (2)  determination  of  internal  stresses  in  a generally  ortho- 
tropic target  caused  by  the  surface  pressure,  and  (3)  determination  of 
I failure  modes  in  the  target  caused  by  the  internal  stresses. 

The  pressure  distribution  under  an  impacting  particle  was  obtained  by 
■ analytically  combining  the  dynamic  solution  to  the  problem  of  impact 

of  bodies  with  the  static  solution  for  the  pressure  between  two  bodies  in 
• contact.  Having  the  time- dependent  surface  pressure,  the  time- 

I dependent  triaxial  stresses  in  targets  made  of  composite  materials 

were  determined  using  finite-element  computer  solution.  Knowing  the 
i internal  triaxial  stresses,  the  failure  modes  and  failure  envelopes  were 

' determined  by  applying  failure  criteria  for  generally  orthotropic  solids 

j to  the  stress  state  in  the  target.  The  assumptions  made  in  the  theory 

f development  were  as  follows:.  (1)  both  the  target  and  projectile  are 

linear  elastic,  (2)  impact  duration  is  long  compared  to  time  intervals 
in  stress-wave  response,  and  (3)  the  projectile  impact  is  normal  to  the 
i target  surface.  Assumption  (2)  was  made  only  in  studies  of  impact 

response  of  flat  composite  plates.  However,  the  studies  on  the  impact 
response  of  composite  cylinders  did  consider  stress-wave  effects. 


46.  Greszczuk,  L.  B. , MECHANICS  OF  PARTICLE  IMPACT, 

McDonnell  Douglas  Astronautics  Company  Report,  MDC,  G2612, 
March  1973. 
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IMPACT  CAUSES  PRESSURE  PRESSURE  CAUSES 


(PRESSURE  IS  TIME-  STRESSES 

DEPENDENT) 


STRESSES  CAUSE 
FAILURE 


l 


FAILURE  MODES  INTERACT 
TO  CAUSE  DAMAGE 


Figure  1.  Essential  Features  of  the  Approach. 

’ IMPACT  BETWEEN  TWO  ORTHOTROPIC  BODIES  OF  REVOLUTION 

j To  enable  determination  of  stresses  and  failure  modes  resulting  from 

an  impact  of  an  arbitrary  body  of  revolution  into  a flat  or  curved  target 
i made  of  composite  materials,  a generalized  solution  has  been  obtained 

for  the  area  of  contact,  maximum  surface  pressure,  pressure  distribu- 
i tion,  and  impact  duration  for  two  arbitrary  bodies  of  revolution.  The 

results  are  readily  applicable  to  special  cases  such  as  impact  of  a 
; sphere  or  cylinder  into  flat,  curved,  and  doubly  curved  semi-infinite 

targets.  The  solutions  to  the  above- noted  problems  were  obtained  by 
I analytically  coupling  the  dynamic  (non- stress-wave)  solution  for  the 

problem  of  impact  of  bodies  with  the  static  solution  for  the  pressure 
between  two  bodies  in  contact.  46,  52,  53 


52.  Timoshenko,  S.  , THEORY  OF  ELASTICITY,  McGraw  Hill  Book 
Company,  New  York,  1934. 

53.  Belajef,  W.  M.  , MEMOIRS  ON  THEORY  OF  STRUCTURES, 
Izdatelstvo  Puti  (St.  Petersburg),  1924.  See  also  Belajef, 
STRENGTH  OF  MATERIALS,  Gosudarstvennoje  Izdatelstvo 
Tekh.  -Teor.  1945  (In  Russian). 


26 


■" a *"1, «vv,*-9jyi5> • 


>r*5 ,',:w3PfTr; -*■ 


Area  of  Contact  Between  Two  Impacting  Bodies 

If  a solid  (or  impactor  designated  by  Subscript  1)  having  the  maximum 
and  minimum  radii  of  curvature  of  Rj^  and  Rj^j  is  pressed  by  a force 
P into  a target  (designated  by  Subscript  2)  having  the  maximum  and 
minimum  radii  of  curvature  Rjjm  an<^  ^2M*  area  contact  will  be 
elliptical,  with  major  and  minor  axes  of  ellipse  being 


a = m 


[xP(kl  + k2>  CR] 

b = -[t  p<kl  + k2>  CR] 


1/3 


(1) 

(2) 


where  is  a term  that  takes  into  account  the  curvature  effect 


b 


klm 


R 


2m 


klM 


2M 


(3) 


kj  and  k£  are  the  parameters  that  take  into  account  the  properties  of  the 
impactor  and  the  targets  and  are  defined  later;  m and  n are  parameters 
that  are  a function  of  Rim,  R-1M»  ®2m  an(*  &2M  The  numerical 

values  for  m,  n,  and  s (which  is  used  later)  are  shown  in  Figure  2. 
These  are  plotted  as  a function  of  9 where 


arc  cos 


|Cr[(RL'R!m)  +(RL'R2m) 

4)co82i } 


(4) 


U-  - J~\  / > 


and-4is  the  angle  between  planes  containing  Rjm  and  R^. 

If  the  two  bodies  are  pressed  together  by  a force  P,  the  combined 
deformation  of  both  bodies  can  be  expressed  in  terms  of  similar  param- 
eters. The  maximum  deformation,  is 


CKj  = S 


2 2 2 *i  ^ ^ 

9ir  P^kj  + k2P  ' 


256  C 
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(5) 
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Figure  2.  Values  of  Parameters  m,  n. 
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The  total  force  P resulting  from  impact  can  be  related  to  the  approach 
velocity  v of  the  two  bodies 

v = V1  + v2  (* 


using  an  approach  similar  to  that  presented  in  References  52  and  46 

r ? -i3/5 

p . n2/5 

4N, 


where 


256  CR 

9S3  rr2  (kj  + k2)2 


N.  = — + — 
1 ml  m2 


and  mj  and  are  the  masses  of  the  impactor  and  the  target. 

The  major  and  minor  axes  of  the  elliptical  area  of  contact  and  the 
maximum  deformation,  arj,  can  now  be  expressed  in  terms  of  impact 
velocity,  properties  of  target,  and  the  impactor  and  their  geometries 
by  substituting  Equation  (7)  into  Equations  (1),  (2),  and  (5): 


m ^<k,  ty  crn2/5 


(™,j 


b.„  2s (k,  +k2)  crn2/5 


QTj  = S 


9r2(k,  » k,)2l  f 4/51[5v2 
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Surface  Pressure  and  Ita  Distribution 


For  the  static  contact  problem,  the  relationship  between  the  maximum 
surface  pressure  acting  at  the  center  of  the  area  of  contact  and  the  force 
P is  given  as  52. 


_ 3P 
o ~ 2ffab 


(13) 


In  the  case  of  impact  problem  the  maximum  surface  pressure  q can  be 
obtained  by  combining  Equations  (7),  (10),  and  (11):. 


(14) 


C "3 

whereas  the  pressure  distribution  is  given  by  3 

J/2 


q - q0  1 


(15) 


where  x and  y are  the  coordinate  axes  in  the  directions  of  the  axes  of 
ellipse  a and  b,  respectively. 


Equations  ( 14)  and  (15)  are  the  final  equations  for  the  maximum  surface 
pressure  and  its  distribution  over  the  area  of  contact  resulting  from  an 
impact  of  isotropic  or  composite  impactor,  in  a form  of  body  of  revolu- 
tion, into  an  isotropic  or  composite  target,  also  being  a body  of  revolu- 
tion. The  terms  entering  these  equations  take  into  account 


1.  Shape  or  curvature  effects  (through  C^). 

c.  Properties  of  impactor  and  the  target  (through  kj  and  k^). 

3.  Impact  or  approach  velocity  (through  v = v^  + v^). 

4.  Masses  of  the  impactor  and  the  target  (through  Nj). 

Equations  (10)  and  (11)  give  the  axes  of  the  ellipse  over  which  the  pres- 
sure acts,  while  Equation  (12)  gives  the  deformation  at  the  center  of  the 
area  of  contact  resulting  from  impact  between  two  bodies. 
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Relationship  Between  kj  and  k2  and  Elastic  Properties 


The  terms  kj  and  k£  appearing  in  various  equations  of  the  previous 
subsections  take  into  account  the  mechanical  properties  of  the  impactor 
and  the  target.  For  the  case  of  the  planar  isotropic  composite  target  (see 
Figure  4),  k£  can  be  obtained  from  Reference  54*  and  is 


Vgl  (<  VA11  A22  + g2>2  - <A12 + g;>21 

2ir  G,  ( Aj  j A^2  - Aj2) 


1/2 


(16) 


where 
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I 


J 

\ 


5 

\ 


A12  “ E1  V12  P 


P = 


1 " vi  * Zvnb 


6 =El/E2 


(17) 


and  E,  G,  and  u are  the  Young’s  modulus,  shear  modulus,  and  Poisson's 
ratio  of  the  target,  while  Subscripts  1 and  2 denote  the  radial  and  thick- 
ness directions,  respectively,  2 being  in  the  direction  of  impact.  For  a 
pianar  isotropic  material  the  properties  in  the  1 plane  are  independent  of 
the  orientation.  It  can  be  readily  shown  that  for  isotropic  materials 
Equation  (16)  reduces  to  the  well-known  equation  52. 


k = 


(18) 


*The  equations  given  in  Reference  54  have  several  errors,  which 
were  corrected  to  arrive  at  Equation  (16). 


54.  Conway,  H.  D. , THE  PRESSURE  DISTRIBUTION  BETWEEN  TWO 
ELASTIC  BODIES  IN  CONTACT,  Zeitschrift  fuer  AngeWaadte 
Mathematik  und  Physik,  Vol.  VII,  1956,  pp.  460-465. 
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If  the  impactor  is  also  made  of  material  that  is  planar  isotropic,  then 
kj  will  be  defined  similar  to  ki,  except  that  kj  will  be  a function  of  the 
elastic  properties  of  the  impactor. 

As  is  readily  apparent  from  Equations  (10),  (11),  (12),  and  (14),  the 
parameters  kj  and  ko  influence  the  shape  of  the  area  of  contact  as  well 
as  the  magnitude  of  uie  surface  pressure.  Although  no  expression  has 
been  found  for  k2  for  generally  orthotropic  solids  and  its  derivation 
appears  to  be  extremely  complex,  it  has  been  noted  (based  primarily 
on  experimental  results)  that  when  an  isotropic  spherical  body  is  pressed 
into  a unidirectional  orthotropic  target,  the  area  of  contact  is  slight 
elliptical,  which  indicates  the  dependence  of  k2  on  inplane  fiber  orienta- 
tion. An  approximate  expression  for  k2  for  a generally  orthotropic 
material  can  be  obtained  from  Equations  (16)  and  (17)  by  making  the 
following  substitutions: 


! 
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4 * 4 

,-l  cos  * , sin  + 1_ 

Et  E„  + 4 IG 
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where  for  a fiber- reinforced  composite,  E^  is  the  inplane  Young's 
modulus  in  the  fiber  direction,  E-p  is  the  inplane  Young's  modulus  in 
the  transverse  direction,  Clt  is  the  inplane  shear  modulus,  i>lt  is  the 
major  inplane  Poisson's  ratio,  and  ♦ is  the  angle  between  the  fiber 
direction  and  any  given  inplane  direction.  Figure  3 shows  the  nondimen- 
sionalized  shape  of  the  area  of  contact  in  a generally  orthotropic  solid 
as  well  as  a circle  of  equivalent  area.  Shown  there  also  are  die  material 
properties  on  which  the  results  were  based.  These  material  properties 
are  for  a unidirectional  composite  consisting  of  Courtaulds  HTS  graphite 
fibers  and  an  epoxy  resin.  The  important  thing  to  note  from  Figure  3 
is  that  for  even  a high  degree  of  inplane  orthotropy  the  area  of  contact 
deviates  only  slightly  from  circular  shape.  As  noted^®,  properties  that 
have  the  greatest  influence  on  k2  are  the  properties  associated  with  the 
thickness  direction,  that  is,  the  direction  of  impact. 

Duration  of  Impact 

The  maximum  pressure  qQ  occurs  at  a time  1/2  tQ,  where  tQ  is  the 
impact  duration.  The  latter  can  be  obtained  using  an  approach  similar 
to  that  presented  in  Reference  52.  From  the  problem  of  impact  of  two 
bodies 

1 ..  2 2.  2M.T  2/5  |«.v 

2 (o  - v ) = - jNNj  or  (21) 
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Figure  3.  Actual  and  Approximate  Shape  oi  the  Area  of  Contact 
Resulting  from  an  Indentation  or  Impact  cf  a Rigid 
Spherical  Body  Into  Generally  Orthotropic  Composite 
Material. 
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where  the  various  terms  have  been  defined  previously.  The  value  2.  94 
comes  from  the  evaluation  of  the  integral  appearing  in  Equation  (25). 
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The  pressure  qj  and  the  approach  velocity  v^  at  any  given  time  tj  can  be 
obtained  by  numerically  integrating  Equation  (25).  Knowing  the  relation- 
ship between  vj,  v,  tj,  and  t0,  the  pressure  q^  and  the  principal  axes  of 
the  area  of  contact  a^  and  b^  at  any  given  time  tj  can  be  obtained  from 
Equations  (10),  (11),  and  (14)  by  replacing  v with  v^. 

INTERNAL  STRESSES  IN  FLAT  MULTILAYER  COMPOSITE  PLATES 
CAUSED  BY  SURFACE  PRESSURE  q RESULTING  FROM  FOREIGN 
OBJECT  IMPACT 

Knowing  the  time- dependent  surface  pressure,  its  distribution,  the  area 
over  which  it  acts,  and  the  variation  of  pressure  and  area  of  contact  as 
a function  of  time,  it  is  now  possible  to  determine  the  internal  stresses 
in  composite  targets  as  a function  of  the  variable  noted  above.  The  internal 
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stresses  in  planar  isotropic  targets  were  determined  using  the  SAAS  III 
finite- element  computer  code  5 5,  whereas  in  generally  orthotropic 
solids  the  internal  stresses  were  determined  using  the  ASAAS  finite- 
element  computer  code  ^ which  was  modified,  as  described  later  on, 
to  handle  the  problem  of  generally  orthotropic  solids. 

Internal  Stresses  in  Planar  Isotropic  Targets 


The  internal  stresses  in  planar  isotropic  targets  subjected  to  impact- 
induced  surface  pressure  were  determined  using  the  SAAS  III  finite- 
element  computer  code,  which  was  available  at  McDonnell  Douglas 
Astronautics  Company  (MDAC)  and  operational  prior  to  initiation  of 
this  program.  An  example  of  a planar  isotropic  material  that  could  be 
analyzed  using  SAAS  III  is  a multilayer  composite  plate  in  which  the 
inplane  properties  of  any  given  layer  (Young's  modulus,  Poisson's  ratio, 
and  shear  modulus)  are  different  from  the  corresponding  properties  in 
the  thickness  direction  of  the  plate.  However,  the  inplane  properties 
of  the  plate,  although  not  direction  dependent,  can  be  different  for  various 
layers,  as  shown  in  Figure  4.  Figures  5 and  6 show  some  typical  results 
on  the  internal  stresses  in  planar  isotropic  targets  subjected  to  surface 
pressure  q resulting  from  foreign  object  impact.  The  stresses  shown 
there  were  normalized  with  respect  to  qQ,  whereas  the  distribution  of 
stresses  in  the  target  was  normalized  with  respect  to  a,  the  radius  of 
the  area  of  contact.  Although  only  the  results  for  normal  and  radial 
stresses  are  shown  in  Figures  5 and  6,  similar  results  have  been 
obtained  for  the  circumferential  and  shear  stresses. 

Internal  Stresses  in  Generally  Orthotropic  Targets 

Although  SAAS  III  can  handle  special  cases  of  multilayer  materials  (see 
Figure  4)  subjected  to  impact- induced  pressure,  it  cannot  handle  the 
case  of  a generally  orthotropic  multilayer  solid  such  as  shown  in  Fig- 
ure 7.  The  three-dimensional  stress  state  in  the  latter  type  of  material 
was  determined  by  modifying  the  ASAAS  finite- element  computer  code, 
which  is  the  most  recent  version  of  SAAS  HI  and  is  applicable  to  mater- 
ials with  polar  orthotropy.  To  become  applicable  to  the  multilayer 
generally  orthotropic  solids  having  orthogonal  orthotropy,  the  code  was 
modified  as  described  below. 


55.  Crose,  J.  G. , and  Jones,  R.  M. , SAAS  HI,  FINITE  ELEMENT 
STRESS  ANALYSIS  OF  AXISYMMETRIC  AND  PLANE  SOLIDS  WITH 
DIFFERENT  ORTHOTROPIC,  TEMPERATURE- DEPENDENT 
MATERIAL  PROPERTIES  IN  TENSION  AND  COMPRESSION, 
TR-0059  (S6816-53)  - The  Aerospace  Corporation,  San  Bernadino, 
California,  June  1971. 

56.  Crose,  J.  G. , ASAAS  - ASYMMETRIC  STRESS  ANALYSIS  OF 
AXISYMMETRIC  SOLIDS  WITH  ORTHOTRO.PIC,  TEMPERATURE- 
DEPENDENT  MATERIAL  PROI  ERTIES  THAT  CAN  VARY  CIRCUM- 
FERENTIALLY, The  Aerospace  Corporation  Report 
TR-0172(S28I6-15)-  1,  December  1971. 
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Figure  4.  Planar  or  Transversely  Isotropic  Multilayer  Material. 


The  system  of  equations  solved  by  SAAS  or  the  more  recent  version 

ASAAS,  is  of  the  form 

n r r t 

N 

' r T 

II  I (a11)  (Cn)  (a11)  dv 

{u}  = E 

/ (a11)  {rn}  dv 

n=l[y  vol 

n=  1 

J vol 

N 

r r t 

1 

+ E 

/ (dn)  {pn}  dA 

n=l 

J area 

where 

(a11)  = matrix  expressing  the  relation  between  the  strain  matrix 

( tn ) of  any  element  n and  the  nodal  displacement  matrix  (u) 
T 

(a  ) = transpose  of (an) 
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gure  5.  Internal  Stress  Distribution  in  Semi  - Infinite  Multilayer  Planar  Isotropic 
Material  Resulting  from  Surface  Pressure  Caused  by  Foreign  Object 
Impact  (E.  /E?  = 0.  35). 


Figure  6.  Internal  Stress  Distribution  in  Semi-Infinite  Multilayer  Planar  Isotropic 
Material  Resulting  from  Surface  Pressure  Caused  by  Foreign  Object 
Impact  (E,  /E?  = 2.  86). 


(Cn)  = matrix  for  the  material  properties  (constitutive  relations) 


(rn)  = thermal  load  matrix 

I- 

(Pn)  = surface  traction  matrix 

The  expressions  on  the  right  side  of  Equation  (28)  represent  the  thermal 
; , loads  (first  term)  and  surface  tractions  (second  term),  Thf  coefficient 

of  (U)  on  the  left  side  of  Equation  (28)  is  the  stiffness  matrix.  The 
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matrix  (Cn)  giving  the  constitutive  relations  for  the  material  is  the 
inversion  of  (A),  which  is  defined  as 


L 

* 

J 

all 

a12 

a13 

0 

0 

0 

l 

a21 

a22 

a23 

0 

0 

0 

\ 

\ 

(A)  =■ 

a31 

a32 

a33 

0 

0 

0 

& 

0 

0 

0 

a44 

0 

0 

! 

* 

i. 

E 

0 

0 

0 

0 

a55 

0 

r 

; i 

0 

0 

0 

0 

0 

*66  . 

(29) 


where,  for  material  such  as  shown  in  Figure  4 (transversely  isotropic), 
1 


all  a33  ~ E 


1 


22 


_1_ 

E_ 


_ J JL 

a12  = a21  a23  ' a32  E. 


a13  = a31  = E, 


a44  = a66  G, 
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Ej  = inplane  Young's  modulus 

E^  - Young's  modulus  associated  with  thickness  direction 
Yj  = inplane  Poisson's  ratio 

Dj2  = Poisson's  ratio  associated  with  Direction  1-2  or  3-2. 
Gj  = inplane  shear  modulus 

G^  = shear  modulus  associated  with  Direction  1-2  or  3-2. 
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n 


In  the  existing  version  of  the  finite- element  computer  code,  the  material 
properties  are  independent  of  the  0- coordinate  in  the  plane  of  the  laminae 
Therefore,  although  the  existing  code  could  handle  some  special  cases 
of  multilayer  composite  plates,  it  was  inadequate  in  handling  the  influ- 
ence o t fiber  orientation  or  of  stacking  sequences.  The  code  was  modi- 
fied to  enable  it  to  handle  the  latter  types  of  problems.  Instead  of  using 
the  constitutive  relationship  given  by  (29),  the  following  was  used  in  the 
governing  equation  (Equation  (28)): 


all 

a12 

a13 

0 

al  5 

0 

a21 

a22 

a23 

0 

a25 

0 
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a32 

a33 

0 

a35 

0 

0 

0 

0 

a44 

0 

a46 

a15 

a25 

a35 

0 

a55 

0 

0 

k. 

0 

0 

a64 

0 

a66 

(30) 


where 


all  = 1^8  1*55  + 2a13  + 3 (an  * a33^J  + 1/2  (ajj  - i33)  Cos  28 
+ 1/8  (ajj  + a^3  - a^5  - 2a^3)  Cos40 

a12  =1/2  (a ^ + *23)  + 1/2  (a^  - a^3)  Cos20 

a15  = [au  + 1/8  (an  + a33  - 2a13  - a55)] 

- 1/8  (a^  + a33  - 2a^3  - a,.,.)  Cos40 
= 1/2  (a33  - a^j)  Sin20  - 1/4  (a^j  + a33  - a^3  - 2aj3)  Sin40 

a22  = a22 

*23  = 1/2  (a12  + a23)  + 1/2  (a23  - aJ2)  Cos20 


a25  = (a23  ' a12}  Sin  ZG 

a33  = 1/8  [a33  + 2aj3+3(ajj  + *33)]  + 1/2  (a.  3 - a^)  Cos20 


+ 1/3  (a^j  + a33  - a^3  - 2a^)  Cos40 
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= 1/2  (a^g  - ajj)  Sin29  +1/4  (a^  + a^g  - agg  - 2ajg)  Sin40 


a44  = 1/2  (a44  4 a66)  4 1/2  (a44  - a66)  Co,2S 
a4^  = ^ /2  (a^g  - a^^)  Sin20 

a55  = [a55  + 1/2  (an  + a^  - 2aJ3  - agg)] 

-1/2  (a^j  + a33  - 2ajg  - agg)  Cos40 

a66  ~ <a44  f a66^  + ^a66  ' a44*  Cos2e 

and  the  terms  ajj,  aj2»  a13»  a21*  a22*  a23*  a31»  a32*  a33* 

are  related  to  the  engineering  elastic  constants  of  a material  in  the  three 
orthogonal  directions.  For  example,  for  material  such  as  shown  in 
Figure  7,  the  various  terms  now  become 


where  El,  E2,  and  E3  are  the  Young’s  moduli  of  the  material  in  the 
three  orthogonal  directions,  each  of  them  being  different;  that  is, 

| E1  * Ez  * E3 

i 

I 

j The  above-described  modification  of  ASAAS  enabled  the  following  types 

of  problems  to  be  analyzed: 

1.  Local  and  overall  (plate  bending)  multiaxial  stresses  resulting 
from  impact  in  multilayer  generally  orthotropic  plates  (or 
solids)  made  of  layers  having  different  elastic  properties. 

2.  Multiaxial  stresses  in  multilayer  composite  plates  (or  solids), 

i with  layers  having  different  fiber  orientations  (fibers  in  any 

given  layer  have  orientation  ± 9 where  0°<  9 <90°). 

3.  Variations  in  local  and  overall  multiaxial  stresses  due  to 
stacking  sequence. 

The  modified  version  of  ASAAS  gives  a three-dimensional  solution  for 
the  internal  stresses  in  generally  orthotropic  solids  having  orthogonal 
orthotropy.  Inasmuch  as  closed- form  rigorous  elasticity  solutions 
could  not  be  found  for  this  type  of  problem,  the  modified  ASAAS  com- 
puter code  was  checked  out  by  applying  it  to  several  simple  problems 
involving  isotropic  and  orthotropic  materials.  The  problems  involved 
determination  of  stresses  and  deflections  in  isotropic  and  orthotropic 
circular  plates  fixed  along  the  outer  boundary  and  subjected  to  uniform 
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surface  pressure.  The  latter  problems  have  been  selected  for  checkout, 
because  closed-form  solutions  are  available  for  them.  Figure  8 shows 
a comparison  of  solutions  for  deflections  in  an  isotropic  circular  plate 
subjected  to  uniform  surface  loading.  Presented  therein  are  the  results 
from  a closed- form  solution  as  well  as  from  SAAS  III  and  ASAAS  com- 
puter solutions  in  which  two  different  grid  sizes  were  employed.  The 
closed-form  solution  for  the  problem  is  given  in  Reference  57.  For  a 
given  grid  size,  both  the  ASAAS  and  SAAS  III  computer  codes  give 
identical  results  for  the  deflection.  Inasmuch  as  SAAS  III  has  been 
checked  out  previously^,  this  gives  a good  indication  of  the  accuracy 
of  ASAAS.  By  decreasing  the  size  of  the  grid  elements,  the  ASAAS 
solution  is  shown  to  approach  the  closed-form  solution.  Further  decrease 
in  the  grid  elements  did  not  influence  the  results  significantly,  as  dis- 
cussed later  in  this  section.  The  slight  discrepancy  (*5%)  between  the 
closed-form  and  ASAAS  solution  using  Model  B grid  is  not  so  much  an 
indication  of  the  error  in  computer  solution  as  an  indication  of  error  in 
the  closed-form  solution.  The  computer  solution  solves  the  three- 
dimensional  elasticity  problem,  whereas  the  closed-form  solution  was 
obtained  for  a two-dimensional  plate  problem  ®7,  an(j  j^s  accuracy  is 
governed  by  the  assumptions  made  in  the  plate  theory. 

The  next  checkout  problem  considered  involved  determination  of 
deflections  and  stresses  in  an  orthotropic  circular  plate  fixed  along  the 
boundary  and  subjected  to  a uniform  surface  loading.  An  approximate 
closed-form  solution  (based  on  plate  theory  formulation)  for  this  problem 
is  given  in  Reference  58.  For  this  sample  problem  the  material  pro- 
perties associated  with  the  axes  of  symmetry  of  material  (p=0  and  p=90°) 
were  taken  as 

E^  = E,p  = E^  = 1.  32  x 10^  psi 
GLT  = GLZ  = GTZ  = °*2  x 10  P31 


V 


j 

I 


LT  = VTL  " 


u - 

LZ 


0.  10 


The  variation  of  the  Young’s  moduli,  shear  modulus,  and  Poisson's 
ratio  as  a function  of  orientation  (Angle  P)  is  shown  in  Figures  9 and  10. 
Figure  11  shows  a comparison  of  deflections  obtained  from  the  closed- 
form  solution  ^8  antj  ASAAS  computer  solutions  using  different- sized 
grid  elements.  As  is  readily  seen,  going  from  Model  B grid  to  Model  C 


57.  Timoshenko,  S.  P. , THEORY  OF  PLATES  AND  SHELLS,  McGraw- 
Hill  Book  Company,  Nev  York,  1940,  p.  60. 

53.  Hearmon,  R.  F.  S. , APPLIED  ANISOTROPIC  ELASTICITY,  Oxford 
University  Press,  Oxford,  1961,  p 126. 
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DEFLECTION,  IN 


Q. 

q h-  T q » 1000  LB/IN? 

a - 2.0  IN. 

H h h = 0.2  IN. 

E • 1.32  x 106  PSI 

v * 0.10 

6 » E/2(l+v)  * 0.60  x 106  PSI 

— CLOSED-FORM  SOLUTION  [57] 

ASMS  SOLUTION  USING  MODEL  A 

00  0 SMS  III  SOLUTION  USING  MODEL  A 
ASMS  SOLUTION  USING  MODEL  B 


Figure  8.  Comparison  of  Solutions  for  Deflections  in  a Circular 
Isotropic  Plate  Subjected  to  Uniform  Surface  Pressure. 


CIRCULAR  PLATE  UNDER 
UNIFORM  LOADING 
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Figure  9.  Variation  of  Young's  Moduli  and  Shear  Modulus 
as  a Function  of  (3. 
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Figure  10.  Variation  of  Poisson's  Ratios  as  a Function  of  (3 
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Comparison  of  Solutions  for  Deflections  in  a Circular 
Orthotropic  Plate  Subjected  to  Uniform  Surface 
Pressure. 
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grid  does  not  influence  the  results  significantly.  Consequently,  Model  B 
grid  was  selected  for  use  in  the  application  phase  of  the  program 
described  in  Section  111.  The  ASAAS  computer  running  times  for  the 
various  grid  models  were  as  follows: 

Model  A:  0.  034  Machine  Units 

Model  B:  0.065  Machine  Units 

Model  C:  0.  095  Machine  Units 

The  fact  that  the  deflections  obtained  from  ASAAS  are  higher  than  the 
closed- form  solution  (the  opposite  was  true  in  the  case  of  isotropic 
plate)  is  attributed  to  the  influence  of  low  shear  moduli  that  were 
assumed  for  the  sample  problem. 

The  radial  and  circumferential  stresses  obtained  from  the  closed- form 
solution  and  ASAAS  computer  solutions  using  different  grid  models  are 
compared  in  Figure  12.  As  was  the  case  with  deflections,  going  from 
grid  Model  B to  grid  Model  C does  not  significantly  influence  the  results. 

The  modified  and  checked-out  code  was  used  in  determining  internal 
triaxial  stresses  in  several  sample  problems  involving  generally  ortho- 
tropic solids,  prior  to  applying  the  code  to  the  problems  of  interest. 
Tables  I and  II  show  some  typical  results  for  the  normal  (trz)  and  radial  (trr) 
stresses  in  a sample  problem  shown  in  Figure  13.  Although  only  <rz  and 
<rr  are  shown  in  Tables  I and  II,  similar  results  have  also  been  obtained 
for  the  circumferential  (erg)  and  shear  stresses  (crz  andogz)  at  various 
points  within  the  target.  The  results  shown  in  Tables  I and  II  are  for 
HTS  graphite- epoxy  composites  having  properties  as  shown  in  Figure  3. 


INTERNAL  STRESSES  IN  MULTILAYER  ORTHOTROPIC  CYLINDERS 
CAUSED  BY  SURFACE  PRESSURE  q RESULTING  FROM  FOREIGN- 
OBJECT  IMPACT 

The  impact- induced  stresses  in  multilayer  orthotropic  cylinders  were 
determined  using  the  CYLINDER  computer  code  ^9.  This  code  uses 
Sander’s  shell  equations  ^0  an(j  an  explicit  finite- difference  integration 


59.  Chane,  H.  L. , ELASTIC  PLASTIC  ASYMMETRIC  TRANSIENT 
RESPONSE  OF  ORTHOTROPIC  CYLINDRICAL  SHELLS  SUBJECTED 
TO  MECHANICAL  AND  THERMAL  LOADS,  McDonnell  Douglas 
Astronautics  Company/ Western  Division  MDC  G 3750  (in  preparation). 

60.  Sanders,  Jr.,  J.  L. , NONLINEAR  THEORIES  FOR  THIN  SHELLS, 
Quarterly  Applied  Mathematics,  Vol.  21,  No.  1,  April  1963,  pp.  21-36. 
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Figure  12.  Comparison  of  Solutions  for  Stresses  in  a Circular 
Crthotropic  Plate  Subjected  to  Uniform  Surface 
Pressure. 
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Figure  13.  Sample  Problem  Used  in  Modified  ASAAS  Code. 
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scheme.  It  accounts  for  orthotropic  elastic  perfectly  plastic  material 
behavior  by  using  a modified  Hill-von  Mises  yield  condition  and  flow 
rule,  and  for  isotropic  elastic  strain-hardening  material  behavior  by 
using  a mechanical  sublayer  model  as  well  as  von  Mises  yield  con- 
dition and  flow  rule.  The  calculation  proceeds  in  incremental  time 
steps  using  an  explicit  integration  scheme.  The  shell  thickness  is 
modeled  as  equal  area  and  equally  spaced  layers  connected  by  infinitely 
thin  bond  layers  of  infinite  transverse  shear  rigidity.  Strain  increments 
are  calculated  using  Kirchoff  s assumptions  for  deformations.  The 
finite-difference  mesh  used  represents  one  fourth  of  the  shell,  employ- 
ing symmetry  conditions  on  a half  circle  and  symmetry  conditions  on  a 
half  length. 

The  CYLINDER  code  was  originally  developed  by  the  Kaman  Sciences 
Corporation  ^2,  63#  The  code  in  operation  at  MDAC  ^9  was  modified  to 
eliminate  various  p oblems  that  were  encountered  in  implementing  the 
code  and  to  improve  its  analysis  capabilities  for  vulnerability  and  hard- 
ness analyses.  Modifications  or  improvements  included: 

• Nonconvergence  in  the  iterative  plasticity  calculations  due  to  a 
poor  initial  guess  of  the  plasticity  proportionality  factor  used 
for  determining  the  plastic  strain  increments  was  eliminated 
by  incorporating  a quadratic  function  to  determine  a much  more 
accurate  initial  value  of  the  proportionality  factor. 

• A set  of  consistent  nonlinear  equilibria  n equations  was  added 

to  replace  linear  equilibrium  equations  that  produced  inaccurate 
results  v/hen  deflections  became  moderately  large.  These 
equations  were  taken  from  the  aforementioned  Sander's 
formulation. 

• Strain- displacement  relations  were  modified  to  include  the 
effect  of  inplane  displacements  on  shell  surface  tangent 
rotations  removing  the  Donnell- Mushtari- Vlasov  Approxima- 
tion originally  in  the  code,  and  the  effect  of  rotations  about 
the  shell  surface  normal  was  also  included. 

• For  isotropic  materials,  kinematic  hardening  (Bauschinger 
effect)  and  isotropic  hardening  materials  models  were 
incorpc  rated. 
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• Viscous  damping  behavior  was  added. 

During  this  program,  additional  modifications  and  improvements  were 
made  in  the  CYLINDER  computer  code,  including  addition  of  an  option 
to  input  different  orthotropic  properties  for  different  layers,  increasing 
of  mesh  size  from  15x15  to  16  x 50  segments,  and  addition  of  new  pres- 
sure input  option  to  facilitate  the  input  of  arbitrary  distributed  spot  loads 
resulting  from  impact.  In  contrast  to  the  plate  solution  described 
previously,  the  solution  obtained  from  the  CYLINDER  code  is  the  two- 
dimensional  stress  solution.  Moreover,  the  CYLINDER  code  gives  the 
dynamic  stresses  caused  by  impact,  including  stress-wave  effects, 
whereas  in  the  flat-plate  analysis  only  the  quasi-dynamic  stresses  can 
be  calculated. 

The  input  and  output  of  the  modified  CYLINDER  code  used  in  this  pro- 
gram are  listed  on  Page  55.  The  model  and  the  grid  used  in  the  code 
are  shown  in  Figure  14.  The  cylinder  wall  can  consist  of  a maximum 
of  eight  layers  made  of  orthotropic  materials.  Each  layer  can  have  a 
different  set  of  elastic  properties.  The  grid  contains  50  elements 
(51  nodes)  in  the  circumferential  direction  (half  the  circumference)  and 
16  elements  (17  nodes)  in  the  axial  direction  (half  the  cylinder  length). 

By  using  this  type  of  grid,  approximately  80%  of  the  available  CDC  6600 
computer  core  storage  is  utilized  (250  Kg  of  the  320  Kg  storage  locations 
are  used).  Use  of  finer  grid,  although  feasible,  is  not  practical  from  the 
standpoint  of  computer  running  time.  Figure  15  shows  the  idealized 
and  actual  area  of  contact  resulting  from  impact  by  a spherical  impactor. 
Shown  also  in  Figure  15  is  the  idealized  and  actual  load/pressure  distri- 
bution at  the  centerline  of  the  area  of  contact.  Figure  16  shows  the 
computer  idealization  and  actual  pressure-time  relationship. 

The  sample  problem  used  for  code  checkout  is  shown  in  Figure  17.  The 
problem  involved  determination  of  stresses  in  a multilayer  cylinder 
subjected  to  time-dependent  surface  pressure  caused  by  a 100-ft/sec 
impact  by  a spherical  impactor.  All  the  pertinent  information  for  the 
sample  problem  is  shown  in  Figure  17.  The  cylinder  consisted  of 
eight  layers.  The  fibers  in  the  three  layer?  on  the  outer  surface  and 
the  three  layers  on  the  inner  surface  were  oriented  in  the  circumferen- 
tial direction,  while  the  fibers  in  the  center  two  layers  were  oriented  in 
the  axial  direction.  The  results  shown  in  Figure  17  are  for  a cylinder 
made  of  HTS  graphite  fibers  and  epoxy  matrix.  For  such  a composite, 
the  mechanical  properties  associated  with  the  principal  directions  are: 

Young's  modulus  in  the  fiber  direction,  E^  = 21. 5 x 10^  psi 

Young's  modulus  in  the  transverse  direction,  F,^,  = 1.5x  10^  psi 

Inplane  shear  modulus,  G^.  = 0.7  x 10^  psi 

Principal  Poisson's  ratio,  = 0.25 
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LISTING  OF  INPUT  AND  OUTPUT 
FOR  THE  CYLINDER  COMPUTER  PROGRAM 

INPUTS  TO  THE  COMPUTER  PROGRAM 

1.  Cylinder  Length  L 

2.  Cylinder  Radius  R 

3.  Cylinder  Thickness  t 

a.  Thicknesses  of  layers  within  the  wall  thickness  (maximum  c*  8 
layers). 

b.  Orthotropic  elastic  properties  of  layers;  two  Young's  moduli 
(Ex,  Eq) , Poisson's  ratio  (vx0)  and  shear  modulus  (Gx0).  The 
properties  of  each  layer  can  be  different. 

c.  Density  of  each  layer 

4.  Pressure-time  History 

5.  Time  steps  for  which  calculations  are  to  be  made. 

6.  Total  response  time  and  edit  time  at  which  printout  is  obtained. 

OUTPUT  OF  THE  CYLINDER  CODE  AS  A FUNCTION  OF  TIME  AND  LOCAuON 

1.  Axial,  circumferential  and  shear  stresses  in  each  layer  (o^,  oQ  and  oXQ) 

2.  Axial,  circumferential  and  shear  strains  in  each  layer  ( cQ  and  eXQ) 

3.  Axial,  circumferential  and  shear  forces  acting  on  the  cylinder  wall 
(Nx,  Nq  and  NX0) 

4.  Axial,  circumferential  and  torsional  moments  acting  on  the  cylinder  wall 
(Mx.  M0  and  MX0) 

5.  Axial,  circumferential  and  radial  displacements  U,  V and  W,  respectively. 

6.  Axial  (U)  and  Radial  (W)  deflections  and  their  variation  at  the  crown  of 
the  cylinder  (6  = 0)  as  a function  of  axial  coordinate  (J  = 1,  I = 1,  2, 
17) 

7.  Circumferential  (V)  and  Radial  (VI)  deflections  and  their  variation  at 

centerline  of  cylinder  (X  = L/2)  and  their  variation  as  a function  of 
circumferential  coordiate,  s (I  = 1,  J = 1,  2 51). 
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jure  17.  Variation  of  Impact- Induced  Stresses  (at  s = 0)  as  a Function 
and  Axial  Coordinate,  x. 
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Tensile  strength  in  the  fiber  direction,  = 204  x 10  psi 

3 

Tensile  strength  in  the  transverse  direction,  cr  = 11.8x10  psi 

3 

Compressive  strength  in  the  fiber  direction,  <r^  = 201  x 10  psi 

3 

Compressive  strength  in  the  transverse  direction,  <r^c  = 34.4  x 10  psi 

3 

Shear  strength  (interlaminar),  t = 14  x 10  psi 

Figure  17  shows  the  circumferential  and  axial  stresses  at  the  inner  and 
outer  surfaces  of  the  cylinder  as  a function  of  axial  (x)  coordinate  and  a 
function  of  time,  t.  The  variation  of  stresses  as  a function  of  time  and 
circumferential  coordinate  is  shown  in  Figure  18.  The  maximum  stresses 
at  x = s = 0 (center  of  contact)  as  a function  of  time  are  shown  in  Fig- 
ure 19.  Figures  20  and  21  show  some  results,  based  on  maximum 
stress  criteria,  on  the  initiation  and  progression  of  failure  as  a function 
i of  time.  For  the  problem  considered,  the  total  contact  duration  was 

I calculated  to  be  tD  = 172  psec.  As  shown  in  Figure  20,  failure  of 

! material  initiates  in  the  transverse  direction  in  the  innermost  layer  of 

the  cylinder  at  t = 7.  92  psec.  The  failure  region  increases  with  increas- 
ing t until  at  t = 9.  84  usee,  failure  initiates  in  another  layer  (Mo.  7). 

As  time  increases  beyond  t = 11.04  psec,  the  damage  zone  grows  and 
failure  progresses  to  other  layers.  At  t = 15.4  psec,  critical  damage 
takes  place  in  the  layer  on  the  outer  surface  of  the  cylinder,  as  shown 
in  Figure  21.  At  that  time  increment,  compression  failure  takes  place 
in  the  direction  of  the  fibers.  The  damage  in  other  layers  at  the  critical 
time  increment  (t  = 15.  4 psec)  is  also  shown  in  Figure  21. 

The  results  for  stresses  and  deformation  were  obtained  at  0.  06-psec 
time  intervals.  Although  the  total  response  time  of  interest  is  172  psec 
(total  contact  duration),  only  the  early-time  response  (t  = 11.28  psec) 
was  obtained  in  the  printout,  primarily  to  minimize  the  computer 
] running  time  in  the  checkout  problem.  The  computer  running  time 

I (total  machine  units)  for  the  early-time  response  (t  = 11.28  psec)  was 

| ~15  minutes.  The  computer  performed  approximately  190  time-step 

| calculations  within  that  time. 

i 

| COUPLING  BETWEEN  IMPACT  PARAMETERS,  TARGET  PROPER  ! f£S, 

» AND  FAILURE  MODES 

The  final  step  in  the  theory  development  was  to  establish  failure  modes 
resulting  from  internal  triaxiai  stresses  caused  by  impact-induced 
surface  pressure,  as  well  as  the  time  sequence  for  occurrence  of  the 
various  failure  modes.  This  was  done  by  applying  the  failure  criterion 
for  generally  orthotropic  solids  to  the  impact- induced  triaxiai  stress 
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Figure  18.  Variation  of  Impact- Induced  Stresses  (at  x = 0)  as  a Function  of  Time 
and  Circumferential  Coordinate,  (s). 
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Figure  19.  Variation  of  Stresses  at  x = 0,  s = 0 as  a Function  of  Time. 
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Figure  21.  Failure  in  Various  Layers  at  t = 15.4  jxsec  When  Critical  Damage  Occurs 
(Failure  of  Layer  Number  1 in  the  Fiber  Direction). 


state  in  the  target.  The  following  failure  criterion,  based  on  distortion- 
energy  theory  *&,  was  used  to  predict  failure  and  failure  envelopes: 


2 2 2 2 2 2 

. *11  a2Z  a 33  *12  13  23 

1 = “l"  +~2_  + 2~  + 2~  +“T"  + ~ 2~ 

f F F F^  F^ 
*11  *22  * 33  12  13  23 


■m 


22  \ f(1  + Zv21  - V23)E1  + (1  + 2U12  ' ”l3)E2 


22/  (2  +1,12  + +U21  + y23)ElE2 


fjjA/  °33  \ (1  + 2u31  "U32)E1  * (I  + 2u13  ~ l'l2)E3  (31) 

^VV33/  |<2+"12+”13'  '2+“31+“32)ElE2|1/2 


/ <r22\ / °33  \ (1  + ZU23  " "21)E3  + (1  * 2u32  ~ U31)E2 

\FW\F33/p+“2i  +"23)  <2t“3l  +‘'3Z>E2E3l1/i 

where  <r'  s are  the  impact- induced  direct  and  shear  stresses,  F's  are 
the  allowable  strength  properties  of  material  associated  with  the  three 
orthogonal  directions,  and  E's  and  v' s are  the  Young's  moduli  and 
Poisson's  ratios  associated  with  Directions  1,  2,  and  3.  The  failure 
criterion  was  incorporated  into  the  SAAS  III  and  ASAAS  as  a subroutine. 
To  take  into  account  differences  in  tensile  and  compressive  strength  in 
any  given  direction,  an  allowance  was  made  in  the  computer  subroutine 
that  allowable  compressive  strength,  Fc,  is  used  if  the  corresponding  a 
is  negative,  and  allowable  tensile  strength,  Ft,  is  used  if  o is  positive. 


64.  Greszczuk,  L.  B.  , FAILURE  CRITERIA  FOR  THREE-DIMENSIONAL 
ORTHOTROPIC  SOLIDS,  Douglas  Aircraft  Company  Report  DAC 
60869,  October  1967. 
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SECTION  III 
THEORY  APPLICATION 


The  analysis  methods  that  were  developed  as  described  in  Section  II 
were  applied  to  study  the  impact  response  and  failure  modes  of  graphite  - 
epoxy  composites  as  influenced  by  fiber  and  matrix  properties,  fiber 
orientation,  stacking  sequence,  plate  thickness,  curvature  effects, 
and  material  combinations  (hybrid  composites  and  interleaved  mate- 
rials). In  conjunction  with  this  work,  data  reported  in  the  literature 
was  compiled  on  properties  of  various  fibers  and  composites,  and 
candidate  composites  were  selected  for  detailed  analytical  studies  to 
establish  how  the  variables  noted  above  affect  the  impact  response. 

CANDIDATE  REINFORCEMENT.  MATRIX  AND  COMPOSITE  MATERIALS 
AND  THEIR  PROPERTIES 

There  exist  numerous  types  of  graphite  fibers  and  other  kinds  of 
reinforcement  materials  for  use  in  composites  and  structural  com- 
ponents made  of  composites.  The  various  reinforcing  fibers,  including 
the  numerous  graphite  and  carbon  fibers  as  well  as  their  properties  and 
other  pertinent  data,  are  listed  in  Table  III.  Table  IV  shows  the  prop- 
erties of  composites  made  with  various  types  of  fibers.  Finally,  Table  V 
shows  the  properties  of  composites  made  with  Thornel  300  graphite  fibers 
and  different  resins.  The  properties  shown  in  Tables  HI  and  IV  were 
obtained  from  various  references  given  in  the  literature  and  data  provided 
by  the  fiber  and  composite  manufacturers  as  well  as  from  MDC  in-house 
programs  on  composites. 

Candidate  Fibers  and  Composites  for  Analytical  Studies 

Because  of  the  profusion  of  graphite  fibers  and  even  more  so  of  com- 
posites, it  was  found  necessary  to  concentrate  on  the  following  three 
main  groups  of  composites: 

1.  Composites  having  ultra-high  Young's  modulus  (such  as 
Thornel  75-epoxy  and  Celion  GY70-epoxy). 

2.  Composites  having  high  modulus  and  moderate  strength  (such 
as  HMS  graphite -epoxy  and  Thornel  50s-epoxy). 

3.  Composites  having  high  strength  and  moderate  modulus  (such 
as  HTS  graphite -epoxy  and  Thornel  300-epoxy). 

Table  VI  gives  the  unidirectional  properties  of  two  graphite -epoxy 
composites  belonging  to  each  group.  For  later  use,  properties  of 
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Table  IV 

PROPERTIES  OF  FIBERS  AND  COMPOSITES 
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Table  V 

PROPERTIES  OF  THORNEL  300  GRAPHITE-FIBER  COMPOSITES 
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(1)  Inhouse  test  results; 

(2)  Test  data  from  Inhouse  programs  and  from  Narmco 

(3)  Normalized  to  fiber  content  of  k « 68% 

(4)  SB  denotes  sandwich  beam  test;  c denotes  coupon  test 
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Table  VI 

PROPERTIES  OF  THE  CANDIDATE  GROUPS  OF  COMPOSITES 
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Calculated  or  estimated  frem  data  for  other  composites  for  which  results  are  shown  In  the  table. 
t.t  assumed  the  same  as  t. 


unidirectional  glass -epoxy  composites  are  also  given  therein.  The 
results  presented  in  Table  VI  are  experimental  values,  except  as 
noted  therein.  Using  the  properties  shown  in  Table  VI,  calculations 
were  made  for  the  properties  of  pseudo-isotropic  laminates  for  use  in 
the  SAAS  computer  code.  Table  VII  shows  the  final  results,  while 
Figure  22  shows  a comparison  of  predicted  tensile  strength  and  Young's 
modulus  of  pseudo-isotropic  composites  with  limited  data  found  in  the 
literature. 

Inasmuch  as  the  properties  of  the  two  composites  within  any  of  the  three 
groups  are  approximately  the  same,*  the  seven  candidate  materials 
were  narrowed  down  to  three  graphite -epoxy  composites  and  glass- 
epoxy.  The  properties  required  for  use  with  SAAS  III  of  the  four  final 
materials  are  given  in  Table  VIII.  The  properties  used  in  ASAAS  were 
those  shown  in  Table  VI. 


Candidate  Resin  Materials 

Most  of  the  data  given  in  the  literature  is  for  composites  made  with 
epoxy  resins  having  Young's  modulus  of  between  0.  5 x 10&  to  0.  6 x 10*5 
psi.  No  data  were  found  for  composites  made  with  low -modulus 
(Er  * 0.  10  x 1C)6  psi)  or  high-modulus  (Er  ~ 1 x 10°  psi)  resins.  Conse- 
quently, to  establish  how  properties  of  resin  affect  impact  behavior  of 
composites,  it  was  found  necessary  to  make  estimates  for  the  various 
properties  of  low-  and  high-modulus  resins.  This  was  done  by  com- 
piling test  data  for  various  epoxy  resins,  plotting  tensile  strength  of 
resin  versus  its  modulus  in  tension,  compressive  strength  of  resin 
versus  its  compressive  modulus,  and  compressive  modulus  versus 
tensile  modulus,  and  extrapolating  the  curves  to  estimate  the  various 
properties  of  the  resins  with  high  and  low  Young's  modulus.  Figure  23 
shows  a plot  of  the  compressive  strength  of  various  resins  versus 
Young's  modulus  in  compression.  The  estimated  properties  of  the 
high-,  low-,  and  ultra -low -modulus  resins  are  given  in  Table  IX, 
together  with  data  for  a typical  epoxy  resin  (Mid  E)  used  in  composites 
for  v/hich  data  are  given  in  the  literature  and  shown  in  Tables  VI,  VII, 
and  VIII.  Except  for  the  resin  designated  as  ULE,  the  strength 
properties  of  the  remaining  resins  were  estimated  as  described  above. 
The  strength  of  the  ULE  matrix  was  arbitrarily  assumed  to  be  high. 
Using  the  properties  given  in  Table  IX,  calculations  v/ere  made  for  the 
properties  of  composites  made  with  Thornel  300-graphite  fibers  and 
the  various  matrix  materials.  Table  X shows  the  final  results.  The 
properties  of  the  Thornel  300  fibers  used  in  these  calculations  were  as 
follows: 


Young's  modulus  of  the  fiber  in  the  fiber  direction, 
E^  = 34  x 10°  psi 


:;:Except  for  the  difference  in  tensile  strength  between  Thornel  753- 
epoxy  and  Celion  GY70-epoxy. 


70 


PSI,  Ex  10‘°  PSI  (THEORY) 
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Figure  22.  Test-Theory  Comparison  of  Tensile  Strength  and  Young's 
Modulus  of  Pseudo-Isotropic  Composites  Made  of  Various 
Fiber-Matrix  Combinations. 
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SHEAR  STRENGTH,  <t  KT  psi  8.6  13.7  17.6 
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♦FIBER  VOLUME  FRACTION  OF  ALL  COMPOSITES  WAS  TAKEN  AS  kf  = 0.66. 


Young's  modulus  of  the  fiber  in  the  transverse  direction, 

EfT  = 3.  1 x 106  psi 

Shear  modulus  of  the  fiber,  = 1.  5 x 10^*  psi 

Poisson's  ratio  of  the  fiber,  = 0. 2 

Tensile  strength  of  the  fiber,  cr^  = 325  x 10^  psi 

3 

Compressive  strength  of  the  fiber,  0^  = 325  x 10  psi 

The  elastic  properties  of  unidirectional  composites  (El,  Ex,  Glt> 
and  vlt)  were  calculated  using  equations  given  in  Reference  65.  The 
tensile  and  compressive  strengths  of  the  composites  in  the  fiber 
direction  were  calculated  from  the  law  of  mixtures.  Transverse 
tensile  strength  of  the  composites,  o"xt*  was  calculated  from 
equations  given  in  Reference  66.  In  calculating  crxt»  it  was  assumed 
that  the  volume  fraction  of  voids  within  the  composites  was  kv  = 2%. 

The  compressive  strength  of  the  composites  in  the  transverse  direction 
was  estimated  from  the  results  shown  in  Figure  24,  which  shows  an 
empirical  relationship  between  the  compressive  Young's  modulus  in 
the  transverse  direction  and  the  transverse  compressive  strength  of 
the  composite.  The  results  shown  in  Figure  24  are  similar  to  the 
strength -modulus  relationship  for  the  resins,  which  is  shown  in 
Figure  23.  Finally  the  interlaminar  shear  strength,  tlt»  was 
estimated  from  the  results  given  in  Reference  67,  where  it  is  shown 
that  a relationship  exists  between  <rxt  and  Tlx* 

Using  the  properties  of  unidirectional  composites  given  in  Table  X,  the 
various  properties  of  pseudo-isotropic  laminates  were  then  calculated. 
The  final  results  ar  _ shown  in  Table  XI.  These  results  are  analogous 
to  the  results  given  in  Table  VIII  except  that  the  variable  here  is  the 
matrix  material,  whereas  in  Table  VIII  the  variable  was  the  fiber 
material. 

INFLUENCE  OF  MATERIAL  PROPERTIES  ON  IMPACT  PARAMETERS 
(SURFACE  PRESSURE.  AREA  OF  CONTACT.  AND  IMPACT  DURATION 
VERSUS  IMPACT  VELOCITY) 

The  analyses  described  in  Section  II  were  applied  to  predict  the  influence 
of  impact  velocity  on  the  surface  pressure,  area  of  contact,  and  impact 

65.  Greszczuk,  L.  B. , THEORETICAL  AND  EXPERIMENTAL  STUDIES 
ON  PROPERTIES  AND  BEHAVIOR  OF  FILAMENTARY  COMPOSITES, 
Proceedings  of  the  2 1st  Annual  Technical  Conference  of  the  Society 
of  Plastics  Industries,  Chicago,  Illinois,  February  1966. 

66.  Greszczuk,  L.  B.  , MICROMECHANICS  FAILURE  CRITERIA  FOR 
COMPOSITES,  Final  Summary  Report  prepared  under  Naval  Air 
Systems  Command  Contract  No.  N00019-72-0221,  May  1973. 

67.  Greszczuk,  L.  B. , MECHANICS  OF  FAILURE  OF  COMPOSITES, 
Final  Summary  Report  prepared  under  Naval  Air  Systems  Command 
Contract  No.  N00019-73-0405,  May  1974. 
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Figure  24.  Correlation  Between  Transverse  Young's  Modulus  and 
Transverse  Compressive  Strength  of  Graphite -Epoxy 
Composites. 
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♦MATERIAL  ASSUMED  TO  HAVE  PSEUDO-ISOTROPIC  LAYUP  IN  THE  PLANE  OF  THE  TARGET.  HF- 
PLANE  TARGET  PROPERTIES  DESIGNATED  BY  r;  PROPERTIES  OF  THE  TARGET  IN  THE  THICKNESS 
DIRECTION  DESIGNATED  BY  z.  IMPACT  DIRECTION  COINCIDES  WITH  z-DIRECTION. 
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duration  for  the  four  materials  described  in  Table  VIII.  In  calculating 
the  required  relationships,  the  target  thickness  in  the  impact  direction 
was  assumed  to  be  large  compared  to  the  impactor  dimensions;  fiber 
layup  in  all  four  composites  was  assumed  to  be  pseudo-isotropic,  and 
the  impactor  was  assumed  to  be  a steel  sphere.  The  properties  of  the 
impactor  were  taken  as  (see  Section  II  for  definition  of  terms): 


m 

= i ; 

Ej  = 30  x 10  psi 

n 

= l ; 

Uj  = 0.  33 

s 

= 2 ; 

pj  = 0.288  lb/in.  3 

R.  = R,, , = 0.  75  in. 
lm  1M 

The  target  was  assumed  to  be  flat  so  that  R2m  = R-2M  • and  its  mass 
was  assumed  to  be  large  compared  to  impactor  mass.  Using  the  equa- 
tions given  in  Section  2 and  the  properties  given  above  and  in  Table  VIII, 
calculations  were  made  of  k2's  tor  the  various  materials.  These  in  turn 
were  used  to  calculate  the  surface  pressure,  area  of  contact,  and  impact 
duration  as  a function  of  impact  velocity  for  the  four  materials  of  interest. 
The  values  of  k2's  (for  the  targets)  are  shown  in  Table  XII.  The  latter 
also  shows  values  of  ki  for  the  impactor  and  kj  + k2  values.  The  varia- 
tion of  tlie  maximum  surface  pressure,  qQ,  radius  of  the  area  of  contact, 
a,  and  duration  of  contact,  t0,  as  a function  of  impact  velocity  is  shown  in 
Figures  25,  26,  and  27,  respectively.  Both  the  radius  of  the  area  of 
contact,  a,  and  the  impact  duration,  tD,  have  been  normalized  with 
respect  to  R,  the  radius  of  the  spherical  impactor.  Thus  the  results 


Table  XU 

VALUES  OF  k2  FOR  DIFFERENT  COMPOSITES 
(BASED  ON  PROPERTIES  DATA  SHOWN  IN  TABJ  E VIII) 


MATERIAL 

MATERIAL 

DESIGNATION 

TARGET 
k2  x TO7 
(in?/lb) 

IMPACTOR 
k1  x 107 
(infv’lb) 

(k-j+kg)  x 107 
(1n?/lb) 

CELION  GY-70- 
EPOXY 

2 

2.324 

0.095 

2.419 

HMS-EP0XY 

3 

1.726 

0.095 

1.821 

THORNEL  300- 
EPOXY 

6 

1.541 

0.095 

1.636 

S-GLASS-EPOXY 

7 

1.177 

0.095 

— 

1.272 
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IMPACT  VELOCITY  (IN. /SEC  ) 


SURFACE  PRESSURE,  qQ  (ksl) 

Figure  25.  Maximum  Surface  Pressure  Versus  Impact  Velocity 


Contact  Duration  Between  Taiget  and  Impactor 
Versus  Impact  Velocity. 
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presented  in  Figures  25,  26,  and  27  apply  to  steel  spherical  impactors 
of  arbitrary  radii.  Moreover,  as  can  readily  be  shown  from  equations 
given  in  Section  II,  for  a spherical  impactor  impacting  a fiat  target, 
qQ  ^ f(Rl),  so  that  the  curves  shown  in  Figure  25  also  apply  to  spherical 
impactors  of  arbitrary  radii,  provided  the  latter  are  made  of  ste°l. 


INFLUENCE  OF  FIBER  PROPERTIES  ON  IMPACT  RESPONSE 

With  the  relationship  between  qQ,  a,  and  v knjwn,  a nondimensional 
solution  was  obtained  (using  SAAS  III  computer  program  described  in 
Section  Hand  the  grid  model  shown  in  Figure  28)  for  the  internal  stresses 
in  targets  made  of  various  composite  materials  described  in  Table  VIII. 

To  determine  which  composite  is  most  resistant  to  impact  and  to  minimize 
the  number  of  variables,  the  layup  in  all  composite  materials  was 
assumed  to  be  pseudo-isotropic.  Moreover,  the  thicKnees  of  the  target 
in  the  direction  of  impact  was  assumed  to  be  large  compared  to  the  radius 
of  the  impactor,  so  that  only  local  stresses  and  deformations  (in  the 
region  of  impact)  had  to  be  considered.  For  a given  qD  and  a,  the  com- 
puter output  consisted  of  radial,  circumferential,  normal,  and  shear 
stresses  at  various  points  within  the  target;  deformations  at  various 
points  within  the  targets;  principal  stresses  at  various  points  within  the 
targets;  and  equivalent  stress  contours  as  obtained  by  applying  the  failure 
criterion  at  various  points  within  the  target.  For  a given  q0  and  a,  the 
failure  envelope  was  also  determined.  Figures  29  and  30  shew  the  com- 
puter printout  of  the  equivalent  stress,  cre,  contours  for  Materials  2 
and  6.  Those  were  used  to  define  the  failure  zone  for  any  given  impact 
condition.  Failure  occurs  at  points  where  <re  ^ 1. 


The  pertinent  results  for  the  impact  response  of  the  three  graphite - 
epoxy  composites  are  presented  in  Tables  XIII  and  XIV.  Table  XIII 
shows  the  pertinent  data  for  the  case  when  qQ  is  kept  constant,  the  lat'er 
being  equal  to  the  maximum  surface  pressure  resulting  from  an  impact 
of  a 0.75-in.  -radius  steel  spherical  impactor.  Table  XIV  shows  the 
results  for  a constant  impact  velocity.  Of  the  three  materials  for  which 
results  are  presented,  Material  6 (Thornel  300-epoxy)  has  minimum  Y. 
Examination  of  results  obtained  to  date  shows  that  qQ  is  extremely  sen- 
sitive to  Ez,  the  Youngs  modulus  in  the  direction  of  impact.  By 
decreasing  Ez,  qQ  and  Y will  also  decrease.  Moreoever,  since  material 
exhibits  yield  behavior  in  compression  it  appears  desirable  to  tafct:  this 
into  account  in  future  work  (through  reduced  Ez)  when  calculating  qc,  and  a. 
Howevei  , to  establish  the  relative  rating  of  different  composites,  the 
assumption  of  elastic  behavior  appears  sufficient  at  this  time.  On  the 
basis  of  results  presented  in  Tables  XIII  and  XIV,  the  ranking  of  the 
three  materials  for  response  to  impact  damage  is  (starting  with  the  best 
material) 


1.  Thornel  300-Epoxy 

2.  HMS-Epoxy 

3.  GY70 -Epoxy 
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Figure  28.  Grid  Used  in  the  Computer  Solution. 
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Figure  30, 


Equivalent  Stress  Contours  for  Thornel  300- Epoxy 

{Material  No.  6)  Subjected  to  q = 100  ksi 
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Table  XIII 

DAMAGE  ZONE  FOR  A CONSTANT  SURFACE  PRESSURE 

(qQ  = 100  KSI) 


MATERIAL 

ki+kgXlO^ 

(in.2/lb) 

CONTACT 
RADIUS,  a 
(IN. ) 

IMPACT 

VELOCITY 

V 

(IN. /SEC) 

h 

a 

h 

a 

Y (IN.) 

i 

2 (GY70-E) 

2.42 

0.089 

105 

1.70 

1.87 

0.158 

3 (HMS-E) 

1.32 

0.067 

59 

1.1*0 

1.1*8 

0.096 

6 (T300-E) 

1.63 

0.06l 

1*8 

1.26 

1.26 

0.08l 

; Yg  = Depth  of  damage  zone 

! Y_  = Half  the  vidth  of  the  damage  zone. 

I K 

i Y = y’yTyT'  = Equivalent  radius  of  the  damage  zone. 

h K 


Table  XIV 

DAMAGE  ZONE  FORA  CONSTANT  IMPACT  VELOCITY 

(v  = 100  IN./SEC) 


MATERIAL 

q xlO-3 
0 

PSI 

a 

(IN.) 

h 

a 

h 

a 

Y 

(IN.) 

2 (GY70-E) 

98 

0,088 

2.11 

2.21 

0.190 

3 (HM3-E) 

12l* 

0.083 

2.06 

1.94 

0.166 

6 (T300-E) 

134 

0.081 

2.22 

1.74 

0.158 

— . -J 
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In  addition  to  the  above,  studies  were  also  conducted  on  local  failures 
resulting  from  impact.  This  included  determination  of  relative  velocities 
to  cause  subsurface  shear  failure  and  its  location,  and  determination  of 
the  velocity  to  cause  tensile  failure  at  the  edge  of  the  area  of  contact 
where  the  tensile  stresses  are  maximum.  The  results  for  the  local 
failures  are  shown  in  Tables  XV  and  XVI,  and  demonstrate  that 
Thornel  300-epoxy  is  the  best  material.  Moreover,  they  show  that  the 
failure  is  initiated  by  subsurface  shear. 

INFLUENCE  OF  MATRIX  PROPERTIES  ON  IMPACT  RESPONSE  OF 
COMPOSITES  REINFORCED  WITH  THORNEL  300  GRAPHITE  FIBERS 

To  establish  if  the  response  to  impact  loading  of  composites  made  with 
Thornel  300  could  be  further  improved  by  changing  the  matrix  properties, 
studies  similar  to  those  described  in  the  previous  section  were  conducted 
on  composites  made  with  resins  described  in  Table  IX.  The  properties 
of  the  latter  composites,  calculated  as  described  in  Section  3,  are  given 
in  Tables  X and  XI. 

The  damage  zones  in  composites  made  with  different  resins  and  subjected 
to  a 100-in./sec  impact  by  a 1.  5 -in.  -diameter  steel  sphere  are  shown  in 
Figure  31.  The  results  shown  there  are  for  a semi-infinite  target.  The 
cross-sectional  areas  of  the  damage  zones  in  the  various  materials  are 
given  in  Table  XVII. 

Thus,  of  the  first  four  materials  listed  in  Table  XVII,  the  damage  zone 
in  material  T300/ME  appears  to  be  the  lowest,  although  the  damage  zones 
in  materials  T300/ULE  and  T300/HE  do  not  differ  significantly  from  the 
lowest  value.  Material  T300/LE  has  the  largest  damage  zone.  Its  shape 
indicates  that  the  low  shear  strength  of  T300/LE  is  the  parameter  that 
governs  the  size  and  shape  of  the  damage  zone  (see  Table  XI).  The  fact 
that  the  damage  zone  in  T300/ULE  is  approximately  the  same  as  in 
T300/HF  indicates  that  the  resin  strength  has  a larger  influence  than 
the  resin  modulus.  To  verify  this  conclusion,  the  impact  response  of 
one  more  material  was  investigated:  a material  with  elastic  properties 
corresponding  to  those  of  T300/TTLE  and  strength  properties  corre- 
sponding to  those  of  T300/ME,  lhe  damage  zone  in  the  latter  material, 
designated  in  Table  XVII  as  conceptual  material  T300/ULE  - T300/ME, 
was  significantly  lower  than  that  of  the  other  materials  shown  therein. 
This  result  further  confirms  the  conclusion  that  resin  strength  rather 
than  modulus  has  the  dominating  influence  on  the  size  of  damage  zone. 

INFLUENCE  OF  TARGET  THICKNESS  ON  IMPACT  RESPONSE  OF 
COMPOSITES 

The  composite  material  selected  for  establishing  how  impact  response 
and  failure  modes  are  influenced  by  thickness  of  the  composite  target 
was  Thornel  300 -epoxy.  The  properties  of  this  composite  are  given  in 
Table  VIII,  wherein  the  material  is  designated  as  Material  6.  Inasmuch 
as  target  thickness  was  the  only  variable  of  interest,  the  composite  was 
assumed  to  have  a pseudo-isotropic  layup  in  the  plane  of  the  target. 
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Table  XV 


RELATIVE  IMF  ACT  VELOCITIES  TO  INITIATE  SHEAR  FAILURE 
AT  A POINT  WITHIN  THE  TARGET 


MATERIAL 

MAXIMUM 
SHEAR  STRESS 

XRZ^o 

LOCATION  OP 
MAXIMUM  ^ 

RELATIVE  IMPACT  VELOCITY 
TO  INITIATE  FAILURE* 

R 

a 

Z 

a 

2 (GY70-E) 

0.35^ 

0.90 

0.175 

1 

3 (HMS-E) 

0.326 

0.90 

0.175 

2.2 

6 (T300-E) 

0.288 

0.90 

0.225 

4.5 

•Normalized  with  respect  to  v2g,  the  velocity  to  cause  shear  failure  in 
material  #2. 

Table  XVI 

RELATIVE  IMPACT  VELOCITY  TO  INITIATE  TENSILE  FAILURES 

<atI=  1) 


I 

! 

i 

vi 

MATERIAL 

MAXIMUM  O 
TENSILE  ~ 
STRESS,  ^0 

RELATIVE  IMPACT 
VELOCITY  TO  INITIATE 
FAILURE* 

1 

i 

2 (GY70-E) 

0.86 

3.2 

j 

3 (HMS-E) 

0.63 

10.7 

i 

i 

6 (T300-E) 

0.45 

59.0 

•Normalized  with  respect  to  v2g. 

Table  XVII 


INFLUENCE  OF  RESIN  ON  SIZE  OF  DAMAGE  ZONE 
IN  COMPOSITES  MADE  WITH  THORNEL  300  FIBERS 


1 

MATERIAL 

CROSS-SECTIONAL  AREA 
OF  DAMAGE  ZONE 
(IN. 2) 

T300/ULE 

0.0133 

T300/LE 

0.0325 

T300/ME 

0.0119 

T300/HE 

0.0122 

Conceptual  Material 
T300/ULE  - T300/ME 

0.0057 
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(b)  T300/LE 

Figure  31.  Damage  Resulting  From  a 100-In.  /Sec  Impact  of  a 1.  5-In.  -Diameter 
Steel  Sphere  Into  Semi -Infinite  Composite  Targets  Made  of 
Thornel  300  Fibers  and  Different  Resins. 
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The  influence  of  target  mass  on  surface  pressure  resulting  fror  foreign 
object  impact  was  determined  from  Equation  (14).  If,  for  a given  set 
of  impact  conditions,  q£  is  the  maximum  surface  pressure  in  a target  of 
finite  thickness  and  qoo  is  the  maximum  surface  pressure  in  a semi- 
infinite target,  then  it  follows  from  Equation  (14)  that  the  relationship 
between  q£,  q^.  and  masses  of  the  target  and  the  impactor  is 


R 


pi- 


where 


(32) 


m£  is  the  effective  mass  of  the  target  and  mj  is  the  mass  of  the  impactor. 
The  effective  mass  of  the  target,  m2,  will,  in  general,  be  greater  than 
mj;  therefore,  \ > 1.  For  the  extreme  case  of  X.  = 1,  evaluation  of 

Equation  (32)  gives 


— = 0.87 

q« 


As  \ increases,  the  ratio  q£/q«> — ♦ 1.  The  above  results  show  that  use 
of  Equation  (14)  to  predict  surface  pressure  in  targets  of  finite  thickness 
gives  conservative  results.  Consequently,  when  studying  the  influence 
of  target  thickness  on  impact  response,  Equation  (14)  was  used  to 
calculate  the  maximum  surface  pressure. 

In  obtaining  the  numerical  results  on  the  plate  thickness  effects,  it  was 
assumed  that  the  plates  were  circular  in  shape  and  were  clamped  along 
the  outer  boundary.  Moreover,  the  impactor  was  assumed  to  be 
spherical  in  shape  and  its  properties  were  assumed  to  be  the  same  as 
those  used  in  the  numerical  examples  given  previously  (Section  III): 


m 

= 1 ; 

Ej  = 30  x 10^  psi 

n 

= 1 ; 

Uj  = 0.  33 

s 

= 2 ; 

p j = 0.  288  lb/in. 

Rlm 

= R...  = 0.75  in. 
1 M 

Some  typical  results  on  the  influence  of  plate  thickness  and  dimensions 
on  the  impact-induced  radial  stresses  are  shown  in  figures  32  and  33. 
The  results  presented  therein  give  the  distribution  of  radial  stresses 
through  the  thickness  of  the  plate  along  the  line  passing  through  the 
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center  of  the  area  of  contact.  To  make  the  results  of  practical  interest, 
both  the  plate  thickness  and  its  diameter  were  nondimensionalized  with 
respect  to  the  radius  of  the  area  of  contact,  a.  Thus,  H and  D appearing 
in  Figures  32  and  33  represent  the  nondimensionalized  thickness  and 
diameter,  respectively.  Figure  32  shows  the  influence  of  plate 
thickness  on  the  distribution  of  radial  stresses  for  a plate  with  a given  D, 
whereas  the  results  shown  in  Figure  33  are  for  a plate  with  constant  H 
and  varying  D.  For  the  numerical  example  given  in  Figure  32,  plate 
thickness  effect  is  seen  to  be  significant  when  the  radius  of  the  area  of 
contact  is  equal  to  the  plate  thickness,  that  is,  when  H = 1.  If  H > 1, 
the  radial  stresses  appear  to  be  governed  by  local  deformation  of  the 
plate,  whereas  if  H < 1,  the  radial  stresses  appear  to  be  governed  by 
plate  bending  deformations. 


The  influence  of  local  deformations  and  plate  bending  deformations  is 
further  illustrated  in  Figures  34a  and  35.  The  deformations  at  the  edge 
of  the  art.?  of  contact  have  been  normalized  with  respect  to  the  maximum 
deformation  at  the  center  of  the  area  of  contact.  Moreover,  as  before, 
both  the  plate  radius,  R,  and  the  plate  thickness,  H,  have  been  non- 
dimensionalized with  respect  to  the  radius  of  the  area  of  contact,  a. 

As  is  readily  apparent  from  the  results  presented  in  Figures  32-35, 
the  dimensions  of  the  impacted  plate  (thickness  and  radius)  have  a 
significant  influence  on  the  stress  state  within  the  plate  and  on  the  plate 
deformation  resulting  from  a foreign  object  impact.  To  establish  how 
the  impact -induced  failure  modes  are  influenced  by  the  dimensions  of 
the  composite  target,  SAAS  III  finite -element  computer  code  was  used. 
The  pertinent  results  for  the  influence  of  plate  thickness  and  plate  radius 
on  the  failure  modes  are  shown  in  Figures  36  and  37.  Figure  36  shows 
the  maximum  critical  surface  pressure,  qjj,  at  which  failure  initiates 
at  a point  in  the  plate.  The  results  plotted  there  are  for  a constant  R 
and  varying  H,  the  nondimensional  plate  radius  and  plate  thickness, 
respectively.  The  results  presented  in  Figure  36  show  that  for 
H <1.5,  the  failure  is  governed  by  bending  stresses  and  initiates  on 
the  bottom  surface  of  the  plate.  For  H £ 1.5,  the  failure  is  governed 
by  local  contact  stresses  and  occurs  on  the  top  surface  of  the  plate. 

As  H increases,  q*  approaches  a constant  value,  indicating  that  the 
failure  consists  of  local  damage.  Figure  37  shows  results  similar  to 
those  presented  in  Figure  36  except  that  the  variable  is  the  nondimen- 
sionalized plate  radius  R. 


As  was  discussed  in  Section  II,  for  a given  set  of  impact  parameters, 
both  the  surface  pressure  and  the  area  of  contact  over  which  it  acts  are 
time  dependent.  For  the  case  of  a semi-infinite  composite  target,  the 
variation  of  the  maximum  surface  pressure  with  time  is  shown  in 
Figure  38.  The  extent  of  damage  in  such  a target  at  pressure  or  time 

intervals  denoted  in  Figure  38  as  a,  b,  c, £,  is  shown  in  Figure  39. 

The  points  in  Figure  38  denoted  by  a,  b,  c f,  correspond  to 

Sections  (a),  (b),  (c) (f)  of  Figure  39.  All  the  results  shown  in 

Figure  39  are  drawn  to  scale.  This  includes  relative  magnitude  of  the 
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Figure  34.  Influence  of  Plate  Thickness  on  Local  and  Plate 
Bending  Deformations  at  the  Edge  of  the  Area  of 
Contact 
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Figure  35.  Influence  of  Plate  Radius  on  Local  and  Plate  Bending 
Deformations  at  the  Edge  of  the  Area  of  Contact. 
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Figure  36.  Influence  of  Plate  ThicKness  on  Maximum  Surface 
Pressure  at  Which  Failure  Initiates  at  a Point  in 
a Plate. 
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Figure  37.  Influence  of  Plate  Radius  on  Maximum  Surface 
Pressure  at  Which  Failure  Initiates  at  a Point 
in  a Plate. 


IMPACTOR  RADIUS,  R)m  = Rlf.  = 0.75  IN. 

IMPACTOR  MATERIAL:  STEEI.{E1  = 30  x TO6  P$I,v^=  0.33,w1  = 0.50P  LB. 

p1  = 0.283  LB/ IN?,  m1«1 -31x10“3LB-SCC2/  I'l. ) 

TARGET  MATERIAL:  THORNEL  300/EP0XY 
TARGET  MASS:  SEMI-INFINITE  (m2  » ^ ) 

IMPACT  VELOCITY,  Vfl  = 1500  IN./SEC 
IMPACT  DURATION,  tQ  = 152y  SEC 

MAXIMUM  AREA  OF  CONTRACT,  aQ  = 0.242  IN.  (AT  t=0.5  t ) 


0 20  40  r.  0 CO 


TIME,  t,  y SEC. 

Figure  38.  Increase  of  Maximum  Surface  Pressure  With  Time. 
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surface  pressure,  areas  over  which  it  acts,  and  the  size  of  the  failure 
zone.  The  results  shown  in  Figure  39  can  also  be  used  for  determining 
the  damage  zone  if  vQ  < 1500  in.  /sec.  For  example,  if  the  impact 
velocity  is  v0  = 270  in.  /sec,  then  from  the  curves  given  in  Figure  25, 
qQ  = 200  ksi.  The  failure  or  damage  zone  for  this  impact  condition  is 
given  in  Figure  39(e). 

The  increase  in  the  damage  zone  in  plates  of  finite  thickness  as  a func- 
tion of  increasing  surface  pressure  (or  correspondingly  increasing 
impact  velocity)  is  shewn  in  Figures  40  and  41.  Figure  40  shows  the 
growth  of  the  damage  zone  for  the  case  when  the  thickness  and  radius  of 
the  plate  are  kept  constant  and  the  impact  velocity  is  increased.  Fig- 
ure 41  shows  the  growth  of  the  damage  zone  (shaded  area)  for  the  case 
where  R,  H,  and  q vary,  but  R/H  remains  constant. 

The  results  presented  in  Figures  39,  40,  and  41  further  illustrate  how 
thickness  of  the  target  and  its  dimensions  influence  the  damage  induced 
by  foreign  objects  impacting  targets  made  of  composite  materials. 

INFLUENCE  OF  PLATE  BOUNDARY  CONDITIONS  ON  IMPACT 
RESPONSE 


The  results  obtained  in  Section  III  were  for  circular  composite  plates 
clamped  along  the  outer  periphery.  Two  other  cases  of  plate  support 
were  investigated:  (1)  plates  with  roller  supports  along  the  outer 
boundary  and  (2)  plates  with  hinge  supports  along  the  outer  boundary. 

In  all  cases,  the  plate  material  (Thornel  300 -epoxy),  fiber  layup  (inplane 
pseudo-isotropic),  and  plate  dimensions  (1.5-in.  radius;  0.  125-in. 
thickness)  were  kept  constant.  Impactor  was  assumed  to  be  a 1.  5-in.  - 
diameter  steel  sphere,  whereas  the  impact  velocities  investigated  were 
v = 30  in.  /sec  and  v = 50  in.  /sec.  At  an  impact  velocity  of  30  in.  /sec, 
the  damage  zone  was  localized  near  the  impact  surface.  The  width,  w, 
and  the  depth,  d,  of  the  damage  zone  for  plates  with  various  edge 
supports  are  shown  in  Table  XVIII.  The  results  on  the  dimensions  of 
damage  zones  resulting  from  50-in.  /sec  impact  are  presented  in 
Table  XIX.  In  the  latter  case,  the  damage  zone  extended  through  the 
thickness  of  the  plate  for  all  three  types  of  plate  boundary  supports. 
Figure  42  shows  the  dimensions  of  the  damage  zones  (resulting  from  a 
50 -in.  / sec  impact  by  a 1.  5 -in.  -diameter  steel  sphere)  ir  composite 
plates  with  three  types  of  boundary  conditions. 


From  the  results  presented  above  it  is  seen  that  for  impact  velocities 
which  cause  local  damage,  the  end  fixity  of  the  plate  has  a negligible 
influence  on  the  size  of  the  damage  zone.  On  the  other  hand,  as  the 
impact  velocity  is  increased  so  that  damage  occurs  on  top  and  bottom 
surfaces  of  the  plate,  the  plate  edge  fixity  has  a significant  influence 
on  the  width  of  the  damage  zone.  For  an  impact  velocity  of  50  in.  /sec, 
the  width  of  damage  zone  was  smallest  in  plates  with  clamped  edges 
and  largest  in  plates  with  roller -supported  edges. 
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(or  Increasing  Impact  Velocity). 


Figure  41.  Increase  in  Damage  Zone  as  a Function  of  Increasing 
Surface  Pressure  and/or  Increasing  Impact  Velocity 
(R/H  = Constant). 
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Table  XVIH 

INFLUENCE  OF  PLATE  BOUNDARY  CONDITIONS  ON  DIMENSIONS 
OF  DAMAGE  ZONE  CAUSED  BY  30  IN.  /SEC  IMPACT 
BY  A 1.  5 -IN. -DIAMETER  STEEL  SPHERE 


EDGE  SUPPORT  OF  THE  PLATE 

w(in. ) 

d(in. ) 

CLAMPED  EDGES 

0.120 

0.050 

ROLLER  SUPPORTED  EDGES 

0.119 

i 

0.01+1+  ; 

HINGE  SUPPORTED  EDGES 

0.119 

0.01+5 

Table  XDC 

INFLUENCE  OF  PLATE  BOUNDARY  CONDITIONS  ON  DIMENSIONS 
OF  DAMAGE  ZONE  CAUSED  BY  50  IN.  /SEC  IMPACT 
BY  A 1.  5-IN.  -DIAMETER  STEEL  IMPACTOR 


EDGE  SUPPORT  OF  THE  PLATE 

WIDTH  OF  DAMAGE  ZONE  (IN.) 

TOP  SURFACE* 

BOTTOM  SURFACE 

CLAMPED  EDGES 

0.217 

0.231 

ROLLER  SUPPORTED  EDGES 

0.1+22 

0.1+60 

HINGE  SUPPORTED  EDGES 

0.305 

0.351* 

^Impacted  surface 


(c)  ROLLER  SUPPORTED 


Figure  42.  Influence  of  Plate  Boundary  Conditions  on  Impact 
Response  ( v=  50  In.  /Sec,  q = 102  ksi). 


INFLUENCE  OF  FIBER  ORIENTATION  AND  LAMINATE  STACKING 
SEQUENCE  ON  IMPACT  RESPONSE 


The  material  properties  used  in  studying  the  influence  of  fiber 
orientation  and  laminate  stacking  sequence  were  taken  to  be  those 
corresponding  to  Thornel  300-epoxy,  which  is  designated  as  Material  6 
in  Table  VIII.  The  fiber  orientations  investigated  included  unidirectional 
composites,  bidirectional  composites  having  0°/90°  or  ± 45°  dispersed 
fiber  layup,  and  tridirectional  composites  having  0°/±60°  fiber  layup 
(inplane  pseudo-isotropic).  In  all  cases  the  target  was  assumed  to  be 
a circular  plate,  0.  125  in.  thick  and  having  a 3-in.  diameter.  The 
plate  was  assun  ed  to  be  fixed  along  its  periphery.  As  before,  the 
impactor  was  a 1 . 5 -in.  ••  diameter  steel  sphere.  Because  of  the  shape 
of  the  plate  and  the  fact  that  impact  was  assumed  to  occur  at  the  center 
of  the  plate,  a material  with  ±45°  fiber  layup  is  indistinguishable  from 
a material  with  0°/90°  fiber  layup.  The  pertinent  results  for  the 
influence  of  fiber  orientation  on  the  impact  response  are  presented  in 
Figures  43  through  46.  In  obtaining  the  results  shown  therein  the 
plate  geometry  and  the  impactor  material  size  and  shape  were  kept 
constant.  The  results  shown  in  Figures  43,  44,  and  45  were  cal:ulated 
for  an  impact  velocity  of  30  in.  /sec.  Thus,  the  only  variable  in 
Figures  43  to  45  was  the  fiber  orientation.  The  cross-hatched  area  in 
each  of  the  figures  gives  the  damage  zone  in  the  material  resulting  from 
a 30-in.  /'sec  impact  by  a 1.  5-in.  diameter  steel  sphere.  Figure  46 
shows  the  damage  /.one  in  one  of  the  materials  corresponding  to  a 
50-in.  /sec  impact.  From  the  results  presented  in  Figures  45  and  46, 
it  is  readily  seen  that  if  a plate  is  subjected  to  30 -in.  /sec  impact,  ail 
the  damage  is  within  tne  plate,  whereas  a 50-in.  /sec  impact  causes  the 
damage  to  extend  to  the  outer  surfaces  of  the  plates. 

Of  the  hree  fiber  layups  investigated,  plates  having  only  unidirectional 
fibers  show  the  largest  damage  zone,  wiiereas  the  smallest  damage  zone 
was  found  to  exist  in  plates  with  bidirectional  fiber  layup.  The  damage 
zone  in  plates  with  tridirectional  fiber  layup  is  somewhat  larger  than  it 
is  in  plates  having  bidirectional  fibers. 

To  establish  how  stacking  sequence  influences  the  impact  response,  the 
following  problems  have  been  investigated: 

1.  Bidirectional  composite  with  layers  dispersed  through  the 
thii  kness  so  that  the  bidirectional  composite  could  be  con- 
sidered as  an  equivalent  homogeneous  orthotropic  material. 

2.  Bidirectional  composite  consisting  of  nine  discrete  layers 
having  fiber  orientations  of  0’  and  90°.  Any  one  given  layer 
was  considered  as  homogeneous  orthotropic,  whereas  the 
composite  is  treated  as  multilayer  heterogeneous  orthotropic. 
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igure  43.  Foreign  Object  Impact -Induced  Damage  in  Unidirectional 
Thornel  300- Epoxy  Composite. 


IMPACT  VELOCITY 
V * 30  IN./SEC 


3.  Bidirectional  composite  consisting  of  three  discrete  homo- 
geneous orthotropic  layers,  two  of  these  having  fibers 
oriented  in  the  0°  direction  and  one  layer  having  fibers 
oriented  at  90°  direction. 

The  results  on  the  size  of  the  damage  zone  in  Problem  1 are  shown  in 
Figure  44.  Figures  47  nd  48  show  the  results  for  Problems  2 and  3. 
Comparison  of  results  shown  in  Figures  44  and  47  shows  that  treating 
the  bidirectional  composite  as  an  equivalent  homogeneous  orthotropic 
material  does  not  allow  for  the  complete  determination  of  the  damage 
zone.  Similar  results  would  be  expected  if  the  composite  plate  for 
which  the  results  are  presente  -t  in  Figure  48  were  treated  as  an 
equivalent  orthotropic  material.  Considering  the  plates  as  multilayer 
heterogeneous  orthotropic,  however,  allows  for  a more  accurate 
determination  of  the  damage  zone.  The  influence  of  the  stacking 
sequence  can  be  seen  by  comparing  the  results  shown  ir>  Figures  47 
and  48.  From  the  results  shown  therein  it  is  seen  that  the  damage 
zone  is  reduced  if  the  layers  are  dispersed  through  the  thickness. 


INFLUFNCE  OF  HYBRIDIZATION  (COMBINING  DIFFERENT  COMPOS- 
ITES) ON  IMPACT  RESPONSE 

To  establish  if  impact  response  of  composites  could  be  improved 
through  hybridization,  analytical  studies  were  conducted  on  material 
combinations  and  constructions  shown  in  Figure  51.  Combinations  of 
materials  considered  included:  bidirectional  fiberglass -epoxy  layers 
combined  with  bidirectional  Thornel  300 -epoxy  layers;  and  bidirectional 
PRD  49-epoxy  layers  combined  with  bidirectional  Thornel  300-cooxy 
layers.  In  all  cases  it  was  assumed  that  the  total  plate  thickness 
remained  0.  125  inch.  As  before,  the  plates  were  assumed  to  have 
1.5-in.  radii  to  be  clamped  along  the  outer  boundary,  and  to  be  sub- 
jected to  impact  by  a 1.5-in.  -diameter  steel  impactor.  The  final 
results  for  the  damage  zones  in  fiberglass-epoxy/Thornel  300-epoxy 
and  PRD  40-epoxy/Thornel  300-epoxy  hybrid  composite  plates  sub- 
jected to  30-in.  /sec  and  50-in.  /sec  impact  arc  shown  in  Figures  52,  53, 
54,  and  55.  For  comparison  purposes.  Figures  52  and  53  show  also 
the  damage  zones  in  Thornel  300-epoxy  plates.  Results  trom  studies 
on  fiberglass-epoxy/Thornel  300-epoxy  hybrid  composites  showed 
the  material  Configurations  (a)  and  (b)  to  be  superior  to  material 
Configuration  (c)  (see  Figures  52  and  53);  consequently  only  the  former 
material  configurations  were  considered  in  investigating  the  impai  l 
response  of  hybrid  composites  consisting  of  PRD  49 -epoxy  and  Thornel 
300-epoxy.  All  the  results  given  in  Figures  52  through  55  were  calcu- 
lated using  the  material  properties  given  in  Tables  IV  and  VI.  The 
pertinent  information  on  the  maximum  width  and  depth  of  the  damage 
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NOTE:  SHAl-ED  AREAS  DENOTE  DAMAGE 
ZONES  RESULTING  FROM  IMPACT 


SECTION  O-C 

Induced  Damage  in  a Bidirectional  Thornel 
z a Nine-Layer  Dispersed  Construction. 


Figure  43.  Foreign  Object  Impact -Induced  Damage  in  a Bidirectional  Thornel 
300-Epoxy  Plate  Having  a Three-Layer  Construction. 


j fL  = FAILURE  IN  FIBER  DIRECTION 

fT  = FAILURE  IN  TRANSVERSE  DIRECTION 

j 

NOTE:  ALL  DIMENSIONS  ARE  IN  INCHES 
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Figure  49.  Impact  Damage  on  the  Top  Surface  {Impacted 

Surface)  of  a Nine -Layer  Bidirectional  Thornel 
300- Epoxy  Plate. 


X14 


LOCAL  FAILURE  (TOP)  TENSILE  FAILURE  (BOTTOM) 


Surface  Damage  Caused  by  Impact  of  a 1.  5 -In.  -Diameter 
Steel  Ball. 


^v-'-vtfW-sWv  .istlg#  -«  ■?  a.  - -<•» 





THORNEL  300-EP0XY 

(a)  (b)  (c) 


Figure  51.  Material  Configurations  for  Studying  the  Influence 
of  Hybridization  on  Impact  Response. 

zones  in  various  hybrid  composite  plates  subjected  to  impact  velocities 
at  30  in.  /sec  and  50  in.  /sec  is  summarized  in  Table  XX  and  XXI, 
respectively. 

Because  the  damage  zone  has  a highly  irregular  shape;  no  well-defined 
trends  could  be  established  for  the  influence  of  hybridization  on  impact 
response.  However,  from  the  results  presented  in  Figures  52  through  55 
and  in  Tables  XX  and  XXI  it  is  readily  seen  that  concentrating  fiberglass 
epoxy  at  the  outer  plate  surfaces  makes  the  composite  more  resistant 
to  impact  damage  than  if  the  fiberglass  is  dispersed  through  the  thick- 
ness. If  the  fiberglass  outer  layers  are  replaced  with  PRD  49-epoxy, 
the  damage  zone  grows  in  radius  but  decreases  in  depth.  For  material 
Configurations  (a)  and  (b),  the  increase  in  impact  velocity  causes  a 
greater  increase  in  the  depth  of  the  damage  zone  than  in  its  radius. 

INFLUENCE  OF  CURVATURE  ON  IMPACT  RESPONSE 

In  studying  the  influence  of  target  curvature  on  impact  response,  it  was 
assumed  that  the  targets  are  cylindrical  and  are  made  of  Thornel  300- 
epoxy,  with  67%  of  the  fibers  (five  layers)  oriented  in  the  circumferential 
direction  and  33%  of  the  fibers  (four  layers)  oriented  in  the  axial  direc- 
tion. The  cross  section  of  cylinder  wall,  showing  the  distribution  of 
material  and  thicknesses  of  die  nine  layers,  is  shown  in  Figure  56. 

Radii  of  the  cylinders  were  taken  as  1.5  in.  and  3.  0 in.  , whereas 
their  thicknesses  were  taken  as  t = 0.  125  in.  As  before,  the  impactor 
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RADIAL  DISTANCE,  IN. 
BIDIRECTIONAL  THORNEL  300-EPOXY 


RADIAL  DISTANCE,  IN. 

HYBRID  COMPOSITE 

(a) 


RADIAL  DISTANCE,  IN. 
HYBRID  COMPOSITE 
(b) 


Figure  52.  Damage  Zone  in  3 -In.  -Diamecer  Glass-Epoxy/Thornel 

300- Epoxy  Hybrid  Composite  Plates  Subjected  to  30-In.  / 
Sec  Impact  by  a 1.  5-In.  -Diameter  Steel  Impactor. 
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Table  XX 

INFLUENCE  OF  COMBINING  MATERIALS  ON  THE  SIZE  OF  DAMAGE  ZONE  RESULTING 
FROM  A 30-IN.  /SEC  IMPACT  BY  A 1.5 -IN.  -DIAMETER  STEEL  SPHERE 
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= maximum  width  of  the  damage  zone  in  the  radial  direction 
= maximum  depth  of  the  damage  zone  in  the  thickness  direction 


Table  XXI 
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£s 


AREA  OF  CONTACT 
SECTION  A-A 


2b 


0.125  IN 


FIBERS  ORIENTED  IN 
CIRCUMFERENTIAL  (e)  DIRECTION 


FIBERS  ORIENTED  IN 
AXIAL  (x)  DIRECTION 


0.01668  IN.  (TYP.  FOR  CIRCS.) 
±0.0104  IN.  (TYP.  FOR  AXIAL) 


I 

t 


CYLINDER  WALL 


Figure  56.  Cylinder  Geometry  and  Material  Distribution 
Through  the  Wall  Thickness  of  the  Cylinder. 
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was  assumed  to  be  a 1.  5-in.  -diameter  steel  sphere.  Target  curvature 
was  found  to  influence  the  shape  and  size  of  the  area  of  contact,  maxi- 
mum pressure,  total  load,  and  contact  duration.  The  pertinent  results 
for  the  influence  of  target  curvature  on  the  impact  parameters  are 
presented  in  Table  XXII. 

Using  data  shown  in  Table  XXII  as  an  input  to  the  previously  described 
CYLINDER  code,  calculations  were  made  for  the  dynamic  stresses  in 
the  cylindrical  target,  which  was  assumed  to  have  simply  supported 
ends.  Because  the  low  impact  velocity  resulted  in  a small  area  of 
contact,  it  was  not  possible  to  represent  the  surface  load  and  the  area 
of  contact  as  described  in  Section  II  without  further  extensive  modifica- 
tions of  the  computer  code.  The  best  that  could  be  done  for  the 
problem  considered  here  (without  exceeding  computer  storage  capacity 
or  exceeding  25  minutes  of  computer  time)  was  to  assume  the  load  to 
be  concentrated  on  one  rectangular  element,  the  area  of  which  was 
the  same  as  the  area  of  the  ellipse.  The  input  to  the  code  and  its 
output  were  the  same  as  described  on  Page  55  of  Section  II.  The 
calculated  stresses  were  used  to  determine  the  failure  location, 
failure  mode,  and  failure  sequence  using  maximum  stress  criteria. 
Final  results  of  these  calculations  are  presented  in  Table  XXIII. 

From  the  results  presented  in  Tables  XXII  and  XXIII,  the  following 
conclusions  can  be  made  on  the  influence  of  curvature  on  impact 
response  (assuming  impact  velocity,  target  properties,  and  impactor 
mass,  shape,  and  size  remain  constant): 

1.  In  the  case  of  a cylinder  impacted  by  a sphere,  the  area  of 
contact  is  elliptical  and  approaches  a circle  as  the  radius  of 
the  cylinder  increases  (approaches  flat  surface). 

2.  Area  of  contact  decreases  with  decreasing  cylinder  radius. 

3.  Maximum  load,  resulting  from  impact,  decreases  with 
decreasing  cylinder  raOius. 

4.  Maximum  surface  pressure  increases  with  decreasing 
cylinder  radius. 

5.  Contact  duration  increases  with  decreasing  cylinder  radius. 

6.  Time  required  to  cause  various  failure  modes  decreases  with 
decreasing  cylinder  radius. 

7.  Time  sequence  in  which  different  failure  modes  take  place 
changes  as  the  cylinder  radius  changes. 
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Table  XXIII 

FAILURE  MODE,  LOCATION  AND  SEQUENCE 
IN  3-INCH-  AND  6 -INCH- DIAMETER  CYLINDERS 


FAILURE  MODE 

TIME  AT  WHICH  FAILURE  TAKES  PLACE 
(lisec) 

3-INCH-DIAMETER 

CYLINDER 

6-INCH-DIAMETER 

CYLINDER 

Transverse  Tensile  Failure  on  the 
Inner  Surface 

1U.3 

23.0 

Compression  Failure  in  the  Fiber 
Direction  on  the  Outer  Surface 

61.1 

6U.7 

Compression  Failure  in  Transverse 
Direction  on  the  Outer  Surface 

78.0 

89.0 

Tensile  Failure  in  Fiber  Direction 
on  the  Inner  Surface 

79-7 

85. ^ 

I 


SECTION  IV 


PESO rjPTION  OF  MATERIAL  AND  SPECIMENS  USED 
IN  THE  TEST  PROGRAM 


Following  the  theory  application  phase  of  the  program,  a test  plan  for 
determining  the  threshold  energy  levels  that  cause  fracture  in  graphite  - 
epoxy  composites  subjected  to  impact  was  prepared  and  submitted  for 
approval  to  the  Project  Officer.  Following  its  approval,  the  test  plan  was 
implemented.  Appropriate  graphite -epoxy  specimens  were  then  fabricated 
and  tested  to  experimentally  establish  the  influence  of  the  following  vari- 
ables on  impact  response:  fiber  and  matrix  properties;  fiber  orientation, 
layup,  and  stacking  sequence;  laminate  thickness;  curvature  effects;  and 
hybridization.  Two  types  of  specimens  were  employed:  flat  plates  and 
cylinders.  The  basic  material  (Thornel  300  - epoxy)  used  in  the  majority 
of  test  specimens  was  characterized  for  its  mechanical  properties. 
Moreover,  appropriate  quality  assurance  tests  were  also  performed  on  all 
materials  used  in  the  test  program. 

DESCRIPTION  OF  MATERIALS  USED  IN  TEST  PROGRAM 


Five  different  reinforcing  fibers  and  three  different  resins  were  used  in 
the  various  tasks  of  the  test  program.  These  materials  and  their  proper- 
ties are  described  in  Table  XXIV.  The  properties  data  shown  therein 
were  obtained  from  various  vendors  and  from  literature. 

QUALITY  ASSURANCE  TESTS  ON  INCOMING  MATERIALS  AND 
SPECIMEN  FABRICATION 


Quality  assurance  tests  were  conducted  on  the  various  prepreg  materials 
to  ensure  that  the  material  was  of  acceptable  quality  and  to  verify  and/or 
establish  the  processing  parameters.  The  quality  control  (QC)  tests 
included  flexure  tests  and  short-beam  shear  tests.  Table  XXV  shows 
the  results  from  the  QC  tests,  as  well  as  other  pertinent  information. 
The  final  cure  cycles  that  were  used  in  various  composites  incorporating 
different  resin  systems  were  as  follows: 

1.  Composites  made  with  5208  epoxy 

a.  Enter  cold  autoclave  under  vacuum. 

b.  Hold  vacuum  and  raise  temperature  to  27  5°F  at  5-9°F/min. 


(3)  Material  Batch  32h 


»;e-*-%fi«ft-^-3** 


s^rr?~K . 
. „ ' & : 


c.  Hold  at  275°F  under  vacuum  for  45  minutes. 

d.  Raise  temperature  to  350*F  at  5-9°F/min  and  pressure 
to  100  psi. 

e.  When  pressure  reaches  20  psi.,  vent  bag  to  atmosphere. 

f.  Hold  at  350°F  and  100  psi  for  2 hours. 

g.  Cool  under  pressure  to  below  150°F  and  remove. 

2.  Composites  made  with  ERLA  4617  epoxy 

a.  Enter  cold  autoclave  under  vacuum. 

b.  Hold  vacuum  throughout  cure. 

c.  Raise  temperature  to  175°F  at  5-9°F/min. 

d.  Hold  at  175°F  for  3 hours. 

e.  Raise  temperature  to  200°F  at  5-9°F/min. 

Hold  at  200°F  for  1 hour. 


f. 

g- 


Raise  temperature  to  250°F  at  5-9°F/min  and  pressure 
to  80  psi. 


h.  Hold  at  250°F  and  80  psi  for  2 hours. 

i.  Raise  temperature  to  300°F  at  5-9°F/min. 

j.  Hold  at  300°F  and  80  psi  for  2 hours. 

3.  Composites  made  with  Polysuifone  matrix. 

a.  Enter  hot  press  at  450°F. 

b.  Apply  contact  pressure  for  10-15  min. 

c.  Relieve  pressure  to  allow  volatiles  to  escape. 

d.  Apply  100  psi  for  30  min. 

e.  Remove  from  press  and  cool  under  deadweight  pressure  of 
approximately  2-5  psi. 
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MECHANICAL  PROPERTIES  OF  UNIDIRECTIONAL  THORNEL  300/ 
5208  COMPOSITES 


Because  most  of  the  test  specimens  were  to  be  fabricated  from  Thornel 
300/5208,  appropriate  tests  were  conducted  on  the  material  to  obtain  the 
mechanical  properties  data.  Figure  57  shows  the  geometry  of  the  various 
test  specimens,  whereas  the  typical  failed  specimens  and  the  final  test 
data  are  presented  in  Figure  58  and  Table  XXVI,  respectively.  As  noted 
in  Table  XXVI,  most  of  the  test  values  a-e  an  average  from  three  tests. 
The  test  data  shown  in  Table  XXVI  show  good  correlation  with  test  data 
given  in  the  literature  and  shown  earlier  in  Table  V of  Section  III. 

DESCRIPTION  OF  TEST  SPECIMENS 

The  various  specimens  that  were  used  in  the  experimental  phase  of  the 
program  are  described  in  Table  XXVII.  The  latter  gives  the  panel  or 
cylinder  designation  (for  later  reference),  fiber,  and  matrix  material 
from  which  a given  panel  or  cylinder  was  made,  and  the  fiber  layup. 

Table  XXVIII  relates  the  panels  and  cylinders  to  the  experimental  vari- 
ables that  were  investigated.  A further  description  of  the  various  speci- 
mens that  were  fabricated  is  given  in  Table  XXIX.  The  latter  gives  the 
specimen's  weight,  its  dimensions,  density  of  composite,  fiber  and  void 
content,  and  the  thickness  per  ply.  The  densities  that  are  given  were 
calculated  by  gravimetric  technique  and  also  from  the  weights  and 
dimensions  of  the  specimens.  The  fiber  contents  in  the  various  speci- 
mens were  determined  by  (1)  a nitric  acid  digestion  in  the  case  of 
graphite-fiber-reinforced  composites,  (2)  nitric  acid  digestion  followed 
by  burnout  in  the  case  of  Thornel -epoxy/glass -epoxy  hybrid  composites, 
and  (3)  nitric  acid  digestion  followed  by  sulphuric  acid  digestion  in  the 
case  of  Thornel  300 -epoxy/Kevlar  49-epoxy  hybrid  composites. 

One  problem  encountered  in  fabricating  plates  was  that  the  plate  thick  - 
nesses  could  not  be  kept  exactly  at  0.  125  in.  because  different  prepregs 
and  fiber  layups  give  different  thicknesses  per  ply  of  cured  prepreg  tape. 
Moreover,  to  prevent  warping,  the  fiber  layup  had  to  be  kept  symmetri- 
cal with  respect  to  the  midplane  of  the  plate.  Consequently,  only  the 
total  number  of  plies  of  prepreg  in  each  plate  was  kept  constant,  except 
for  the  plates  in  which  the  plate  thickness  was  a variable. 


■f  15°  e)0‘  O'  0'  10*  0* 

TENSION  FLEXURE  COMPRESSION 


Figure  57.  Test  Specimens  for  Measuring  the  Mechanical  Properties 
of  Thornel  300-5208  Unidirectional  Composite. 
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Table  XXVI 

AVERAGE  MECHANICAL  PROPERTIES  FOR  THORNEL  300-5208 
UNIDIRECTIONAL  COMPOSITES  (MEASURED  DATA) 

DIRECTION 

PROPERTY  VALUE 

1 

Fiber  Content  by  Volume  (3)**  k = 58.5% 

Void  Content  (3)  ky  = 0.85% 

Composite  Density  p = 0.0557  LB/IN.3 

Thickness  per  ply  t = 0.00578  IN. 

AXIAL 

Tensile  Strength  (3)  au  = 196.7  x 103  PSI 
Tensile  Modulus  (4)  Eu  = 19.4  x 106  PSI 
Tensile  Strain  at  Failure  (3)  eLt  = 1.01% 

Compressive  Strength  (6),  oLc  = 203  x 103  PSI 

Flexuru  Strength  (3),  = 210.5  x 103  PSI 

Flexure  Modulus  (2),  ELf  = 20.8  x 10^  PSI 

Interlaminar  Shear  Strength  (3), 


t = 15.89  x 103  PSI 


TRANSVERSE 


45°  (WITH 
RESPECT  TO 
FIBER 
DIRECTION) 


Tensile  Strength  (3), 

Tensile  Modulus  (3), 

Tensile  Strain  at  Failure  (2), 
Minor  Poisson's  Ratio  (3), 

Compressive  Strength  (3), 
Compressive  Modulus  (2), 
Compressive  Strain  at  Failure  (2), 


an  = 6.23  x 103  PSI 
En  = 1.31  x 106  PSI 
eJt  = 0.50% 

Vy^  = 0.0136 

c'Tc  = 34.03  x 103  PSI 
EJc  = 1.45  x 106  PSI 
^ * 2.67% 


Tensile  Strength  (3), 

Tensile  Modulus  (3), 

Poisson's  Ratio  (3), 

Tensile  Strain  at  Failure  (2), 


a45t  * 11,94  x 103  PSI 
E45t  = 1.90  x 106  PSI 

v45  = 0.263 
e45  = °-70X 


Inplane  Shear  Strength  (3), 
Inplane  Shear  Modulus  (3), 


tlt  = 16.87  x 10J  PSI 
Glj  = 0.754  x 106  PSI 


* Calculated  from  data  on  tensile  specimens  with  45°  fiber  orientation. 
**  Denotes  number  of  tests  on  which  average  is  based. 
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Table  XXVIII 

SPECIMENS  FOR  INVESTIGATING  THE  TEST  VARIABLES 


VARIABLE 


TEST  SPECIMENS 
FOR 

INVESTIGATING  THE  VARIABLE 


FIBER  PROPERTIES 
MATRIX  PROPERTIES 
FIBER  LAYUP 
STACKING  SEQUENCE 
THICKNESS 
HYBRIDIZATION 
CURVATURE 


1.1a,  1.1b,  1.1c 
1.2a,  1.2b,  1.1b 
1.3a,  1.3b,  1.3c,  1.1b 
l.Ua,  l.Ub,  1.1b 
2.1a,  1.1b,  2.1c 
1.5a,  1.5b,  1.5c,  1.5d 
2.2a,  2.2b,  1.1b 


?7g>,  -'  >j*"» -*« j <£  *■>  C7  ■*  *v‘ 


Table  XXIX 

DESCRIPTION  OF  COMPOSITE  TEST  SPECIMENS 


DENSI TY  OF  MATERIAL 


FIBER 
CONTENT 
BY  VOLUME 


AVERAGE 

VOID  THICKNESS 
CONTENT  PER  PLY 


PLATE  SIZE  , CONTENT  VOID 

(WIDTH  Y LENGTH  X THICKNESS)  PLATE  VEIOiT  LB  /IN.  LB  /IN.  BY  VOLUME  CONTENT 
PLATE [INj (LB  ) UJ. [2] (<)  (2)  (»)  (2) 

1.1(a)  1 10.016  x 10.017  x 0.132L  0.809  0.0608 

1.1(a)  2 9.985  X 9.989  X 0.1317  0.801  0.0611  0.0611  61.7  0.12 

1.1(a)  3 9-986  x 10.02L  x 0.0981  0.601  0.0611 

1.1(B)  1 10.000  x 10.000  x 0.1391  0.772  0.0555 

1.1(b)  2 10.000  x 10.010  x 0.1395  0.772  0.0552  0.0556  58.0  0-93 

1.1(b)  3 10.015  x 10.015  x 0.1391  0.772  0.0553 

1.1(c)  1 10.008  x 10. OIL  x 0.1632  0.395  0.0SL6 

l.i(c)  2 iO.OlL  x 10.007  x 0.168’*  0.898  0.0531  0.05L1  55.9  0.81 

1.1(c)  3 10.00L  x 11.010  x 0.1690  0.896  0.0529 

1.2(a)  1 10.011  x It. 021  x O.lOLl  0.610  0.0582 

1.2(a)  2 10.016  X 10.021  x 0.1130  0 616  0.0568  0.0551  51.1  0.61 

1.2(a)  3 10.015  X 10.020  x 0.1116  0.628  0.0559 

1.2(b)  1 10.072  x 9.972  x 0.1229  0.623  0.0505 

1.2(b)  2 10.018  X 9.989  X 0.1211  0.618  0.0511  0.0510  50.7  8.50 

1.2(b)  i 10.022  x 9.981  x 0.1228  0.622  0.0507 

1.3(a'  1 12.017  X 12.012  x 0.1366  1.100  0.0557 

1.30.)  2 12.011  X 12.011  X 0.1397  1.12o  0.0555  0.0557  58.5  O.85 

1.3' a)  3 12.036  X 12.033  x 0.1103  1.120  0.0551 

1.3(b)  1 10.015  x 10.008  x 0.1L10  0.775  0.0518 

1.3(b)  2 10.012x10.012x0.1111  0.791  0.0517  0.0552  5o.2  0.92 

4.3(b)  3 10.016  x 10.016  x 0.1130  0.790  0.0550 

1.3(c)  1 9.996  x 10.010  x 0.1370  0.765  0.0512 

1.3(e)  2 10.005x10.002x0.1127  0.790  0.0551  0.0558  58.1  0.70 

1.3(e)  3 10.017  X 10.012  x 0.1119  0.786  0.0552 

1.1(a)  1 10.007  x 10.018  x 0.1321  0.737  0.0555 

1.1(a)  2 10.015  x 10.016  x 0.1332  0.739  0.0553  0.0561  60. b 0.81 

1.1(a)  3 10.020  x 10.007  x 0.1322  0.737  0.0555 

1.1(b)  1 10.013  x 10.0i5  x 0.1316  0.751  0.0555 

1.1(b)  2 10.015  X 10.012  X 0.1352  0.752  0.0551  0.0559  59.1  0.65 

1.1(b)  3 10.001  x 10.015  x 0.1351  0.751  0.0555 

1.5(a)  1 9.993  x 9.962  x O.llOO  0.885  O.063L 

1.4(a)  2 10.058  x 10.020  x 0.1391  0.900  0.06Ll  0.06L1  27.1/32.1  1.20 

1.5(a)  3 10.032  1 10.202  x 0.1352  0.910  0.0668  (3) 

1.5(b)  1 9.995  x 10.003  x 0.1389  0.892  0.0612 

1.5(b)  2 9.96o  x 10.002  x O.UlO  0.890  0.0612  0.0610  28.0/31.7  1.25 

1.5(b)  3 9-996  x 9.998  x 0.1391  0.690  0.0639  (3) 

1.5(c)  1 10.002  x 10.012  x 0.1510  0.789  0.0521 

1.5(c)  2 10.000  X 10.000  x 0.1511  0.79)  0.0513  0.0522  31.0/30.7  1.96 

1.5(c)  3 )o. 030  x 10.009  x 0.1560  v.,90  0.0501  (3) 

1.5(d)  1 9.995  x 10.005  x 0.1529  0.793  0.0519 

1.5(d)  2 10.005  x 10.027  x 0.1510  0.795  0.0525  0.0521  28.3/35.0  0.90 

1.5(d)  3 10.008  x 10.017  x 0.1531  0.799  0.0X20  (3) 

2.1(a)  1 10.820  X 10.256  x 0.066*.  0.101  0.05H 

2.1(a)  2 10.028  x 10.012  x 0.0662  0.369  0.0551  0.056'  63.5  1.03 

2.1(a)  3 10.012  x 10.012  x O.C651  0.363  0.0551 

2.1(c)  1 10.000  x 10.015  x 0.2675  1.521  0.0567 

2.1(c)  2 10.001  x 10.012  x 0.2707  1.527  0.0562  0.0557  5e.5  1.07 

2.1(c)  3 10.016  x 10.013  x 0.2685  1 .* 21  0.0565 

AP-1  9.961  x 10.185  x 0.1268  0 67  3 0.0521 

AP-2  10.075x10.200x0.2592  1.397  0.0525 

2.2(a)  1 8.076"  O.D.  X 3.985  x 0.138  0.283  0.0557 

2.2(a)  2 3.072"  O.D.  x 1.007  x 0.113  0.277  0.0526 

2.2(a)  3 3.076"  O.D.  x 1.000  x 0.138  0.280  0.0550 

2.2(b)  1 I Speclben.  contained 

2.2(b)  2 ? surface  irregularities 

2.2(b)  3 ( and  were  not  tested 

(1)  Froo  weights  and  diornsions  of  plates  or  cylinders. 

(2)  Froa  gr.vinrtrlc  test  and  acid  digestlo*  . 

(3)  Fir3t  end  second  r.uabcrs  denote  percentages  of  graphite  and  E-Glass  or  Kevlar-19  fibers  respectively. 
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SECTION  V 

EXPERIMENTAL  STUDIES  ON  IMPACT  RESPONSE  OF  COMPOSITES 


To  establish  threshold  energy  levels  required  to  cause  brittle  fracture  in 
graphite -epoxy  composites  and  to  establish  the  influence  of  variables 
discussed  previously  on  the  impact  response,  ball-drop  impact  tests  were 
conducted  on  the  various  composite  panels  and  cylinders.  Specimens  were 
instrumented  with  strain  gages  and  manganin  gages  to  enable  he  impact- 
induced  strains  and  surface  pressure  to  be  experimentally  determined. 

A number  of  techniques  were  employed  to  detect  damage  in  composite 
plates,  including  visual  observations,  dye  penetrants.  X-rays,  through- 
transmission,  and  coating  of  specimens  with  powder  prior  to  impact.  Use 
of  surface  coatings  and  dye  penetrants  to  detect  surface  damage  was  fairly 
successful.  Internal  damage  within  composite  specimens  could  not  be 
detected  by  conventional  techniques. 

TEST  TECHNIQUE  AND  SPECIMEN  INSTRUMENTATION 

Because  of  the  low  predicted  velocities  required  to  initiate  damage 
(v  = 30-50  in.  /sec),  the  impact  tests  were  conducted  using  ball-drop  tests 
(employing  a 1.5-in.  -diameter  steel  sphere)  rather  than  low-pressure 
gas -gun  tests  as  planned  initially.  The  test  setup  is  shown  schematically 
in  Figure  59,  whereas  Figure  60  shows  the  actual  test  setup,  each  test 
setups  have  previously  been  employed  at  MDAC.  Impact  velicities  ranging 
up  to  v = 300  in.  /sec  were  achieved.  Moreover,  ball-drop  impact  tests 
have  also  been  used  by  other  investigators  studying  the  impact  response  of 
graphite -epoxy  composites. 

The  flat  plate  specimens  and  the  cylindrical  specimens  were  supported  as 
shown  in  Figure  61.  Several  plate  specimens  were  also  tested  using 
6-in.  -diameter  support  rings.  Two  of  the  completed  cylindrical  specimens 
ready  for  testing  are  shown  in  Figure  62. 

Two  types  of  instrumentation  were  employed  on  several  of  the  plate  and 
cylindrical  specimens:  (1)  strain  gages  to  measure  strain  versus  time 
and  (2)  manganin  gages  to  measure  pressure  versus  time.  Both  types  of 
instrumentation  were  used  on  the  flat  plates  and  on  the  cylinders.  The 
instrumentation  for  a given  impact  location  in  plates  was  as  follows: 

1.  One  manganin  gage  at  the  impact  point  (r  = 0). 

2,  One  biaxial  strain  gage  on  the  surface  opposite  the  impact 
surface  (r  = 0). 
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Figure  59.  Test  Setup  for  Low- Velocity  Impact. 
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Figure  60.  Test  Setup  and  Instrumentation. 


t (VARIABLE) 


Figure  62.  Cylindrical  Test  Specimens. 
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3.  One  strain  gage  midway  between  the  impact  site  and  the  edge  of 
the  support  (r  = R/2). 

4.  One  strain  gage  at  a quarter  point  from  the  impact  site  to  the 
edge  of  the  support  (r  = R/4). 

where  R is  the  radius  of  the  support  rings  (R  = 1.5  in.  ).  The  instrumen- 
tation for  a given  cylinder  was  as  follows: 

1.  One  manganin  gage  at  the  impact  point. 

2.  One  biaxial  gage  on  the  surface  opposite  the  impact  surface. 

3.  One  strain  gage  at  180°  from  the  impact  point. 

Table  XXX  describes  in  detail  the  instrumentation  that  was  used  in  the 
various  specimens. 


The  strain-versus-time  traces  for  a given  impact  test  were  obtained 
using  an  oscilloscope.  The  technique  for  recording  the  output  from  the 
instrumented  specimens  was  verified  using  a 1/8-in.  -thick  aluminum 
plate.  The  instrumentation  on  the  latter  plate  was  the  same  as  on  the 
all-composite  plates.  Typical  outputs  of  the  various  gages  are  shown  in 
the  oscilloscope  traces  of  Figure  63.  Figure  63(a)  shows  two  oscilloscope 
traces.  The  top  trace  is  for  the  manganin  pressure  gage  located  at 
r = 0;  the  vertical  sensitivity  is  1 x 10"^  V /cm.  The  bottom  trace  is  for 
element  No.  1 of  the  biaxial  strain  gage  located  at  r = 0.  The  vertical 
sensitivity  is  1000  pin.  /in.  /cm.  Figure  63(b)  shows  three  oscilloscope 
traces.  The  top  trace  is  for  element  No.  2 of  the  biaxial  strain  gage 
located  at  r = 0;  the  vertical  sensitivity  is  1000  pin.  /in.  /cm.  The  middle 
trace  is  for  the  strain  gage  at  R/4;  the  vertical  sensitivity  is  500  pin.  /in.  / 
cm.  The  bottom  trace  is  for  the  strain  gage  at  R/2;  the  vertical  sensi- 
tivity is  500  pin.  /in.  /cm.  The  horizontal  sensitivities  were  200  psec/cm 
for  all  traces.  As  can  be  seen,  for  a 1.  5- in. -diameter  steel  sphere 
impacting  at  100  in,  /sec,  the  duration  of  impact  was  approximately  1 mil- 
lisecond, with  a peak  strain  reading  of  approximately  3300  pin. 


Those  plates  and  cylinders  which  were  not  instrumented  were  tested  only 
under  impact  using  1.  5-in.  -diameter  steel  sphere  as  the  impactor.  For 
a given  plate  the  distance  between  impact  sites  was  ~3  inches.  The 
techniques  that  were  employed  for  detecting  impact  damage  are  described 
later  on. 


RESULTS  FROM  NONINSTRUM  ENT  ED  IMPACT  TESTS 


The  test  setup  described  in  the  previous  paragraphs  was  used  to  study  the 
influence  that  the  variables  described  in  Table  XXVIII  had  on  the  impact 
response  of  composites.  The  impact  velocities  selected  on  the  basis  of 
theorectical  studies  were  30  in. /sec,  50  in. /sec,  and  100  in. /sec.  Initial 
impact  tests  performed  on  the  various  specimens  at  the  above  velocities 
did  not  show  any  visible  damage.  Moreover,  even  at  150-in.  /sec  impacts, 
no  damage  could  be  observed.  X-rays  were  then  taken  of  the  various 


Table  XXX 

SPECIMEN  INSTRUMENTATION* 


Plate  or 
Cylinder 


1.1(b)-! 


l.l(b)~2 


1.1(c)-! 


l.l(c)-2 


1.3(c)-! 


2.2(a) 


2.2(b) 


Location 

On  Plate  Or  Front  Face 
Cylindev  ( Impacted) 


Back 

Face 

r=0 

r=R/k 
(or  0=90) 

r-R/2 

(or  6=l80°) 

BSG 

USG 

USG 

BSG 

USG 

USG 

BSG 

1 

SG 

USG 

BSG 

USG 

USG 

BSG 

USG 

USG 

BSG 

USG 

- 

USG 

BSG 

USG 

USG 

- 

- 

- 

(l)  MG  denotes  manganin  gage 
BSG  denotes  biaxial  gage 
USG  denotes  uniaxial  gage 


(2)  Gage  was  on  outer  surface  of  cylinder 
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Figure  63.  Oscilloscope  Traces  From  Instrumented  Aluminum  Plate 
Impacted  at  100  In.  /Sec  With  a 1.  5 -In.  -Diameter  Steel 
Sphere. 
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specimens  to  see  if  impact  damage,  such  as  cracking  or  delaminations, 
could  be  detected.  The  latter  technique  proved  unsuccessful.  Another 
technique  that  was  tried  in  an  attempt  to  establish  at  w at  impact  velo- 
cities damage  initiated,  consisted  of  sprinkling  white  power  on  the  front- 
faces  (impacted)  and  backfaces  of  the  specimens  before  impact.  The 
technique  was  quite  successful.  Upon  impact,  a "footprint"  was  left  on 
the  impacted  surface,  showing  an  area  proportional  and  concentric  to  the 
area  of  contact.  Moreover,  when  failure  such  as  transverse  cracking 
took  place  on  the  backface  of  the  specimens,  this  could  readily  be  seen. 

All  subsequent  testing  of  graphite -epoxy  specimens  involved  use  of  powder 
coating.  The  specimens  that  were  tested  initially  (without  any  prepara- 
tions for  determining  damage)  were  retested  by  impacting  them  at  another 
location  on  the  plate.  Failures  were  observed  to  occur  at  impact  veloci- 
ties of  100  in.  /sec,  whereas  during  the  initial  testing  of  unprepared  speci- 
mens no  damage  could  be  detected  visually  or  by  using  x-rays,  even  at 
velocities  of  150  in.  /sec. 

Impact  tests  performed  on  a number  of  specimens  at  velocities  of  30  in.  / 
sec  and  50  in. /sec  showed  no  surface  damage;  consequently,  the  majority 
of  the  specimens  were  subjected  to  impact  velocities  equal  to  or  greater 
than  100  in.  /sec.  Testing  at  lower  velocities  was  conducted  only  in  the 
case  of  specimens  exhibiting  significant  damage  at  100  in.  /sec.  It  was 
found  that  a number  of  specimens  did  indeed  fail  at  impact  velocities  of 
50  in.  / sec. 

Influence  of  Fiber  Properties  on  Impact  Response 

The  specimens  used  in  establishing  the  influence  of  fiber  properties  on 
impact  response  are  designated  as  1.  la,  1.  lb,  and  1.  lc.  The  fiber  layup 
in  all  three  specimens  was  pseudo-isotropic,  whereas  the  matrix  was 
Na-mco  5208  epoxy.  The  reinforcing  fibers  in  the  three  specimens  were 
Celior.  GY-70  in  1.  la;  Thornel  300  in  1.  lb,  and  Modmor  II  in  1.  lc.  Results 
from  the  impact  tests  performed  at  velocities  ranging  from  50  in.  /sec  to 
160  in.  /sec  are  presented  in  Table  XXXI.  None  of  the  three  materials 
failed  at  an  impact  velocity  of  50  in.  /sec.  All  the  impacts  at  100  in.  /sec 
caused  significant  damage  in  Celion  GY-70/5208.  One  of  the  1.  lb  spec- 
imens impacted  at  100  in.  /sec  also  showed  cracking  in  the  transverse 
direction.  Only  one  out  of  three  of  the  1.  lc  specimens  failed  under 
100 -in.  /sec  impact.  Upon  increasing  the  impact  velocity  to  128  in.  / 
sec,  the  latter  specimen  (1.  lc)  did  not  show  any  damage.  Further 
increase  in  velocity  to  160  in.  /sec  caused  cracking  on  the  backface  of 
the  plate.  Figure  64  shows  typical  failures  resulting  from  impact  on  the 
three  groups  of  specimens.  In  all  cases  cracking  occurred  in  the  trans- 
verse direction.  In  the  case  of  Celion  GY-70/5208,  the  failure  was  at 
least  three  layers  deep,  as  can  readily  be  seen  from  Figure  64..  where 
the  failure  pattern  shown  follows  the  fiber  pattern  in  adjacent  layers. 
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Of  the  three  materials  for  which  the  results  are  shown  in  Table  XXXI  and 
in  Figure  64,  the  rating  of  materials  with  respect  to  resistance  to  impact 
damage  is  as  follows  (starting  with  the  best  material): 


v = 100  IN. /SEC 
SPEC.  § 1.1b- 


v = 100  IN . /SEC 
SPEC.  # 1. lb- 


1,  la  - CELION  GY70/5208 
1.1b  - THORNEL  300/5208 
1.1c  - MODMOR  11/5208 


Figure  64.  Continued. 
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INFLUENCE  OF  FIBER  PROPERTIES  ON  EXPERIMENTALLY 
OBSERVED  IMPACT  DAMAGE 
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If  two  numbers  are  the  same  the  area  proportional  to  the  e.rea  of  contact  was  circular;  if  two  numbers 
are  different  the  area  was  elliptical. 

Numbers  in  brackets  denote  width  of  damage  zone. 


2.  Thornel  300/5208 

3.  Celion  GY-70/5208 

The  Modmor  11/5208  specimens  were  approximately  30%  thicker  than  the 
other  specimens  (see  Table  XXIX).  If  all  specimens  were  of  the  same 
thickness,  more  extensive  damage  would  be  expected  in  the  Modmor 
11/5208  specimens  than  is  shown  in  Table  XXXI  or  in  Figure  64. 

Influence  of  Matrix  Properties  on  Impact  Response 

Specimens  1.2a,  1.2b,  1.1b,  and  AP  of  Tables  XXVIII  and  XXIX  were  used 
for  establishing  the  influence  of  matrix  properties  on  impact  response. 

As  before,  the  fiber  pattern  in  Specimens  1.2a,  1.2b,  and  1.  lb  was 
pseudo-isotropic,  and  the  reinforcing  fibers  used  in  these  specimens 
were  Thornel  300.  Specimens  AP  were  added  to  the  test  group  at  the  sug- 
gestion of  the  Project  Officer.  These  specimens  were  also  made  with 
Thornel  300  fibers;  however,  the  layup  was  not  pseudo-isotropic  (see 
Table  XXVII).  The  results  from  the  impact  tests  are  shown  ’n  Table 
XXXII,  whereas  Figures  65,  66,  and  67  show  typical  failure  in  the  various 
specimens.  Of  the  four  groups  of  specimens,  those  made  with  a poly- 
sulfone  matrix  are  unquestionably  most  resistant  to  impact  damage.  The 
four  groups  of  materials  are  rated  as  follows  with  respect  to  resistance  to 
impact  damage: 

1.  Thornel  300/polysulfone 

2.  Thornel  300/5208 

3.  Thornel  300/4617 

4.  Thornel  300/APCO 

Influence  of  Fiber  Orientation  on  Impact  Response 

Composite  plates  (consisting  of  Thornel  300/5208)  with  four  different 
fiber  orientations  were  used  in  experimental  investigation  of  the  influence 
of  the  Latter  parameters  on  impact  response.  The  fiber  orientations  in 
the  various  specimens  were  as  follows: 

1.  Unidirectional  fibers  in  Specimen  i.3a 

2.  1:1  Bidirectional  fibers  in  Specimen  1.3b 

3.  2:1  Bidirectional  fibers  in  Specimen  1.3c 

4.  r ridirectional  fibers  in  Specimen  1.  lb 
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v = 100  IN. /SEC 
SPEC.  # APCO  - 


Figure  65.  Continued. 
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v = 100  IN./SEC  V = 192  IN./SEC 

SPEC.  # 1.2  b-  SPEC  # 1. 26- 


Figure  66.  Impact  Damage  in  Thornel  300-Polysulfone  as  a Function 
of  Impact  Velocity. 
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v = 250  IN. /SEC 
SPEC.  # 1.2b- 

OVERALL  VIEW  OF 
FRONT  FACE 


v = 250  IN. /SEC 
SPEC.  # 1.2b- 

ENLARGED  VIEW 
OF  IMPACT  SITE 


NOTE:  NO  DAMAGE  WAS  VISIBLE  ON  THE 
BACK  FACE  OF  THE  SPECIMEN 


Figure  67.  Damage  on  the  Top  Face  in  Thornel  300-Polysulfone 
Composite  Plate. 
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INFLUENCE  OF  RESIN  PROPERTIES  ON  EXPERIMENTALLY 
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The  pertinent  data  on  the  impact  damage  in  the  specimens  described  above 
are  presented  in  Table  XXXIII.  Figure  68  shows  some  typical  failure 
modes  in  the  various  specimens.  From  the  results  presented  in 
Table  XXXIII  and  Figure  68,  the  rating  of  the  four  groups  of  specimens 
for  resistance  to  impact  damage  is  as  follows  (starting  with  the  best 
material): 

1.  1:1  Bidirectional  Thornel  300/5208 

2.  Tridirectional  Thornel  300/5208 

3.  2:1  Bidirectional  Thornel  300/5208 

4.  Unidirectional  Thornel  300/5208 

As  is  readily  seen  from  Figure  68,  in  the  case  of  unidirectional  speci- 
mens, the  crack  formed  on  impact  progressed  through  the  width  of  the 
specimen,  and  only  the  support  ring  stopped  further  crack  growth.  In 
the  case  of  plates  with  multidirectional  fibers,  the  reinforcement  normal 
to  the  crack  direction  acted  as  the  crack  stopper  and  thus  limited  crack 
growth. 

Influence  of  Stacking  Sequence  on  Impact  Response 

The  test  results  for  the  influence  of  the  stacking  sequence  on  impact 
response  are  presented  in  Table  XXXIV,  which  also  contains  a description 
of  the  stacking  sequences  in  the  various  groups  of  specimens.  Typical 
results  for  the  failure  modes  within  groups  of  specimens  are  presented 
in  Figure  69.  Based  on  data  presented  in  Tables  XXXI  and  XXXIV,  the 
rating  of  specimens  with  respect  to  resistance  to  impact  damage  is  as 
follows: 


1.  Specimens  1.  lb  with  stacking  sequence  10, -60, +60] ^)g 

2.  Specimens  1.4a  with  stacking  sequence  {(+60, -60)^  (0)^]g 

3.  Specimens  1.4b  with  stacking  sequence  1(0)^  (+60, -60)^]g 
Influence  of  Plate  Thickness  on  Impact  Response 

Plates  of  three  different  thicknesses  were  used  to  experimentally  investi- 
gate the  influence  of  the  latter  parameter  on  impact  response.  The  nom- 
inal plate  thickensses  were  0.066  in.,  0.  14  in.,  and  0.27  in.  Most  of 
the  specimens  were  made  with  Thornel  300/5208,  with  a pseudo-isotropic 
fiber  layup,  l'0,  60, -60)nls,  where  n was  selected  to  give  the  desired 
thickness.  One  0.  26-in.  plate  made  with  Thornel  300/APCO  resin  was 
also  tested  for  impact  response  in  order  to  compare  the  impact  response 
data  with  the  corresponding  data  obtained  on  a thinner  plate  described 
previously  (see  Table  XXXII).  The  pertinent  results  from  the  impact 
tests  are  presented  in  Table  XXXV,  whereas  Figures  70,  71,  and  72 
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;D!  10°  IN./SEC 
SPEC.  # 1.3c- 


1.3a  - UNIDIRECTIONAL  THORNEL  300/5208 
1.3b  - 1:1  BI-DIRECTIONAL  THORNEL  300/5208 
1.3c  - 2:1  BI-DIRFCTIONAL  TH0PNEL  300/5208 


Figure  68.  Continued. 
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INFLUENCE  OF  PLATE  THICKNESS  ON  EXPERIMENTALLY 
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v = 200  IN. /SEC 
SPEC.  # 2.1c- 


Figure  71.  Continued 
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v = 250  IN. /SEC 
SPEC.  # APC0- 


Figure  72.  Impact  Damage  in  Thornei  30Q-APCO  Resin 
Composite  Plate  (Back  Face). 
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show  some  typical  failure  modes  in  plates  of  various  thicknesses.  The 
results  for  the  0,  14-in.  -thick  plate  were  presented  in  Table  XXXT  and 
Figure  64.  As  one  would  expect,  increasing  the  thickness  of  the  plates 
causes  increased  resistance  to  impact  damage. 

To  establish  if  impact-induced  damage  could  be  measured  independently 
of  visual  observations,  two  types  of  through-transmission  techniques  were 
used  on  the  0.  066 -in.  -thick  plate:  (1)  ultrasonic  attenuation  and  (2)  meas- 
urement of  longitudinal  velocity  through  the  thickness  of  the  plate.  The 
experimental  results  from  these  tests  are  presented  in  Figures  73  and 
74.  From  the  results  shown  in  Figure  73,  it  is  apparent  that  an  impact 
velocity  of  ^100  in.  / sec  causes  change  in  the  material — probably  delami- 
nation within  the  plate — since  the  transverse  cracking  was  observed  to 
take  place  at  an  impact  velocity  of  75  in. /sec  (see  Table  XXXV.  The 
step -type  curve  for  the  longitudinal-wave  velocity  versus  impact  velocity 
(Figure  74)  is  indicative  of  the  different  failure  modes  that  take  place  as 
the  target  is  impacted  at  progressively  higher  velocities.  Additional 
experimental  studies  are  required  in  order  to  understand  and  interpret 
the  results  from  the  two  test  techniques  and  to  accurately  assess  their 
value  in  the  impact  problems  involving  composites. 

Influence  of  Hybridization  (Combining  Different  Materials)  on  Impact 
Response 

To  establish  if  impact  resistance  of  graphite -epoxy  could  be  improved 
through  hybridization,  experimental  studies  were  conducted  on  the  following 
two  different  combinations  of  materials: 


1.  Bidirectional  Thornel  300/5208  combined  with  bidirectional 
fiberglass /5208. 

2.  Bidirectional  Thornel  300/5208  combined  with  bidirectional 
Kevlar  49/5208. 


As  described  in  Table  XXVII,  for  either  combination,  the  materials  were 
either  dispersed  uniformly  through  the  thickness,  or  the  second-phase 
material  was  added  on  the  top  and  bottom  surfaces  of  graphite-epoxy, 
giving  a three-layer  construction.  The  experimental  results  for  impact 
response  of  hybrid  composites  are  presented  in  Table  XXXVI.  Figures  75 
and  76  show  typical  failure  modes  in  the  Thornel  300-fiberglass-5208 
hybrid  composite  plates.  Impact  velocities  as  low  as  50  in.  / sec  were 
sufficient  to  cause  failures  in  the  bottom  layer.  At  an  impact  velocity  of 
100  in.  / sec,  failure  in  at  least  two  adjacent  layers  was  observed  to  take 
place  as  shown  by  the  orthogonal  crack  pattern  corresponding  to  the  fiber 
directions  in  the  two  layers.  To  enhance  the  crack  visibility  and  to  ensure 
that  no  cracks  existed  before  testing,  dye  penetrants  were  applied  to  the 
surface  opposite  the  impacted  surface.  White  powder  coatings  employed 
for  graphite -epoxy  plate  could  not  be  used  for  hybrid  composites  because 
there  was  insufficient  contrast  between  the  white  color  of  the  powder,  and 
the  yellow  color  of  the  fiberglass  and  Kevlar  49.  The  use  of  dye  penetrants 
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v = 100  IN. /SEC 
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Figure  76.  Continued. 
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(3)  D denotes  layers  of  bidirectional  fiberglass  or  Kevlar  49  dispersed  between  layers  of  bidirectional 
graphite;  3L  denotes  three-layer  construction  with  outer  two  layers  being  bidirectional  fiberglass 
or  Kevlar  49  and  center  layer  being  bidirectional  graphite. 


v = 100  IN. /SEC 
SPEC.  # 1. 53- 


Figure  75.  Impact  Damage  in  Thornel  300-Fiberglass-5208  Epoxy 
Hybrid  Composite  Having  Dispersed  Construction. 
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was  quite  successful,  as  can  be  seen  from  Figure  76.  In  the  case  of  hy- 
brid composite  plates  made  with  fiberglass,  no  cracks  were  detected  in 
the  specimens  before  impact  testing.  On  the  contrary,  hybrid  composite 
plates  made  with  Kevlar  49  exhibited  extensive  cracking,  transverse  to  the 
fiber  direction,  before  testing.  Figure  77  shows  the  crack  pattern  in  one  of 
the  plates.  The  surface  on  which  these  cracks  were  present  was  adjacent  to 
the  caul  plate  during  specimen  fabrication.  The  cracks  on  the  opposite  plate 
surface,  which  was  adjacent  to  the  bleeder  cloth  during  fabrication,  were  not 
as  numerous  as  those  shown  in  Figure  76.  Consequently,  when  testing  the 
plates  made  with  Kevlar  49-Thornel  300,  the  latter  surface  was  placed 
opposite  the  impact  surface.  Figures  78  and  79  show  some  typical  failure 
modes  in  hybrid  composites  made  with  Kevlar  49  and  Thornel  300.  Fig- 
ure 78a  shows  the  crack  pattern  in  the  plate  before  testing.  Although  use  of 
dye  penetrant  enhanced  the  visibility  of  the  cracks,  they  could  not  be  photo- 
graphed readily.  Consequently,  to  permit  the  initial  crack  pattern  to  be 
photographed,  ink  traces  were  then  marked  on  the  surface,  as  shown  in 
Figure  78b. 

An  approach  similar  to  that  described  above  was  used  in  identifying  which 
cracks  formed  during  impact  in  the  three-layer  Thornel  300-Kevlar 
49-5208  epoxy  composite.  Typical  results  for  the  failure  modes  in  the 
latter  material  are  presented  in  Figure  79. 

From  the  results  presented  in  Table  XXXVI  and  Figures  75  through  79, 
hybridization  does  not  appear  to  increase  the  threshold  velocity  at  which 
impact  induces  damage.  Part  of  the  reason  why  impact  response  of 
hybrid  composites  was  so  poor  may  be  the  presence  of  initial  cracks  in 
the  as  - fabricated  specimen.  Although  use  of  dye  penetrants  allowed 
identification  of  the  cracks  open  to  the  plate  surfaces,  there  is  also  a pos- 
sibility that  the  specimens  could  have  contained  subsurface  cracks  into 
which  the  dye  penetrant  could  not  penetrate.  The  results  presented  in 
Table  XXXVI  seem  to  indicate  that  for  resistance  to  impact  damage, 
dispersed  construction  is  superior  to  three-layer  construction,  and  hybrid 
composites  consisting  of  Thornel  300-Kevlar  49-5028  epoxy  are  some- 
what more  resistant  to  impact  damage  than  are  those  made  of  Thornel 
300-Kevlar  49- 5208  epoxy. 

Influence  of  Curvature  on  Impact  Response 

Three  types  of  specimens  were  used  in  studying  the  influence  of  curvature 
on  impact  response:  3-in.  - and  6-in. -diameter  cylinders,  and  flat  plates 
(R  =«).  The  fiber  layup  in  the  various  specimens  was  2:1  bidirectional. 

In  the  case  of  cylinders,  twice  as  many  fibers  were  oriented  in  the  hoop 
direction  as  in  the  axial  direction.  Most  of  the  impact  tests  were  con- 
ducted on  cylinders  whose  ends  were  potted  into  metal  rings,  as  shown  in 
Figure  61.  A 3-in.  -diameter  cylinder  with  free  ends  was  also  tested  to 
establish  how  the  end  fixity  influences  impact  response.  The  results  from 
impact  tests  on  cylinders  are  described  in  Table  XXXVII.  Impact  tests 
at  100  in.  /sec  on  a 3-in.  -diameter  cylinder  with  free  ends  showed  no 
damage  on  the  inner  surface  of  the  cylinder.  However,  when  a 3-in.  - 
diameter  cylinder  with  end  rings  (clamped  ends)  was  subjected  to  a 
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Figure  76,  Impact  Damage  in  Three- Layer  Thornel  300-Fiberglass- 
5208  Epoxy  Hybrid  Composite  Plate. 
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v = 100  IN. /SEC 
SPEC.  # 1.5b- 


Flgure  76.  Continued. 
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Figure  77.  Resin  Surface  Cracks  in  Thornel  300 -Kevlar  49  Hybrid 

Composite  Plate  Prior  to  Testing  (Cracks  are  in  Kevlar  49). 
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100-in. /sec  impact,  axial  cr  oks  were  observed  on  the  inner  surface  of 
the  cylinder.  In  the  case  of  cne  6-in.  -diameter  cylinder  pro'  ided  with 
metal  end  rings  on  both  ends,  impact  caused  circumferential  cracks  to 
form. 

The  difference  in  failure  modes  between  the  3-in.  - and  6-in.  -diameter 
cylinders  (axial  versus  circumferential  cracks)  was  caused  by  the  fact 
that  the  inner  layer  in  the  3 -in.  -diameter  cylinder  consisted  of  axial 
fibers,  whereas  the  inner  layer  in  the  case  of  the  6-in.  -diameter  cyl- 
inder consisted  of  hoop  fibers.  The  difference  in  construction  of  the 
two  cylinders  resulted  from  a fabrication  error.  The  shapes  of  the  areas 
of  contact  in  the  3-in.  - and  6- in.  -diameter  cylinders  are  shown  in  Fig- 
ure 80. 

Figure  81  shows  the  impact  damage  on  the  inner  surface  of  the  6-in.  -dia- 
meter cylinders.  Damage  on  the  inner  surface  of  the  3-in.  -diameter 
cylinders  could  not  be  photographed  readily  because  the  openings  in  the 
end  rings  were  small  (see  Figure  63).  Because  of  the  error  made  during 
the  fabrication  of  the  3-in.  -diameter  cylinders,  no  data  are  available  for 
a direct  comparison  of  the  impact  damage  versus  cylinder  diameter. 
However,  from  the  comparison  of  test  results  on  the  6-in.  -diameter 
cylinder  and  on  the  flat  plate  (see  Table  XXXV,  Specimen  1.  3c),  it  can 
be  seen  that  curvature  limits  the  damage  caused  by  impact.  Whereas  for 
a given  impact  velocity  the  average  crack  'ength  in  flat  plates  was  « 1 in., 
the  average  crack  length  in  the  6-in.  -diameter  cylinders  was  only 
«0.  37  in.  The  average  crack  length  in  the  3 -in.  -diameter  cylinders  would 
be  expected  to  be  even  shorter. 

Some  additional  results  on  the  influence  of  variables  discussed  in  this 
section  and  obtained  from  instrumented  impact  tests  are  presented  in  the 
following  paragraphs. 

RESULTS  FROM  INSTRUMENTED  IMPACT  TEbiS 

Attempts  to  measure  the  impact  response  using  specimens  instrumented 
with  strain  gages  and  manganin  gages  were  only  partially  successful.  The 
fact  that  the  areas  of  contact  resulting  from  impact  are  small  in  compari- 
son with  the  size  of  manganin  gages  created  difficulties  in  interpreting 
the  results.  To  overcome  this  problem  an  attempt  was  made  to  calibrate 
the  manganin  gages  by  applying  known  static  load  (through  the  impactor) 
and  measuring  the  resultant  strain.  This  turned  out  to  be  a reasonable 
solution;  however,  as  the  static  load  was  increased  another  problem  was 
encountered  - the  output  from  the  manganin  gage  changed  sign.  This 
occurred  because  manganin  gages  are  sensitive  both  to  surface  pressure 
normal  to  the  plane  of  the  gage  and  to  strains  in  the  direction  of  the  plane 
of  the  gage.  Thus,  w’hen  the  applied  static  load  was  low,  the  strains  due 
to  plate  bending  were  small  and  the  output  from  the  manganin  gage  was 
governed  by  applied  load.  As  the  surface  load  increased,  the  strains  due 
to  plate  bending  increased,  so  that  the  manganin  gage  output  was  governed 
by  both  the  normal  pressure  and  the  inplane  strains.  Therefore,  direct 
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measurements  of  the  surface  pressure  caused  by  impact  could  not  be  made 
without  undertaking  an  extensive  experimental  program. 

Difficulties  were  also  encountered,  -vith  the  strain-gage  instrumentation. 

The  problem  of  gage  size  compared  co  the  area  of  contact  was  minimized 
by  using  small  gages  (gage  length  of  «0.  62  in.).  However,  it  was  found 
that  the  gage  output  was  extremely  sensitive  to  the  impact  site.  Moreover, 
because  the  strain  gages  could  not  be  placed  directly  within  the  projected 
area  of  contact  resulting  from  impact,  the  experimental  results  show  only 
an  average  response  in  the  impact  region. 

Another  factor,  which  may  be  important  but  whose  effect  could  not  be 
assessed  readily,  is  the  influence  of  the  gage  itself.  Presence  of  the 
manganin  gages  at  the  impact  site  may  affect  the  impact  response;  whereas 
the  presence  of  the  strain  gages  at  the  backface  of  the  plate  could  influence 
the  failure  mode.  The  influence  of  both  types  of  instrumentation  could  not 
be  established  without  undertaking  a more  extensive  program  than  that 
discussed  in  Section  IV. 

Impact  Test  Data  on  Instrumented  Composite  Plates 

Typical  results  showing  oscilloscope  traces  of  strain  versus  time  versus 
velocity  are  shown  in  Figures  82  through  86.  Figure  82  shows  the  strain 
versus  time  in  Modmor  11/5208  pseudo- isotropic  plate  (Specimen  1.  lc) 
subjected  to  impact  velocities  of  30  in.  /sec,  50  in.  /sec,  and  100  in.  /sec. 
The  results  shown  there  are  on  the  backface  of  the  plate  (opposite  the 
impacted  surface)  at  r = 0,  that  is,  at  the  impact  point.  Figure  83  shows 
the  radial  strains  on  the  backface  of  the  plate  in  Specimen  1.1c  at  a 
distance  of  r = R/2  and  r - R/4  from  the  impact  point,  where  R is  the 
inner  radius  of  the  support  ring.  The  two  curves  shown  in  each  portion  of 
Figure  82  represent  strains  from  the  two  biaxial  gages.  The  differences 
in  strain  measured  by  the  tv.  ->  gages  are  due  to  the  sensitivity  of  the  gage 
location  with  respect  to  the  center  of  impact.  Comparison  of  strains  at 
r = 0,  r = R/4,  and  r = R/2  (Figures  82  and  83)  shows  that  the  radial 
strain  decreases  rapidly  with  increasing  distance  from  the  center  of  con- 
tact. At  r = R/4  the  radial  strain  is  14%,  while  at  r = R/2  the  radial  strain 
is  approximately  1.8%  of  the  strain  at  the  impact  point. 

Figure  84  shows  results  similar  to  those  shown  in  Figure  82  for  the  pseudo- 
isotropic plate  made  of  Thornel  300/5208  (Specimen  1.  lb).  The  influence 
of  repeated  impacts  of  constant  velocity  at  a given  point  on  the  strain-time 
response  in  Specimen  1.  lb  is  shown  in  Figure  85.  The  latter  figure  shows 
a flattening  of  the  strain-time  curve  with  increasing  impact,  which  is  caused 
by  decrease  in  modulus  of  the  plate  due  to  local  failure.  The  results  pre- 
sented in  Figures  84  and  85  are  for  two  different  specimens,  which  were 
made  identically  and  were  processed  at  the  same  time.  However,  com- 
parison of  strains  in  the  two  plates  shows  significant  differences.  From 
Figure  84b  the  maximum  strain  resulting  from  100-in.  / sec  impact  was 
3100  pin,  whereas  for  the  plate  for  which  the  results  re  shown  in  Fig- 
ure 85,  the  maximum  strain  on  the  first  impact  was  ^4000  pin.  The  dif- 
ference in  strains  is  attributed  to  the  gage  location  with  respect  to  the 
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v = 30  IN. /SEC 
h = 500  n IN. /IN. /CM 
1 « 200  (iSEC./CM 
r = 0 
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v = 50  IN. /SEC 
h - 1000  fiIN./IN./CM 
1 * 200  j»SEC. /CM 
r = 0 
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v = 100  IN. /SEC 

h = 1000  (jJN./IN./CM 

1 = 200  .uSEC./CM 
^ - n 


v = IMPACT  VELOCITY 
n = VERTICAL  SCALE 
I * HORIZONTAL  SCALE 
r = IMPACT  LOCATION 
(r  = 0 DENOTES  CENTER 
OF  IMPACT) 


MOTE:  THE  TWO  CURVES  SHOWN  IN  EACH 

PORTION  OF  THE  FIGURE  WERE 
OBTAINED  FROM  the  BIAXIAL 
GAGES  MOUNTED  ON  THE  SURFACE 
OPPOSITE  THE  IMPACT  SITE. 

THE  UPPER  CURVES  GIVE  STRAINS 
IN  THE  FIBER  DIRECTION;  LOWER 
CURVES  ARE  AT  90s. 


Figure  82. 


Impact- Induced  Strains  in  Pseudo- Isotropic  Modmor  11-5208 
Composite  Plate  (Specimen  1.1c)  at  r - 0. 
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v = 30  IN. /SEC 
h = 100  nIN./IN./Ctl 
I = 200  ^SEC./CM 
r = R/4  (TOP  CURVE) 
r = R/2  (BOTTOM  CURVE) 


v = SO  IN. /SEC 
h = 250  IN. /IN. /CM 
1 = 200  SEC. /CM 
r = R/4  (TOP  CURVE) 
r = P/2  (BOTTOM  CURVE) 


R/4  AND  R/2  LOCATIONS  ARE 
QUARTER  AND  HALF  THE 
DISTANCE  FROM  THE  IMPACT 
POINT  TO  EDGE  OF  SUPPORT 
RING  RESPECTIVELY. 


NOTE:  THE  GAGE  GRID  WAS 

NORMAL  TO  THE  FIBER 
DIRECTION.  SEE 
FIGURE  81  FOR 
ADDITIONAL 
INFORMATION. 


v = 100  IN. /SEC 
li  = 500  |iIN./IN./CM 
I = 200  hlSEC/CM 
r = R/4  (BOTTOM  CURVE) 
r = R/2  (TOP  CURVf) 


Figure  83.  Impact- Induced  Strains  in  Pseudo-Isotropic  Modmor  11-5208 
Composite  Plate  (Specimen  1.1c)  at  r = R/4  and  r - R/2. 
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v = 30  IN. /SEC 
h = 500  mIN./IN. /CM 
1 = 200  mSEC/CM 
r = 0 


v = 50  IN. /SEC 
h - 1000  (ilN./IN./CM 
1 = 200  hiSEC/CM 
r = 0 


1 1^- 

v = 100  IN. /SEC 
h * 2000  ,iIN./IN./CM 
1 = 200  mSEC/ CM 
r = 0 


Figure  84.  Impact- Induced  Strains  in  Pseudo-Isotropic  Thornel  300-5208 
Composite  Plate  (Specimen  1.1b)  at  r = 0. 
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IMPACT  n 

n = 1000  IN. /IN. /CM  (TOP  CURVE) 
h = 500  IN. /IN. /CM  (BOTTOM  CURVE) 


IMPACT  *2 

h = 1000  nIH./lN./CM  (TOP  CURVE) 
h = 500  |xIN./IN./CM  (BOTTOM  CURVE) 


R? 


IMPACT  =3 
h 2000  mIN./IN. /CM  (TOP  CURVE) 


IMPACT  ft 5 
h --  2000  |iIN./IN./CM  (TOP  CURVE) 


h = 1000  t,H./IN./CM  (BOTTOM  CURVE  n = 2500  H IN./IN./CM  (BOTTOM  CURVE) 

NOTE:  IN  ALL  CASES  v = 100  IN. /SEC;  1 = 200  pSEC/CM;  r = 0 

Figure  85.  Influence  of  Repeated  Impacts  on  Strain- Time  Response 
(Thornel  300-5203  Pseudo-Isotropic  Laminate 
Designated  1.  lb). 
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h = 2000  nIN./IN./Cfl 

1 = 200  fiSEC/CM 
r = 0 

NOTE:  TOP  CURVE  IS  FOR 

DIRECTION  1;  BOTTOM 
CURVE  IS  FOR 
DIRECTION  2.  VOLUME 
FRACTION  OF  FIBERS 
IN  DIRECTION  2 WAS 
TWICE  AS  HIGH  AS  IN 
DIRECTION  1. 


h = 500  nIN./IN./CM  FOR 
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impact  point  and  to  the  influence  of  the  manganin  gage.  The  specimen  for 
which  data  are  presented  in  Figure  85  contained  a manganin  gage  at  the 
impact  point,  while  the  specimen  for  which  data  are  presented  in  Fig- 
ure 84  did  not.  The  impact  test  results  on  an  instrumented  composite 
plate  having  2:1  bidirectional  layup  at  Thornel  300/5208  are  presented 
in  Figure  86.  The  test  data  obtained  from  impact  tests  on  instrumented 
plates  are  presented  in  Table  XXXVIII. 

Impact  Tests  on  Instrumented  Composite  Cylinders 

Tests  similar  to  those  described  in  the  preceding  section  were  conducted 
on  the  3-in.  - and  6-in.  -diameter  cylinders.  The  6-in.  -diameter  cylinder 
was  instrumented  with  a manganin  gage  at  the  impact  site  on  the  outer 
surface  and  biaxial  strain  gages  on  the  inner  surface.  Because  of  the 
problems  discussed  at  the  beginning  of  this  section,  direct  measurement 
of  the  surface  pressure  resulting  from  impact  could  not  be  made.  More- 
over, the  strain  gages  located  at  180°  from  the  impact  site  did  not  show 
any  output  because  of  the  lew  stresses  and  strains  at  that  location.  Con- 
sequently, the  only  strains  measured  in  the  cylindrical  specimens  were 
the  strains  at  the  impact  site  on  the  inner  surface  of  the  cylinder.  Typical 
results  for  the  axial  and  circumferential  strains  in  the  6-in.  - and  3-in.  - 
diameter  cylinders  are  presented  in  Figures  87  and  88.  Comparison  of 
the  ratio  of  axial  to  circumferential  strains  in  the  two  cylinders  showed 
this  ratio  to  be  significantly  different.  In  trying  to  establish  the  reason 
for  this  difference,  it  was  discovered  that  an  error  had  been  in  the  fabri- 
cation of  cylinders.  The  construction  of  the  two  cylinders  was  as  follows 
(see  Table  XVII): 

1.  Construction  in  3-in. -diameter  cylinder:  [0,  90,  90]g 

2.  Construction  in  6-in.  -diameter  cylinder:  [90,  90,0)^]g 

Because  of  this  error,  the  influence  of  the  cylinder  radius  on  the  impact 
response  could  not  be  assessed.  However,  the  influence  of  curvature  can 
be  assessed  by  comparing  the  data  on  a 2:1  bidirectional  plate  (see  Fig- 
ure 86)  with  the  data  on  the  6-in.  -diameter  cylinder.  For  an  impact 
velocity  of  100  in.  /sec,  comparison  of  results  shows  the  following: 


1.  Maximum  strain  in  the  cylinder  was  6800  pin. /in.  , whereas  in 
the  flat  plate  the  maximum  strain  was  6500  pin.  /in. 

2.  The  ratio  of  the  axial  to  hoop  strain  in  the  cylinder  was  2.26, 
whereas  the  ratio  of  maximum  to  minimum  strain  in  the  plate 
was  1.  55. 


3.  The  contact  duration  in  flat  plate  was  longer  than  in  the  cylinder 
(1400  psec  vs  1200  psec). 

4.  The  area  of  contact  resulting  from  impact  was  elliptical  in  the 
case  of  the  cylinder  and  nearly  circular  in  the  case  of  the  flat 
plate. 
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From  the  above  results  and  the  results  presented  previously  in  this 
section,  it  appears  that  curvature  has  an  influence  on  impact  response. 
Moreover,  it  would  be  expected  that  the  curvature  effects  would  increase 
with  decreasing  cylinder  radius. 
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Figure  87.  Impact- Induced  Strains  in  6- In.  - Diameter,  2:1  Bidirectional 
Thome!  300-5208  Composite  Cylinder. 
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h = 4000  jiIN./IN./CM 
1 = 200  pIN./IN./CM 
v = 100  IN. /SEC 
r = 0 


NOTE:  MATERIAL  LAYUP  IN 
THE  3- INCH -DIAMETER 
CYLINDER  WAS  DIFFERENT 
THAN  IN  THE  6-INCH- 
DIAMETER  CYLINDER. 


Figure  88.  Impact -Induced  Strains  in  3-In.  -Diameter 

3:1  Bidirectional  Thornel  300-5E08  Composite 
Cylinder. 
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SECTION  VI 

EVALUATION  OF  RESULTS,  CONCLUSIONS. 
A NDRECOMM  ENDATIONS 


Inasmuch  as  the  analytical  studies  on  impact  response  of  composite  plates 
were  conducted  prior  to  fabrication  and  testing  of  the  various  plates  and, 
as  described  in  Section  IV,  the  thicknesses  of  the  fabricated  plates  dif- 
fered from  the  thicknesses  used  in  the  analysis,  a direct  test-theory 
comparison  of  the  threshold  impact  velocities  to  initiate  damage  could 
not  be  made  for  all  the  plates  without  making  numerous  computer  runs 
for  the  various  cases.  Only  in  certain  cases  could  a test-theory  com- 
parison be  made  without  making  additional  computer  runs.  The  results 
of  the  latter  comparison  are  used  to  verify  the  theory.  Further  theory 
verification  is  obtained  by  comparing  the  theoretically  predicted  and 
expel  imentally  observed  trends  in  impact  response  of  composite  plates 
made  with  different  fibers  and  resins  and  having  different  fiber  orienta- 
tions and  layups,  different  thicknesses,  and  different  geometries. 

TEST-THEORY  COMPARISON  OF  THE  THRESHOLD  IMPACT  VELOCITY 
TO  INITIATE  DAMAGE 

As  noted  above,  the  plate  thicknesses  used  in  the  analytical  studies  were 
0.  125  in.  , whereas,  because  of  the  fabrication  constraints  discussed  in 
Section  IV,  the  actual  thicknesses  of  the  fabricated  specimens  were  some- 
what higher  (see  Table  XXIX).  Figure  36  shows  that  even  small  varia- 
tions in  the  plate  thickness  have  a significant  influence  on  the  surface 
pressure,  q*,  and  thereby  on  the  impact  velocity  at  which  failure  initiates 
at  a point  in  a plate.  Moreover,  as  shown  in  Figure  37,  changing  the 
radius  of  the  support  ring  does  not  have  a significant  effect  on  q*.  The 
analytical  results  shown  in  Figures  36  and  37  are  for  pseudo-isotropic 
composite  plates  made  with  Thornel  300  fibers.  Additional  analytical 
results  on  impact  damage  in  this  material  are  presented  in  Figures  25, 

26,  27,  45,  and  46.  With  the  information  presented  in  those  figures, 
it  is  possible  to  estimate  the  threshold  impact  velocity  to  cause  damage 
in  the  tridirectional  Thornel  300/5208  composite  platep.  For  a 0.  125-in. - 
thick  plate  it  was  shown  in  Figure  46  that  significant  damage  occurred 
on  the  bottom  surface  of  the  plate  when  the  plate  was  impacted  a: 

50  in.  /sec.  From  the  curve  labeled  No.  6 in  Figure  25,  the  maximum 
surface  pressure  corresponding  to  this  impact  condition  is  q = 102  ksi. 

The  average  thickness  of  the  plate  that  was  tested  under  impact  was 
t = 0.  1392  in.  (see  data  for  Specimen  1.  lb  in  Table  XXIX).  The  thres- 
hold impact  velocity  to  cause  significant  damage  in  the  0.  1392-in.  plate 
can  be  estimated  from  the  results  shown  in  Figures  36,  2 5,  and  89. 

Figure  89  is  a replot  of  the  portion  of  the  curve  snown  in  Figure  36, 
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igure  8°.  Influence  of  Plate  Thickness  on  Impact- Induced  Pressure 
Which  Initiates  Bending  Failure  at  the  Back  Face  of  the 
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which  gives  failure  on  the  bottom  surface  of  the  plate  due  to  bending 
stresses.  From  Figure  89,  the  value  of  nondimensional  thickness,  H, 
corresponding  to  t - 0.  125  in.  , v - 50  in.  /sec,  qQ  = 102  ksi  is  H = 2.  10. 
Taking  the  ratio  of  the  actual  plate  thickness,  t'  = 0.  1392  in. , to  the 
plate  thickness  used  in  the  analysis,  t = 0.  125  in.  , and  multiplying  the 
above-given  value  of  H by  the  thickness  ratio  gives  the  value  of  H'  for 
the  thicker  plate 


The  maximum  surface  pressure  to  cause  damage  for  the  case  when 
H'  = 2.34  is  (from  Figure  89) 


q*  = 128  ksi 

whereas  from  Figure  25  the  threshold  velocity  to  initiate  significant 
damage  is 


v = 88  io/sec 

The  extent  of  damage  on  the  bottom  surface  of  the  plate  can  be  estimated 
in  a similar  manner  using  the  results  presented  in  Figures  27  and  46.  The 
final  estimated  width  of  the  damaged  region  on  the  bottom  surface  of  the 
plates  is  w = 0.  224  in.  The  experimental  results  presented  in  Table  XXXI 
j for  Specimen  1.  lb  show  that,  out  of  four  impact  tests  conducted  at 

1 100  in.  /sec,  failure  was  observed  in  one  of  the  tests.  Moreover,  when 

the  impact  velocity  was  increased  to  128  in.  /sec,  two  of  the  specimens 
did  not  show  any  damage,  whereas  the  third  specimen  showed  extensive 
damage  at  that  impact  velocity. 

t 

| Although  no  results  such  as  shown  in  Figures  25,  26,  27,  36,  37,  45,  and 

j 46  were  generated  cor  other  plate  specimens,  nevertheless  the  results 

• presented  in  those  figures  can  be  used  to  obtain  a rough  estimate  of  the 

j impact  velocities  required  to  cause  damage  in  Specimens  1.  lc  and  1.  3a. 

| Applicability  of  the  above  results  to  Specimen  1.  lc  is  partially  justified 

■ inasmuch  as  the  properties  of  Modmor  II  fibers  do  not  differ  significantly 

f from  the  Thornel  300  fibers  (see  Tables  III  and  IV).  In  the  case  of 

I Specimen  i,3a.  unidirectional  Thornel  300-5208  composite,  results 

I similar  to  those  presented  in  Figure  45  are  given  in  Figure  43;  conse- 

I quently.  Figure  89  is  used  only  to  correct  for  the  thickness  effect. 

j In  the  case  of  Modmor  11-5208  composite  plates  (Specimen  1.  lc),  the 

| average  plate  thickness  was  0.  1669  in.  Using  an  approach  similar  to 

| that  described  above  gives  the  following: 

! H'  - 2.  80 
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For  the  unidirectional  Thornel  300/5208  composite  (Specimen  1.3a)  the 
predicted  velocity  to  initiate  significant  damage  in  0.  125-in.  -thick  plate 
was  v -30  in.  /sec,  which  corresponds  to  q*  = 84  ksi  and  H = 1.92. 
Estimate  of  the  threshold  impact  velocity  in  the  actual  test  specimen 
(t  = 0.  139  in. ) gives 

H*  = 2.  14 

q*  = 108  ksi 
v = 58  in.  /sec 

Impact  velocities  at  which  damage  was  observed  in  Specimens  1.  lc  were 
100  in. /sec  and  160  in. /sec,  whereas  in  the  case  of  unidirectional 
Thornel  3 J0-5208  composites  the  impact  velocity  that  caused  damage  was 
50  in.  /sec  <V  < 100  in.  /sec,  the  damage  extending  through  the  width  of 
the  plate  for  the  100-in.  /sec  impact  velocity  (see  Tables  XXXI  and 
XXXIII). 

The  results  presented  above  show  fair  correlation  between  tests  and 
theory  if  one  recognizes  that  the  theoretical  results  discussed  here  were 
obtained  in  an  approximate  manner.  Moreover,  no  account  was  made  in 
the  theory  for  the  influence  of  the  plate  deformation  on  the  surface 
pressure  or  for  the  fact  that  the  material  properties  used  in  the  initial 
theoretical  analyses  may  have  been  different  from  the  actual  material 
properties  within  the  various  plate  specimens.  Use  of  the  properties 
data  from  actual  test  specimens  would  be  expected  to  improve  test-theory 
correlation.  The  properties  data  needed  is  in  the  direction  of  the  thick- 
ness of  the  specimen.  For  example,  the  Young's  modulus  in  the  thick- 
ness direction  of  the  plate  has  a significant  influence  on  the  impact- 
induced  surface  pressure,  which  in  turn  affects  the  internal  stresses, 
the  failure  mode,  and  the  extent  of  failure. 

As  to  the  influence  c f the  plate  flexibility,  it  too  would  influence  the 
magnitude  and  distribution  of  the  surface  pressure  as  well  as  the  contact 
duration  resulting  from  impact.  The  contact  duration  is  shown  to 
increase  with  decreasing  Ez,  the  Young's  modulus  in  the  direction  of 
plate  thickness  \see  Table  VIII  and  Figure  26).  The  plate  flexibility 
could  be  looked  upon  as  decreasing  the  Young's  modulus  of  the  material, 
Ez,  to  some  effective  modulus,  Ez,  wh-»re  E'z  < Ez.  Thus,  E~  of  the 
plate  would  be  expected  to  decrease  in  comparison  to  Ez,  as  the  flexural 
rigidity  of  the  plate  increased,  which  would  cause  an  increase  in  the 
contact  duration.  The  experimental  results  for  instrumented  specimens 
(see  Figures  8 2 through  86)  do  indicate  that  the  contact  durations  are 
significantly  higher  than  the  theoretically  predicted  values,  which 
assumed  the  plate  bending  effects  to  have  a negligible  influence.  Whereas 
for  an  impact  velocity  of  100  in.  / sec,  the  theoretically  predicted  con- 
tact durations  were  between  240  psec  and  300  psec,  the  experimentally 
measured  contact  durations  were  of  the  order  of  600  psec.  The  fact  that 
the  flexibility  of  the  specimen  influences  the  contact  duration  follows 
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from  the  comparison  of  contact  durations  in  flat  plates  and  cylinders. 

For  a composite  having  bidirectional  layup,  the  contact  duration  resulting 
from  impact  of  a 1.  5 -in.  -diameter  steel  sphere  into  a cylinder  was 
~600  psec,  whereas  for  the  same  impact  condition  into  flat  plate  the 
impact  duration  was  700  psec. 

The  results  predicted  analytically  and  observed  experimentally  are 
briefly  compared  in  sections  that  follow. 

Influence  of  Fiber  Properties  on  Impact  Response 

Both  test  and  theory  agree  on  the  influence  of  fiber  properties  on  impact 
response.  The  resistance  of  impact  damage  increases  with  increasing 
fiber  strength  and  decreasing  fiber  modulus.  Impact  damage  is  more 
sensitive  to  strength  than  to  modulus. 

Influence  of  Matrix  Properties  on  Impact  Response 

Beth  tests  and  theory  show  that  for  a given  reinforcing  fiber  and  a given 
fiber  layup,  the  resistance  of  the  composite  to  impact  damage  increases 
as  the  Young's  modulus  of  the  matrix  decreases  and  its  tensile  strength 
increases.  Composites  made  with  a polysulfone  matrix  were  shown  to 
be  more  resistant  to  impact  damage  than  those  made  with  5208  epoxy 
or  ERLA  4617.  The  fact  that  the  composites  made  with  polysulfone  were 
of  poorer  quality  than  the  other  composites  (higher  void  content  and 
lower  stiffness  and  strength)  may  have  altered  the  mode  of  failure. 


Influence  of  Target  Thickness  on  Impact  Response 

Both  tests  and  theory  show  that  by  increasing  the  target  thickness  the 
velocities  required  to  cause  damage  by  a given  impactor  also  increase. 
Theory  shows  that  for  thick  targets,  local  or  subsurface  damage  at  the 
impact  site  will  occur  first.  Moreover,  the  theory  predicts  when  the 
damage  will  be  local  and  when  it  will  be  overall.  No  experimental  work 
was  conducted  to  define  the  type  of  subsurface  and  local  damage  as  a 
function  of  target  thickness  and  impact  velocity. 

Influence  of  Fiber  Orientation  on  Impact  Response 


Both  theoretical  and  experimental  results  show  that  a composite  having 
a 1:1  bidirectional  layup  is  more  resistant  to  impeat  damage  than  a 
composite  with  a tridirectional  fiber  layup.  The  latter  in  turn  is  shown 
to  be  significantly  superior  to  a composite  with  unidirectional  fiber  lay- 
up, as  one  would  expect.  Both  theory  and  experiments  show  that  the  area 
of  contact  in  uni-directional  composites  is  elliptical. 
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Both  theory  and  experiment  show  that  for  a given  composite  and  fiber 
layup,  uniform  and  complete  dispersion  of  layers  through  the  specimen 
thickness  gives  a more  impact-damage- resistant  material  than  does  a 
construction  in  which  the  fibers  are  not  dispersed  but  rather  are  placed 
in  distinct  layers.  For  example,  for  a 2:1  bidirectional  layup  in  which 
twice  as  many  fibers  are  in  Direction  x as  in  Direction  y,  complete 
dispersion  of  the  x fibers  between  the  y fibers  would  give  a more  impact- 
damage- resistant  construction  than  would  a construction  in  which  all  the 
x fibers  are  placed  in  the  two  outer  layers  and  the  y fibers  are  placed  in 
one  central  layer. 

Influence  of  Hybridization  on  Impact  Response 

Neither  theoretical  results  nor  experimental  results  indicate  any  dra- 
matic improvements  in  resistance  to  impact  damage  due  to  hybridization 
(combining  different  composite  materials).  In  fact,  in  some  cases  the 
experimental  threshold  velocities  to  initiate  damage  in  hybrid  composites 
were  lower  than  the  threshold  velocities  required  to  initiate  damage  in 
all  graphite- epoxy  composites.  In  the  case  of  graphite-fiberglass- resin 
hybrid  composites,  theory  shows  that  a three-layer  construction  where- 
by the  bidirectional  fiberglass  layers  are  placed  on  both  sides  of  a 
bidirectional  graphite- epoxy  layer  is  more  impact  resistant  than  a con- 
struction in  which  the  bidirectional  fiberglass  and  bidirectional  graphite 
layers  are  interspersed.  The  experimental  results  were  not  sufficient 
to  permit  an  accurate  assessment  of  which  type  construction  is  more 
efficient.  Test  results  on  hybrid  composites  consisting  of  Thorr.el  300- 
fiberglass- resin  and  Thornel  300-Kevlar  49-resin  show  that  for  a given 
impact  velocity  the  backface  damage  in  the  latter  composite  is  less 
extensive  than  in  the  former  one.  This  appears  to  be  in  agreement  with 
theoretical  predictions. 

Influence  of  Curvature  on  Impact  Response 

Both  theory  and  experiments  show  that  the  area  of  contact  resulting  from 
an  impact  of  a sphere  into  a cylinder  will  be  elliptical,  with  the  major 
axis  of  the  ellipse  coinciding  with  the  cylinder  axis.  The  theoretically 
predicted  ratio  of  the  axial-to-circumferential  strains  at  the  backface 
of  the  impact  site  was  ~1.  5,  whereas  the  tests  show  this  ratio  to  be  «2.  3. 
(Additional  results  for  the  curvature  effects  are  discussed  in  Section  III). 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  general,  the  theoretically  predicted  trends  for  the  influence  of  fiber 
properties,  matrix  properties,  fiber  layup,  stacking  sequence,  thick- 
ness, and  curvature  effects  show  fair  agreement  with  the  test  results. 
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Further  theory  refinement  appears  desirable,  however.  Among  the 
important  items  that  need  to  be  done  to  refine  the  theory  are  the  follow- 
ing two: 


1.  Determination  of  how  target  flexibility  influences  the  surface 
pressure,  internal  stresses,  and  contact  duration. 

2.  Determination  of  the  influence  of  local  compressive  yielding 
under  the  impactor  on  the  impact- induced  surface  pressure. 


Studies  reuuired  to  make  a more  accurate  assessment  of  the  reliability 
of  the  theory  should  include: 


3.  Complete  characterization  of  materials,  including  determinr.- 
tion  of  inplane  and  through-the-thickness  properties  of 
materials. 

4.  Experimental  studies  on  internal  failure  modes  in  composite 
targets  as  influenced  by  impact  velocity  and  target  thickness. 

This  should  include  development  of  experimental  techniques 
for  detecting  surface  and  subsurface  damage. 

5.  Additional  experimental  studies  using  instrumented  specimens. 

In  order  to  establish  which  type  of  impact  damage  is  critical,  it  is  recom- 
mended that  theoretical  and  experimental  work  be  conducted  on  speci- 
mens impacted  at  various  velocities.  This  should  include  determination 
of  reduction  in  static  and  fatigue  strength  and  elastic  properties  as  a 
function  of  impact  velocity.  Studies  are  also  recommended  on  repair 
techniques  for  impact- damaged  composite  materials,  including  a 
definition  of  when  it  is  feasible  to  repair  and  when  it  is  necessary  to 
replace  a composite  structural  component  damaged  by  impact.  Finally, 
it  is  recommended  that  the  refined  theory  for  predicting  the  impact 
response  of  composites  be  applied  and  implemented  through  fabrication 
and  testing  of  composite  materials  or  combinations  of  materials  that 
are  sfcjown  to  be  resistant  to  impact  damage. 
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LIST  OF  SYMBOLS 


a 


Major  radius  of  the  area  of  contact,  in. 


A..,  A,?,  A??  Terms  which  are  a function  of  elastic  properties  of 
material  (see  Equation  17). 


u 


Minor  radius  of  the  area  of  contact,  in. 


Term  which  takes  into  account  curvature  effect  (see 
Equation  3). 

Diameter  of  a circular  plate,  in. 

Nondimensionai  diameter  of  a circular  plate  D = ^ 
Young's  modulus 

Allowable  strength  of  material,  psi 
Shear  modulus,  psi 
Thickness  of  target,  in. 

Nondimensionai  thickness  of  a circular  plate  (H  = 

Stiffness  parameter  for  the  impactor  and  the  target, 
respectively,  psi-* 


L 

in,  n,  s 


mj,  m^ 


Length  of  cylinder,  in. 

Parameters  which  are  functions  of  the  principal  radii 
of  curvature  of  the  impactor  and  the  target  (see 
Figure  2) 

Mass  of  the  impactor  and  the  target,  respectively, 
lb-  sec^  /in. 


N,  Nj  See  Equations  (8)  and  (9) 

q Surface  pressure  resulting  from  impact,  psi 
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LIST  OF  SYMBOLS  (Continued) 

Maximum  surface  pressure  resulting  from  impact,  psi 

Radial  coordinate  or  radius  of  plate 

Nondimensional  radius  of  the  circular  plate  (R  = — ) 

(or  radius  of  a cylinder,  in. ) 

Principal  radii  of  curvature  of  the  impactor 

Principal  radii  of  curvature  of  the  target 

Time,  sec 

Contact  duration,  sec 
Velocity  of  the  impactor,  in.  /sec 
Velocity  of  the  target,  in.  /sec 
Approach  velocity  (v  = v^  + v^),  in.  /sec 
Average  radius  of  the  damage  zone,  in. 

Half  width  of  the  damage  zone,  in. 

Depth  of  damage  zone,  in. 

Displacement  under  impactor,  in. 

Maximum  displacement  under  an  impactor,  in. 
Derivative  of  a with  respect  to  t 
Strain,  in.  /in. 

Poisson’s  ratio 

3 

Density  of  material,  lb/in. 

Normal  stress,  psi 
Shear  stress,  psi 
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LIST  OF  SYMBOLS  (Continued) 


SUBSCRIPTS 
x,  y,  z 

6,  r,  z 

L,  T,  Z 

t,  c 
1,  2,  3 


Denote  properties  in  rectangular  coordinates 

Denote  properties  in  cylindrical  or  polar  coordinates 

Denote  fiber,  transverse  and  through- the-thickness 
(transverse)  direction  in  filamentary  composite 

Denote  tensile  and  compressive  properties 

Properties  in  principal  coordinates  of  the  material. 

If  used  with  k,  m,  or  R subscripts,  1 and  2 refer  to 
the  impactor  and  the  arget,  respectively. 


